Displaying 2 results from an estimated 2 matches for "svcom".
Did you mean:
secom
2012 Sep 04
3
Comparing Von Bertalanffy Growth Curves
...2L, 2012L, 2012L, 2012L, 2012L)), .Names = c("Age", "MM", "Year"), class = "data.frame", row.names = c(NA, -37L))
?
In case it's helpful here is all my code before this point:
str(LMB)
MMi=as.integer(MM)
Yearf=as.factor(Year)
Agei=as.integer(Age)
( svCom <- vbStarts(MMi~Agei,data=LMB))
( svGen <- lapply(svCom,rep,2) )
vbGen <- MMi~Linf[Yearf]*(1-exp(-K[Yearf]*(Age-t0[Yearf]))+error)
fitGen <- nls(vbGen,data=LMB,start=svGen)
?
Thank you, April
2012 Sep 10
0
More help need on Von Bertalanffy Growth Curves
...", header=T)
attach(BLG)
names(BLG)
dput(BLG,"")
$ Von Bertalanffy Growth Model$
source("http://www.rforge.net/NCStats/InstallNCStats.R")
source("http://www.rforge.net/FSA/InstallFSA.R")
library(FSA)
library(NCStats)
library(nlstools)
# provide starting values
svCom <- list(Linf=199.8132, K=0.5631745, t0=-1.960147)
svCom <- vbStarts(Length~Age,data=BLG)
# repeat the starting values for each Cohort
svGen <- lapply(svCom, rep, length(unique(Cohort)))
# define a new variable, uniquely numbered for each Cohort
Cohortf <- as.factor(Cohort)
# fit a si...