Displaying 1 result from an estimated 1 matches for "lag_estimation".
2008 Aug 11
3
Peoblem with nls and try
...od2)
c<-coeff["c"]
A<-coeff["A"]
B<-coeff["B"]
e<-coeff["e"]
k<-coeff["k"]
summary(nlmod2)
plot(P1)
lines(P1[,1],nlmod2,lwd=2,col=2)
lag1 <- e-(2/k)
print(lag1)
abline(v=lag1,col=4)
print("nlmod2")
{
pdf(file= "Lag_estimation.pdf" )
plot(P1[,1],predict(nlmod2),xlab='Time',ylab='ThT normalizate',main="Lag estimetion",type='l', col=1)
lines(P1, col=8)
abline(v=lag1,col=4)
lines(P1[,1],nlmod2,lwd=2,col=2)
}
The loops don't` function. the program stop after the first cycl...