Displaying 1 result from an estimated 1 matches for "fglm".
Did you mean:
fgl
2011 Jul 14
1
glm() scale parameters and predicted Values
...y <- list(yr=yr,n=n)
require(utils)
ny <- expand.grid(ny)
y = rgamma(100, shape=1.5, rate = 1, scale = 2)
Gdata <- cbind(ny,y)
Gdata2<- Gdata
Gdata$x1 <- cos((3.14*yr)/365.25)
Gdata$x2 <- sin((3.14*yr)/365.25)
#Fitting Generalized Linear Models
Gdata <- split(Gdata,Gdata$n)
FGLM <- lapply(Gdata, function(x){
????????????? m <- as.numeric(x$y)
????????????? x1 <- m <- as.numeric(x$x1)
????????????? x2 <- m <- as.numeric(x$x2)
????????????? summary(glm(m~1+x1+x2, family=Gamma),dispersion=NULL)
?????????????? })
#Save the results of the estimated parameter...