Displaying 1 result from an estimated 1 matches for "raieslossa".
Did you mean:
raieslossb
2004 Jul 18
2
bootstrap and nls
...del) 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))
RaiesLossB.nls<-nls(Solb~a/(1+b*Tps)^c,start=c(a=33,b=1.5,c=0.5))
# bootstrap
library(boot)
dif.param<-function(data,i){
RaiesLossA.nls<-nls(SolA[i]~a/(1+b*Tps[i])^c,start=c(a=31,b=0.5,c=0.6))
RaiesLossB.nls<-nls(Solb[i]~a/(1+...