Comments (3)
comment from @RickMethot on 2016-04-20:
From Rick:
matrix Q_setup(1,Nfleet,1,5)
1 do power,
2 env-var,
3 extra sd,
4 devtype(<0=mirror, 0=float_nobiasadj 1=float_biasadj, 2=parm_nobiasadj, 3=rand, 4=randwalk, 5=set parameter equal to the scaling coefficient so can then have a prior)
5 Q offset
fundamentally, do power and devtype could be replaced by a link designator;
and the Q offset need not be mandatory; instead it is just one possible added parameter depending on the link function that is chosen
Proposaed New Q_setup
1: fleet index number - by specifying this we do not need input for all fleets, only those that have survey data
2: link type
3: 0/1 to select extra sd parameter
4: 0/1 for biasadj or not
5: 0/1/2 to float q: where 0 means not float (old option 2), 1 means float in phase 1 only, 2 means float continuously (old option 5)
So, the q parameter will always show up in the list of parameters as long parameter lines
As a long parameter, it can invoke additional parameters to implement: block/trend, env linkage, and mean-reverting random walk (invokes 2 parameters for the stddev of q and rho for autocorrelation)
The family of possible link functions needs flexibility because the survey obs could be a time series of an environmental index, which could have + and - values
and not be linear in relation to the model factor which could be abundance, dev vector, or a derived process parameter
Link types, where "x" is a model quantity as selected by Svy_units:
<=-1: mirror
1: proportional - y=qx
2: proportional with offset - y=q(x-a)
3: proportional with offset and power: y=q*(x-a)^c
4: logistic
QUESTIONS:
-
can q use any of the time-varying options and also float? I think it can if the time-varying component of q is separated from the component that can float or have a prior.
e.g. y = q1q2x
where q1 can float or have a prior
and q2 starts with a value of 1.0 in the first year and can be time-varying after that -
Is it feasible to allow the non-q parameters, e.g. a and c in list above or the env link parameter, to be time-varyin
from ss3-source-code.
comment from @RickMethot on 2016-04-22:
progress on conversion of Q parms
- test with Q as parm, Qlink as power fxn, and float Q
- still need to test with env effect on Q, and to test mirror
- random devs approach currently is disabled
- still need to update ss_write.tpl to output the new format
- the new approach will create a q for any fleet with survey obs, even with the float option
- still need to revise ss_readctl_330 to be able to read the new format
- then need to create block/trend approach and dev approach like that for mgparm and selparm
New format is:
//Q_setup for 3.30
// 1: link type
// 2: extra input for link, i.e. mirror fleet
// 3: 0/1 to select extra sd parameter
// 4: 0/1 for biasadj or not
// 5: 0/1 to float
but need to add another field with fleet id so that not all fleets will need to be read
the above info populates a fleet x 5 matrix with parameter index for each type of parameter needed
Q_setup_parms(1,Nfleet,1,5) // index of first parameter for: 1=base q; 2=env; 3=block/trend; 4=dev; 5=extrastd;
current link functions are:
// Link types
// 1 simple q, 1 parm
// 2 mirror simple q, 1 mirrored parameter
// 3 q and power, 2 parm
from ss3-source-code.
comment from @RickMethot on 2020-08-25:
consolidate and rename this issue
from ss3-source-code.
Related Issues (20)
- [Feature]: add discard body weights to wtatage.ss
- [Refactor]: create function for t=f(y,s)
- Extending time-blocks into forecast years HOT 1
- refactor Fmethod controls and expand Fmethod 4 capabilities
- [Feature]: 10+ multiyear average F_std_scaling; with change to log(ratio) approach HOT 3
- [Feature]: Revise recdev_bias adjustment to be a customizable list of line segments
- [Bug]: Possible bug relating to the estimation of Q_extra_SD when doing NUTS MCMC HOT 10
- [Feature]: Mean recruitment adjustment HOT 10
- [Feature]: Add 2-parm catchability (Q) option with offset and slope HOT 2
- [Feature]: add link to artifact of exe build to pull request discussion?
- [Bug]: "inf" value in bootstrap of mean size-at-age data for ages with sample size 0
- [Feature]: Cohort growth deviation added to Linf (or L2) HOT 1
- Add more detail for "catches" warning HOT 1
- [Feature]: Augment Mgmt_Quant report with recruitment values HOT 2
- [Feature]: One step ahead (OSA) residuals for composition data
- [Bug]: Time_Series output of forecast recruitment in a seasonal model HOT 1
- revise SPR_SERIES in Report.sso to clarify HOT 5
- revise RECR_DIST output to work better for models without 50/50 sex ratio at recruitment? HOT 2
- [Feature]: consider adding more stock-recruit functions HOT 1
- [Bug]: rho parameters not impacting semi-parametric (2D-AR) selectivity when multiple sigma values are used HOT 6
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from ss3-source-code.