thanks Laurie ... the problem now is this line:<br><br>print("OK to here3")<br>ctrl<-fwdControl(data.frame(year=yr.TAC,max=c(NA,2),val=c(mean(TAC),NA), quantity=c("catch","f")))<br><br>Year =2011<br>
[1] "OK to here1"<br>[1] "OK to here2"<br>=== mean TAC 32935.08 <br>[1] "OK to here3"<br>Error in `[<-`(`*tmp*`, 1, ac(min(ctrl@target[, "year"])), value = NA) : <br> error in evaluating the argument 'j' in selecting a method for function '[<-': Error in Summary.factor(c(1L, 1L), na.rm = FALSE) : <br>
min not meaningful for factors<br><br>Sorry for the inconvenience<br><br>Best,<br>Luis<br><br><br><div class="gmail_quote">On Wed, Mar 21, 2012 at 5:06 PM, Laurie <span dir="ltr"><<a href="mailto:lauriekell@googlemail.com" target="_blank">lauriekell@googlemail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Arial;word-spacing:0px"><span style="font-family:Monospace;font-size:13px;line-height:17px;text-align:-webkit-left;white-space:pre-wrap">
<pre style="font-family:Monospace;font-size:10pt!important;outline-style:none;outline-width:initial;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-width:initial;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.3">
<span style="color:blue">This is wrong
target
</span>
Target
year quantity min val max
1 2011 f NA 0.4550558 NA
2 2012 f NA 0.3200000 NA
3 2012 catch 0 NA 40000
4 2012 f NA <font color="#ff0000"><a href="tel:2.0000000" value="+4520000000" target="_blank">2.0000000</a></font> NA
min val max
1 NA 0.45506 NA
2 NA 0.32000 NA
3 0.00000 NA 40000.00000
4 NA NA <font color="#ff0000">2.00000</font>
Try this! I have added an argument to fwd maxF that sets the upper bound on F
print("OK to here1")
## project MP to TAC year with the HCR and geom mean recr.
c.=fwdControl(data.frame(val=c(rep(0.2,2),NA),min=c(NA,NA,0),max=c(NA,NA,40000),year=c(2011:2012,2012),quantity=c(rep("f",2),"catch")))
MP<-fwd(MP,ctrl=c.,sr=list(model="mean",params=FLPar(1)),sr.residuals=SRrs,maxF=2) ## HERE IS THE PROBLEM
print("OK to here2")
</pre>
</span></span><br>
Laureie<div><br>
On 03/21/2012 05:15 PM, Luis Ridao wrote:
</div><div><div><blockquote type="cite">here it comes ... sorry about it
<div><br>
</div>
<div>best,</div>
<div>Luis<br>
<br>
<div class="gmail_quote">On Wed, Mar 21, 2012 at 4:11 PM, Laurie
<span dir="ltr"><<a href="mailto:lauriekell@googlemail.com" target="_blank">lauriekell@googlemail.com</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> There is no funcs.R
file<span><font color="#888888"><br>
<br>
Laurie</font></span>
<div>
<div><br>
<br>
On 03/20/2012 01:09 PM, Luis Ridao wrote:
<blockquote type="cite">
<div><br>
</div>
I have attached the input files and the script (the
workspace is too large >60 MB)
<div><br>
<div><br>
</div>
<div>Best,</div>
<div>Luis<br>
<br>
<div class="gmail_quote">On Tue, Mar 20, 2012 at
11:57 AM, Laurie <span dir="ltr"><<a href="mailto:lauriekell@googlemail.com" target="_blank">lauriekell@googlemail.com</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> Send
me the objects and code and I will try and
run it.<br>
<br>
Laurie
<div>
<div><br>
<br>
On 03/20/2012 12:21 PM, Luis Ridao
wrote: </div>
</div>
<blockquote type="cite">
<div>
<div>FLR-help,
<div><br>
</div>
<div>I'm having problems with a MSE.</div>
<div>I copy-paste the line which
gives the error:</div>
<div><br>
</div>
<div>...</div>
<div>...</div>
<div>
<div> <b> print("OK to here1")</b></div>
<div><b> ## project MP to TAC
year with the HCR and geom
mean recr.</b></div>
<div><b>
MP<-fwd(MP,ctrl=target,sr=list(model="mean",params=FLPar(1)),sr.residuals=SRrs)</b></div>
</div>
<div>...</div>
<div> ...</div>
<div><br>
</div>
<div>
<div>Year =2011</div>
<div>[1] "OK to here1"</div>
<div><b>Error in .local(object,
ctrl, ...) : 6</b></div>
</div>
<div><br>
</div>
<div><br>
</div>
<div>My target is :</div>
<div>
<div>> <b>target</b></div>
<div><br>
</div>
<div>Target</div>
<div> year min val max
quantity season area unit spp
fleet metier rel.year</div>
<div>1 2011 NA 0.4550558 NA
f NA NA NA NA
NA NA NA</div>
<div>2 2012 NA 0.3200000 NA
f NA NA NA NA
NA NA NA</div>
<div>3 2012 0 NA 40000
catch NA NA NA NA
NA NA NA</div>
<div>4 2012 NA <a href="tel:2.0000000" value="+4520000000" target="_blank">2.0000000</a>
NA f NA NA NA
NA NA NA NA</div>
<div> rel.season rel.area
rel.unit</div>
<div> 1 NA NA
NA</div>
<div>2 NA NA
NA</div>
<div>3 NA NA
NA</div>
<div>4 NA NA
NA</div>
<div> </div>
<div> min val max
</div>
<div> 1 NA 4.5506e-01
NA</div>
<div> 2 NA 3.2000e-01
NA</div>
<div> 3 0.0000e+00 NA
4.0000e+04</div>
<div> 4 NA NA
2.0000e+00</div>
</div>
<div><br>
</div>
<div>the '<b>sr</b>' object :</div>
<div><br>
</div>
<div>
<div>>
list(model="mean",params=FLPar(1))</div>
<div>$model</div>
<div>[1] "mean"</div>
<div><br>
</div>
<div>$params</div>
<div>An object of class "FLPar"</div>
<div>params</div>
<div>a </div>
<div>1 </div>
<div>units: NA </div>
</div>
<div><br>
</div>
<div>and the Sr residuals:</div>
<div><br>
</div>
<div>
<div>> <b>SRrs</b></div>
<div>An object of class "FLQuant"</div>
<div>iters: 1000 </div>
<div><br>
</div>
<div>, , unit = unique, season =
all, area = unique</div>
<div><br>
</div>
<div> year</div>
<div>quant 2011 2012 2013
</div>
<div> all 27355(0) 27355(0)
27355(0)</div>
<div><br>
</div>
<div>units: NA </div>
</div>
<div><br>
</div>
<div>R version :
R.2.14.2</div>
<div>
<div>FLash Version: 2.5.0</div>
</div>
<div>FLAdvice Version: 1.0</div>
<div>
<div>FLXSA Version: 2.5</div>
</div>
<div>
<div>FLAssess Version: 2.5</div>
</div>
<div>FLCore Version: 2.5</div>
<div><br>
</div>
<div>Running on Ubuntu 10.04.04</div>
<div><br>
</div>
<div>Thanks in advance</div>
<div><br>
</div>
<br>
<fieldset></fieldset>
<br>
</div>
</div>
<pre>_______________________________________________
flr-list mailing list
<a href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a>
</pre>
</blockquote>
<br>
</div>
<br>
_______________________________________________<br>
flr-list mailing list<br>
<a href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
flr-list mailing list
<a href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a></pre>
</blockquote>
<br>
</div>
</div>
</div>
<br>
_______________________________________________<br>
flr-list mailing list<br>
<a href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a><br>
<br>
</blockquote>
</div>
<br>
</div>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
flr-list mailing list
<a href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a></pre>
</blockquote>
<br>
</div></div></div>
<br>_______________________________________________<br>
flr-list mailing list<br>
<a href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a><br>
<br></blockquote></div><br>