[R] help

Thierno 12 Sow thiernosow12 at gmail.com
Mon Feb 1 18:03:46 CET 2010


>
> Dear Madam, Sir,
> Dear Madam, Sir,
> I am a student in statistics at the Université Rennes 1 (France), I've used
> the polr function in R to estimate an ordered probit model (consumer
> preference) but the strange thing is that *I did not get the probabilities
> which can tell me whether the variables are significant or not*. My model
> has a dep var with 6 levels (0 to 5) and 19 explanatory variables (both
> dummy and continues). Below is the code that I used and the result that I
> got. Also I enclosed a part of the data used in the model for further help.
> I am waiting for your help to fix this issue, many thanks.
> *******Code****
> m <- sqlQuery(channel = 1, select * from [Feuil2$])
> m$Depvar <- as.factor(m$Depvar)
> library(MASS, pos=4)
> prm <- polr(Depvar ~ var1 + var2 + var3 + var4 + var5 +var6 + var7 + var8 +
>
>   var9, method="probit", data=m, Hess=TRUE)
> summary(prm)
>
> *****Output*****
> > m <- sqlQuery(channel = 1, select * from [Feuil2$])
> > m$Depvar <- as.factor(m$Depvar)
> > library(MASS, pos=4)
> > prm <- polr(Depvar ~ var1 + var2 + var3 + var4 + var5 +var6 + var7 + var8
> +
> +   var9, method="probit", data=m, Hess=TRUE)
> > summary(prm)
> Call:
> polr(formula = Depvar ~ var1 + var2 + var3 + var4 + var5 + var6 +
>     var7 + var8 + var9, data = m, Hess = TRUE, method = "probit")
> Coefficients:
>            Value Std. Error    t value               *how to get P-values
> ?
> *var1  0.04454384 0.10390217  0.4287094                *?
> *var2  0.22172411 0.10221326  2.1692305                *?*
> var3 -0.19858516 0.09909994 -2.0038879                *?*
> var4  0.48806601 0.10375119  4.7041967                *?*
> var5  0.01845784 0.10018124  0.1842444                *?*
> var6 -0.37819633 0.10170134 -3.7186955                *?*
> var7  0.06047232 0.03234284  1.8697284                *?*
> var8  0.39855835 0.19211257  2.0746084                *?*
> var9  0.16119371 0.04692679  3.4350044                *?*
> Intercepts:
>     Value   Std. Error t value
> 0|1 -0.2063  0.2931    -0.7040
> 1|2  0.3162  0.2903     1.0889
> 2|3  0.7557  0.2900     2.6059
> 3|4  1.2209  0.2913     4.1915
> 4|5  1.7293  0.2944     5.8736
> Residual Deviance: 1667.097
> AIC: 1695.097
>  Thank you for your time
> sincerely
> Thierno, Université de Rennes 1 France
>


More information about the R-help mailing list