2015-12-18 7 views
5

ich einen Datenrahmen haben, die etwa wie folgt aussieht:Erstellen eines Two-Mode-Frequenzmatrix in R

CASENO Var1 Var2 Resp1 Resp2 
1   1  0  1  1 
2   0  0  0  0 
3   1  1  1  1 
4   1  1  0  1 
5   1  0  1  0 

Es gibt mehr als 400 Variablen im Datensatz. Dies ist nur ein Beispiel. Ich muss eine einfache Frequenzmatrix in R erstellen (ausgenommen die Fallzahlen), aber die table Funktion funktioniert nicht. Insbesondere möchte ich einen Teil der Spalten tabellarisch darstellen, um eine Zwei-Modus-Matrix von Frequenzen zu erstellen. Die Tabelle sollte wie folgt aussehen:

 Var1 Var2 
Resp1 3  1 
Resp2 3  2 

In Stata, lautet der Befehl:

gen var = 1 if Var1==1 
replace var= 2 if Var2==1 

gen resp = 1 if Resp1==1 
replace resp = 2 if Resp2==1 

tab var resp 

Antwort

5

man sollte für eine beliebige Anzahl von Var arbeiten & Resp:

d <- structure(list(CASENO = 1:5, Var1 = c(1L, 0L, 1L, 1L, 1L), Var2 = c(0L, 0L, 1L, 1L, 0L), Resp1 = c(1L, 0L, 1L, 0L, 1L), Resp2 = c(1L, 0L, 1L, 1L, 0L)), .Names = c("CASENO", "Var1", "Var2", "Resp1", "Resp2"), class = "data.frame", row.names = c(NA, -5L)) 

m <- as.matrix(d[,-1]) 
m2 <- t(m) %*% m 
rnames <- grepl('Resp',rownames((m2))) 
cnames <- grepl('Var',colnames((m2))) 
m2[rnames,cnames] 

[UPDATE] Eine elegantere Version, in dem Kommentar von G.Grothendieck bereitgestellt:

m <- as.matrix(d[,-1]) 
cn <- colnames(m); 
crossprod(m[, grep("Resp", cn)], m[, grep("Var", cn)]) 
+2

Eine weitere Vereinfachung wäre 'm <- as.matrix (d)' da der greps nie sowieso die erste Spalte entsprechen. –

+0

Danke! Das ist so hilfreich. Wie würde ich mit dem Befehl crossprod Spaltennummern statt Spaltennamen referenzieren? – jj987246

+0

@ jj987246, verwende einfach Vektoren, die Spaltennummern enthalten, z.B. 'crossprod (m [, 1: 4], m [, 5: 8])' –

4

Ich bin sicher, es gibt einen anderen Weg, aber man könnte dies tun: Diese

library(reshape2) 
library(plyr) 

df1 <- melt(df[,-1],id=1:2) 
ddply(df1,.(variable),summarize, 
     Var1 = sum(value==1&Var1==1), 
     Var2 = sum(value==1&Var2==1)) 

# variable Var1 Var2 
# 1 Resp1 3 1 
# 2 Resp2 3 2 
3

Hier ist ein Ansatz mit xtabs.

# get names of non "variables" 
not_vars <- c("Resp1", "Resp2", "CASENO") 

# get names of "variables" 
vars <- as.matrix(d[,!names(d) %in% not_vars]) 

# if you have many more than 2 response variables, this could get unwieldy 
result <- rbind(
    xtabs(vars ~ Resp1, data=d, exclude=0), 
    xtabs(vars ~ Resp2, data=d, exclude=0)) 

# give resulting table appropriate row names.  
rownames(result) <- c("Resp1", "Resp2") 
#  Var1 Var2 
#Resp1 3 1 
#Resp2 3 2 

Beispieldaten:

d <- read.table(text=" 
CASENO Var1 Var2 Resp1 Resp2 
1   1  0  1  1 
2   0  0  0  0 
3   1  1  1  1 
4   1  1  0  1 
5   1  0  1  0", header=TRUE)