Displaying 1 result from an estimated 1 matches for "tidal1".
Did you mean:
tidal
2011 Jul 11
1
GLS - Plotting Graphs with 95% conf interval
...SE=TRUE argument.
Here is my code so far:
f1<-formula(MaxNASC40_50~hu3+flcmax+TidalFlag)
vf1Exp<-varExp(form=~hu3)
B1D<-gls(f1,correlation=corGaus(form=Lat~Lon, nugget=TRUE),weights=vf1Exp ,
data=ocean)
ochu3<-seq(from=2.9,to=4,length=120)
ocflc<-seq(from=0,to=0.8,length=120)
tidal1<-rep(c("1"),length=120)
mydata1<-data.frame(TidalFlag=factor(tidal1),hu3=ochu3,flcmax=ocflc)
lineNASC1<-predict(B1D,newdata=mydata1,type="response")
lineNASC1<-as.data.frame(lineNASC1)
plot(ocean$MaxNASC40_50[ocean$TidalFlag==1]~ocean$flcmax[ocean$TidalFlag==1)
lin...