Displaying 4 results from an estimated 4 matches for "tissuedata".
2009 Nov 12
0
writing selfStart models that can deal with treatment effects
...ecayAndConstantInflowModel,
expDecayAndConstantInflowModelInit, c("aL", "aN", "T0"))
Ignoring the treatment effects, the routines seem to work find:
-----------------------------------
> getInitial(Count ~ SSexpDecayAndConstantInflow(B0, Time, aL, aN, T0), data =
tissueData)
aL aN T0
4.600144e-02 1.000000e-04 1.082172e+03
> nls(Count ~ SSexpDecayAndConstantInflow(B0, Time, aL, aN, T0), data =
tissueData)
Nonlinear regression model
model: Count ~ SSexpDecayAndConstantInflow(B0, Time, aL, aN, T0)
data: tissueData
aL...
2009 Oct 15
0
Setting random effects within a category using nlme
...al. 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 296.4755 0.09270158 0.09542535
Naive/...
2009 Oct 17
1
custom selfStart model works with getInitial but not nls
...s)
}
##Now create new model with selfStart attributes
const.PBMC.tcell.modelSS<- selfStart(model = const.PBMC.tcell.model,
initial=const.PBMC.tcell.selfStart)
##Test routines using getInitial -- This works
> getInitial(Count ~ const.PBMC.tcell.modelSS(B0, Time,aL, aN, T0), data =
> tissueData)
[1] 0.05720924
$aL
[1] 0.05720924
$aN
[1] 0.1981895
$T0
[1] 1360.292
##Now try to use the SS model -- this doesn't work
> nls(Count ~ const.PBMC.tcell.modelSS(B0, Time,aL, aN, T0), data = tissueData)
Error in numericDeriv(form[[3L]], names(ind), env) :
Missing value or an infinity pro...
2009 Oct 15
2
Proper syntax for using varConstPower in nlme
...al. 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/IngLN 2114.1868 0.05298126 0.12795589...