Displaying 1 result from an estimated 1 matches for "lamintr1".
Did you mean:
laminar
2012 May 08
1
Error with psi value for 'segmented' package for R
...with an error message and I can't find solution. It is connected with psi value, and the error says:
Error in seg.glm.fit(y, XREG, Z, PSI, weights, offs, opz) :
(Some) estimated psi out of its range
This is the code I am using:
library("segmented")
curva<-read.table("lamintr1.txt", header=T)
attach(curva)
fit.glm<-glm(gpp~temp, weight=NULL, family=gaussian)
plot(temp,gpp,xlab=expression("temp"),
ylab="gpp",pch=15,cex=0.8,xlim=c(0,50), ylim=c(0,40))
o1<-glm(gpp ~ temp, weight=NULL, family=gaussian)
os1<-segmented(o1, seg.Z=~temp, psi=15,...