search for: glm1_glm

...13-233-5101 University of Leeds fax (44) 113-233-5090 Leeds LS2 9JT, England e-mail: J.T.Kent@leeds.ac.uk Input: x_1:4 # regressor variable y_c(2,6,7,8) # response binomial counts n_rep(10,4) # number of binomial trials ym_cbind(y,n-y) # response variable as a matrix glm1_glm(ym~x,binomial) # fit a generalized linear model f_fitted(glm1) rp1_(y-n*f)/sqrt(n*f*(1-f)) # direct calculation of pearson residuals rp2_residuals(glm1,type="pearson") # should be pearson residuals rd_residuals(glm1) # deviance residuals cbind(rp1,rp2,rd) # note second column is wrong, in...