Displaying 1 result from an estimated 1 matches for "rainfallgps".
2011 Dec 12
3
shorter way of coding
...alues for for y (y1, y2, ...y300) which means I will have to write the code below 300 times to generate the 300 different fitted values for y. Is there a short way of doing that ?
Many thanks in advance
Mintewab
library(mgcv)
dat <- read.table("e:/minti's laptop/C/GBG/allround_survey/rainfallGPS.csv", header=T, sep=",")
b<-gam(y1~s(x1, x2, k=100),data=dat)
vis.gam(b)
fitted(b)