Displaying 1 result from an estimated 1 matches for "alval".
Did you mean:
alva
2009 Oct 17
1
custom selfStart model works with getInitial but not nls
...t1, "Count"]);
T2 = mean(data[data$Time==t2, "Count"]);
if(T0 < T2){ ##increase -- doesn't work
stop(paste("Error in const.PBMC.tcell.start: T0 < T2 for data: ", data[1,
]));
}
##Estimate aL based on exponential decline from t=0 to t=24
aLVal = -(log(T1) - log(T0))/(t1-t0);
##Estimate aNVal based on final value
aNVal = aLVal*T2/B0;
values = list(aLVal, aNVal, T0);
names(values) <- mCall[c("aL", "aN", "T0")]; #mimic syntax used by P&B
return(values)
}
##Now create new model with se...