On 2012-07-02 02:37, Jennifer Kaiser wrote:
Hi,
I have a problem with the df.
I read in a big csv file.
Tabelle <-
read.csv("C:\\Users\\Public\\Documents\\Bachelorarbeit\\eingabe8_durchnummeriert.csv" ,
header = T , sep=";")
then I try this:
ygamma <- glm(Tabelle$sb_ek_ber ~1+ Tabelle$FAHRL_C + Tabelle$NUTZKREIS +
Tabelle$schw_drittel_c , family = Gamma)
anova(ygamma, test="Chisq")
Analysis of Deviance Table
Model: Gamma, link: inverse
Response: Tabelle$sb_ek_ber
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev Pr(>Chi)
NULL 1236805 35451551
Tabelle$FAHRL_C 1 33987 1236804 35417564 0.0018493 **
Tabelle$NUTZKREIS 1 48903 1236803 35368661 0.0001880 ***
Tabelle$schw_drittel_c 1 47328 1236802 35321334 0.0002388 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
str(Tabelle)
'data.frame': 1236806 obs. of 9 variables:
$ Alter_Jüngster_C_inkl_AlterNutz: int 1 1 1 1 1 1 1 1 1 1 ...
$ ALTERKAU_C : int 1 2 2 1 3 3 3 4 1 1 ...
$ FAHRL_C : int 1 2 1 3 4 3 3 1 5 1 ...
$ NUTZKREIS : int 1 2 2 2 2 2 2 1 1 2 ...
$ RKL_U12 : int 1 1 1 2 3 4 4 3 5 6 ...
$ SF_Sonder_aufgefüllt : int 1 2 3 4 4 4 4 5 6 7 ...
$ schw_drittel_c : int 1 2 3 4 3 3 3 3 1 1 ...
$ sb_ek_ber : num 0.001 0.001 0.001 0.001 0.001 0.001
0.001 0.001 0.001 0.001 ...
$ JE_gewichtet : num 0.384 3.952 3.952 2.81 3.952 ...
I don't understand why the df are always 1.
You probably intended all the variables that are of
type "integer" (e.g. FAHRL_C) to be _factors_. My guess
is that, for ease of data entry, you coded these with
integers 1-7.
You'll have to tell R that you want factors:
Tabelle$FAHRL_C <- factor(Tabelle$FAHRL_C)
etc.
Peter Ehlers
it would be great if you could help me.
[[alternative HTML version deleted]]
______________________________________________
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.