wie durch jede Zeile in r

Question

type<-c("111","111","111","222","222","222","222","333","333","333","444","444","444") 

location<-c("aaa","bbb","ccc","aaa","bbb","ddd","eee","aaa","ccc","ddd","ccc","ddd","fff") 

datafr1 <- data.frame(type,location) 

gemeinsamen Wert proportional erhalten möchte ich das Ergebnis wie folgt aus. (Brauchen nicht ein Format so sein.)

  222  333 444 
111  0.66% 0.66% 0.33% 

und 222, 333, 444 auch. Finden Sie die gemeinsamen Werte und erhalten Sie dann Anteil der Schnittmenge.

Können Sie mir helfen? !!

Answer 1

Wir können auf den table Ausgang von ‚datafr1‘ ein Kreuzprodukt tun, dann teilen jede Zeile von dem Maximalwert der Ausgabe für die rbindlist am Ende

m1 <- tcrossprod(table(datafr1)) 
m1/matrixStats::rowMaxs(m1) 
#  type 
#type  111  222  333  444 
# 111 1.0000000 0.6666667 0.6666667 0.3333333 
# 222 0.5000000 1.0000000 0.5000000 0.2500000 
# 333 0.6666667 0.6666667 1.0000000 0.6666667 
# 444 0.3333333 0.3333333 0.6666667 1.0000000 

Answer 2

I nur data.table verwendet werde. Sie können dies durch eine do.call auf rbind oder etwas ähnliches ersetzen, wenn Sie nicht data.table verwenden möchten.

typevals <- unique(datafr1$type) 

funky <- function(typeval) { 
    uniqueloc <- unique(datafr1[type == typeval,]$location) 
    props <- sapply(typevals, 
        function(x) { 
        length(intersect(uniqueloc, datafr1[type == x,]$location))/length(uniqueloc) 
        }) 
    props <- round(props, 2) 
    df <- data.frame(t(props)) 
    names(df) <- typevals 
    df <- cbind(data.frame(type = typeval), df) 
} 
rbindlist(lapply(typevals, funky)) 
# type 111 222 333 444 
# 1: 111 1.00 0.67 0.67 0.33 
# 2: 222 0.50 1.00 0.50 0.25 
# 3: 333 0.67 0.67 1.00 0.67 
# 4: 444 0.33 0.33 0.67 1.00