search for: 2df

Dear fellows, I'm trying to extract the AIC statistic from a GLM model with quasipoisson link. The formula I'm referring to is AIC = -2(maximum loglik) + 2df * phi with phi the overdispersion parameter, as reported in: Peng et al., Model choice in time series studies os air pollution and mortality. J R Stat Soc A, 2006; 162: pag 190. Unfortunately, the function logLik doesn't work for a quasipoisson link. Do you know a fast method to extract t...

...put (shown at the bottom). I have two questions: 1) What exactly are the p-values in the Pr(>|z|) representing? I understand that the coefficients have standard errors, etc., but i'm not sure how the p-value there is calculated. 2) At the bottom, where it shows the results of an LRT with 2df, i don't quite understand what model the ratio is being tested against. If the current model has two variables (x1 and x2), and those are the extra degrees of freedom, then the baseline should then have 0 variables, but that's not really a Cox model? thanks for any help. Brian > sum...

...-> number of animais alive (the response) n -> total number of animals To test proportion of alive, make d/n. In GLIM: $yvar d$ $error binomial n$ $fit +h$ scale deviance = 25.730 (change = -9.138) at cycle 4 d.f. = 15 (change = -2) factor h is significant by chisq? with 2df fron tables = 5.99 $disp e$ estimate se parameter 1 -0.1054 0.2055 1 2 0.7985 0.2961 h(2) 3 0.08827 0.26764 h(3) Scale parameter taken as 1.000 In R: > modelo.glex24.1 <- glm((d/n)~h,family=binomial) Warning message: non-integer #successes in a binomial glm! in:...

...An example of what I'm looking for would be the "K-squared" statistic that tests for normality (D'Agostino and Pearson 1973, Biometrika 60: 613, also in Zar, 1996, Biostatistical Analysis, p89). The expected distribution of the test statistic is approximately chi-squared with 2df, regardless of values of estimated parameters or sample size (provided sample size is sufficiently large). Petr's reply: >Maybe it is easier and sufficient to use the Kolmogorov - Smirnov >test, that is implemented as ks.test() in R. However, I am not able >to check your reference...

Hi all, Does anyone have an R implementation of the queplot (K.?M. Goldberg and B.?Iglewicz. Bivariate extensions of the boxplot. Technometrics, 34(3):pp. 307?320, 1992)? I'm struggling with the estimation of the asymmetry parameters. Hadley -- Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/

....3165 house3 1.58 1 530.37 < 2e-16 *** house4 0.76 1 254.75 < 2e-16 *** Residuals 8.86 2974 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 GVIF Df GVIF^(1/2Df) as.factor(parkfee) 1.570 2 1.119 fare 1.562 1 1.250 as.factor(ohouse) 1.063 3 1.010 tripMIN 1.110 1 1.054 as.factor(Devt2) 2.378 3 1.155 cityarea 4.188 1 2.046 coreDistance 7.037 1...

Hello, I want to use the quantile function so I read the doc but I don't understand with this > qchisq(seq(0.05,0.95,by=0.05),df=(length(don)-1)) [1] 62667.11 62795.62 62882.42 62951.47 63010.74 63064.00 63113.39 63160.27 63205.65 63250.33 63295.04 63340.48 63387.48 63437.03 63490.53 63550.14 63619.68 [18] 63707.24 63837.16 Can you help me please?