Displaying 2 results from an estimated 2 matches for "gdata2".
Did you mean:
data2
2011 Jul 11
3
fitdistr() Error
...#Grouped vector
n <- c(1:100)
yr <-c(1:100)
ny <- list(yr=yr,n=n)
require(utils)
ny <- expand.grid(ny)
y = rgamma(1000, shape=1.5, rate = 1, scale = 2)
Gdata <- cbind(ny,y)
#MLE Estimation of Gamma Distribution Parameters
library(MASS)
#Generate starting values
y <- as.numeric(Gdata2$y)
me? <- mean(y)
sde <- sd(y)
sh? <- sqrt(me/sde)
sc? <- sqrt(sde)/me
Gdata <-split(Gdata,Gdata$n)
parm <- lapply(Gdata, function(x){
????????? y <- as.numeric(x$y)
?????????? fitdistr(y ,"gamma",list(shape=sh.mom, scale=sc),
#method = c("Nelder-Mead", &quo...
2011 Jul 14
1
glm() scale parameters and predicted Values
...te any help.
?
?
?
#Call required R packages
require(plyr)?
require(stats)
require(fitdistrplus)
require(MASS)
#Grouped vector
n <- c(1:10)
yr <-c(1:10)
ny <- 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 &l...