Displaying 2 results from an estimated 2 matches for "solb".
Did you mean:
sol
2004 Jul 18
2
bootstrap and nls
Hi,
I am trying to bootstrap the difference between each parameters among two non linear regression (distributed loss model) as
following:
# data.frame
> Raies[1:10,]
Tps SolA Solb
1 0 32.97 35.92
2 0 32.01 31.35
3 1 21.73 22.03
4 1 23.73 18.53
5 2 19.68 18.28
6 2 18.56 16.79
7 3 18.79 15.61
8 3 17.60 13.43
9 4 14.83 12.76
10 4 17.33 14.91
etc...
# non linear model (work well)
RaiesLossA.nls<-nls(SolA~a/(1+b*Tps)^c,start=c(a=32,b=0.5,c=0.6))...
2004 Mar 12
0
Basic questions on nls and bootstrap
...C0*exp(-k*Tpst) with Ct the concentration of pollutant at time t, C0 is the
initial concentration (i.e. t=0), Tps is the time in days.
The table containing data is called «tabMika»
Here, you can find my R program:
> tabMika<-read.delim("RMika.txt")
> tabMika
Tps SolA Solb
1 0 32.97 35.92
2 0 32.01 31.35
3 1 21.73 22.03
4 1 23.73 18.53
5 2 19.68 18.28
6 2 18.56 16.79
# and continue like that until 29 days.
> library(nls)
> attach(tabMika)
> Expon<-function(Tps,parm){
+ C0<-parm[1]
+ k<-parm[2]
+ }
> DegSA.nls<-nls(SolA...