Dividir una variable según compartimientos de otra variable

Question

He continua de datos "A", binario de datos categóricos "O", género/sexo y la edad de varios de los participantes en un estudio.

Un modelo lineal en la I muestra que no hay correlación entre la Una y la edad. Ahora me gustaría grupo en grupos por edad y ver si hay una diferencia entre los grupos. Yo sé acerca de 'hist' y 'split' en R, pero estos no hacen lo que yo necesito.

(1) ¿Cómo puedo dividir/dividir en grupos según su edad (de 18 a 27 28 a 37, etc.).

(2) una Vez que he hecho eso, me gustaría utilizar una $\chi^2$ prueba?

(3) también Podría probar O en los mismos grupos, utilizando la cuenta?

Answer 1

> A   <- round(rnorm(100, 100, 15), 2)       # generate some data
> age <- sample(18:65, 100, replace=TRUE)
> sex <- factor(sample(0:1, 100, replace=TRUE), labels=c("f", "m"))

# 1) bin age into 4 groups of similar size
> ageFac <- cut(age, breaks=quantile(age, probs=seq(from=0, to=1, by=0.25)),
+               include.lowest=TRUE)

> head(ageFac)
[1] (26,36.5] (26,36.5] (36.5,47] [18,26]   [18,26]   [18,26]  
Levels: [18,26] (26,36.5] (36.5,47] (47,65]

> table(ageFac)   # check group size
ageFac
[18,26] (26,36.5] (36.5,47]   (47,65] 
     27        23        26        24

# 2) test continuous DV in age-groups
> anova(lm(A ~ ageFac))
Analysis of Variance Table
Response: A
          Df  Sum Sq Mean Sq F value Pr(>F)
ageFac     3    15.8   5.272  0.0229 0.9953
Residuals 96 22099.2 230.200               

# 3) chi^2-test for equal distributions of sex in age-groups    
> addmargins(table(sex, ageFac))
     ageFac
sex   [18,26] (26,36.5] (36.5,47] (47,65] Sum
  f        11        10        12      11  44
  m        16        13        14      13  56
  Sum      27        23        26      24 100

> chisq.test(table(sex, ageFac))
        Pearson's Chi-squared test
data:  table(sex, ageFac) 
X-squared = 0.2006, df = 3, p-value = 0.9775