Displaying 4 results from an estimated 4 matches for "pbmc".
Did you mean:
bmc
2009 Oct 17
1
custom selfStart model works with getInitial but not nls
...al to the selfStart model, I get results. However, if I try
usint it with nls or nlsList, these routines complain about a lack of initial
conditions.
If someone could point out what I'm doing wrong, I'd greatly appreciate it.
Details:
## Nonlinear model I want to fit to the data
const.PBMC.tcell.model <- function(B0, t, aL, aN, T0){
Tb0 = B0;
x = exp(-log(aL) + log(T0*aL+(-1+exp(t * aL))*Tb0 * aN) - t * aL);
return(x);
}
##Define selfStart routine
const.PBMC.tcell.selfStart<- function(mCall, LHS, data){
t0 = 0;
t1 = 24;
t2 = 48;
##Get B0 Value
B0...
2009 Oct 15
0
Setting random effects within a category using nlme
...hus giving me 8 datasets in total. The dataset is
in a groupedData object with
formula= Count ~ Time|Type/Tissue
and there are three basic model parameters that I am trying to fit: T0, aN,
and aL.
Calling nlsList gives
------------------------------------------
>nlsList(Count ~ quad.PBMC.model(aL, aN, T0), data = tissueData, start =
list(T0 = 1000, aN = exp(-2), aL = exp(-2)))
Call:
Model: Count ~ quad.PBMC.model(aL, aN, T0) | Type/Tissue
Data: tissueData
Coefficients:
T0 aN aL
Naive/CLN 1530.0088 0.26508876 0.04730525
Naive/ILN...
2009 Oct 15
2
Proper syntax for using varConstPower in nlme
...hus giving me 8 datasets in total. The dataset is in a
groupedData object with
formula= Count ~ Time|Type/Tissue
and there are three basic model parameters that I am trying to fit: T0, aN,
and aL.
Calling nlsList gives
------------------------------------------
>nlsList(Count ~ quad.PBMC.model(aL, aN, T0), data = tissueData, start =
list(T0 = 1000, aN = exp(-2), aL = exp(-2)))
Call:
Model: Count ~ quad.PBMC.model(aL, aN, T0) | Type/Tissue
Data: tissueData
Coefficients:
T0 aN aL
Naive/CLN 1530.0088 0.26508876 0.04730525
.
.
.
Memory/I...
2008 Sep 05
0
help for color parameter
...how can I specify the color I want in this situation and code. Basically, I would like to have 10 different colors of line and pch(line and pch will be in same color) for each combination of ID and DID. I expected first 10 colors in R for this plot initially.
I appreciate your instructions.
dogpbmc=read.table('E:\\MPI\\4924\\PBMC-PKPD\\mpbmc.tab', skip=1, header=T)
dogpbmc=dogpbmc[dogpbmc$MDV!=1,]
par(mfrow=c(2,2),mar=c(5,5,4,4))
# E vs t for simultaneous modeling
temp.col=0
for (i in unique(dogpbmc$ID)){
j=dogpbmc$ID==i
plot(dogpbmc$TIME[j], dogpbmc$DV[j], type='n', xlim=c(0,...