Hi, Below is the code that is giving me this error message:
#=================================================================== #The code below results in error: Error: NA/NaN/Inf in foreign function #call (arg 1) #=================================================================== A=rep(c(-1,1),8);B=rep(c(-1,-1,1,1),4);C=rep(c(-1,-1,-1,-1,1,1,1,1),2) D=c(rep(-1,8),rep(1,8)) y=c(38.01697,42.04596,43.89677,47.78144,37.34869,41.90744,44.01468,47.96458,37.97263,41.83482,44.01936,48.12203,38.73267,42.92054,44.16163,48.88900) bcols=c(1,2) dcols=c(3,4) locmodel=glm(y~A+B) mlcoefs=locmodel$coef res=locmodel$residuals ressq=res^2 dispmodel=glm(ressq~C+D,family=Gamma(link = "log")) #=================================================================== #The code below results in error message: # Warning message: In glm.fit(x = X, y = Y, weights = weights, start = # start, etastart = #etastart, : algorithm did not converge #=================================================================== A=rep(c(-1,1),8) B=rep(c(-1,-1,1,1),4) C=rep(c(-1,-1,-1,-1,1,1,1,1),2) D=c(rep(-1,8),rep(1,8)) y=c(37.93499,42.19333,44.02388,48.04979,37.52459,41.90613,43.44424,47.81314, 38.04942,42.06124,44.00249,47.89033,39.00372,40.47911,43.78126,48.55228) locmodel=glm(y~A+B) mlcoefs=locmodel$coef res=locmodel$residuals ressq=res^2 dispmodel=glm(ressq~C+D,family=Gamma(link = "log")) #=============================== I researched on this error and it can be caused by missing values, which is not the case for this example. I would really appreciate any help that can be extended. Rainier Sabangan BGSU MSAS Student ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.