Displaying 2 results from an estimated 2 matches for "degsa".
Did you mean:
degs
2004 Apr 09
1
bootstrap function coefficients
..., 10
estimations of the parameters C0 and k from the function: SolA =
C0*exp(-k*Tps).
# First, I fit the data and all is OK
> tabMika<-read.delim("RMika.txt")
> library(nls)
> attach(tabMika)
> Expon<-function(Tps,parm){
+ C0<-parm[1]
+ k<-parm[2]
+ }
> DegSA.nls<-nls(SolA~C0*exp(-k*Tps),start=c(C0=35, k=1),tabMika)
> summary(DegSA.nls)
Formula: SolA ~ C0 * exp(-k * Tps)
Parameters:
Estimate Std. Error t value Pr(>|t|)
C0 25.682104 1.092113 23.52 < 2e-16 ***
k 0.087356 0.007582 11.52 6.36e-13 ***
---
Signif. codes: 0 `***...
2004 Mar 12
0
Basic questions on nls and bootstrap
...a
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~C0*exp(-k*Tps),start=c(C0=35, k=1),tabMika)
> summary(DegSA.nls)
Formula: SolA ~ C0 * exp(-k * Tps)
Parameters:
Estimate Std. Error t value Pr(>|t|)
C0 25.682104 1.092113 23.52 < 2e-16 ***
k 0.087356 0.007582 11.52 6.36e-13 ***
---
Signif. codes: 0 `***...