<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
OK<br>
<br>
I will have a look, but busy running stuff at the moment.<br>
<br>
This might be of interest. <br>
<br>
Laurie<br>
<br>
On 03/22/2012 12:33 PM, Luis Ridao wrote:
<blockquote
cite="mid:CAL09-Ef3v2VQ8HWreZJ83RoUfFmV3_kLLWQ8EWWGK32qr_G02A@mail.gmail.com"
type="cite">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 moz-do-not-send="true"
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 moz-do-not-send="true" 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
moz-do-not-send="true"
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
moz-do-not-send="true"
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
moz-do-not-send="true"
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 moz-do-not-send="true" href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a>
<a moz-do-not-send="true" 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 moz-do-not-send="true"
href="mailto:flr-list@flr-project.org"
target="_blank">flr-list@flr-project.org</a><br>
<a moz-do-not-send="true"
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 moz-do-not-send="true" href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a>
<a moz-do-not-send="true" 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 moz-do-not-send="true"
href="mailto:flr-list@flr-project.org"
target="_blank">flr-list@flr-project.org</a><br>
<a moz-do-not-send="true"
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 moz-do-not-send="true" href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a>
<a moz-do-not-send="true" 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 moz-do-not-send="true"
href="mailto:flr-list@flr-project.org" target="_blank">flr-list@flr-project.org</a><br>
<a moz-do-not-send="true"
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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
flr-list mailing list
<a class="moz-txt-link-abbreviated" href="mailto:flr-list@flr-project.org">flr-list@flr-project.org</a>
<a class="moz-txt-link-freetext" href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/flr-list</a>
</pre>
</blockquote>
<br>
</body>
</html>