Wie eine Schleife durch eine matrix Zeilen von i:n und Spalten

Ich bin noch ein Anfänger in r. Ich versuche eine Schleife durch eine matrix und berechnen von Summen looping aus jeder einzelnen Zeile(1:5), bis der Schwellenwert (z.B. 14) in dieser Zeile erreicht ist. Dann gehen Sie durch alle Spalten (2:6). Was ich will ist, dass meine matrix(b) ist gefüllt mit den Nummern(91:95), wo die Zeilen haben den Schwellenwert erreicht. Ich habe versucht, diese. Hat jemand eine Idee, was ich falsch mache.

a <- matrix(c(91:95,1:25), ncol=6)    
b <- matrix(NA,ncol=6,nrow=5)


for(i in 1:nrow(a)) {
    for(j in 2:ncol(a)) {
        for(n in 1:5) {
            if(sum(a[seq(i:n),j])>14) {
                b[i,j] <- a[n,1]
                break
            }
        }
    }
}

Dies ist mein Ergebnis.

      [,1] [,2] [,3] [,4] [,5] [,6]
[1,]   NA   95   93   92   91   91
[2,]   NA   NA   94   91   91   91
[3,]   NA   NA   91   91   91   91
[4,]   NA   NA   91   91   91   91
[5,]   NA   91   91   91   91   91

Mein gewünschtes Ergebnis sollte wie folgt Aussehen. Die Ergebnisse sollten auf den Startpunkt(i) und zeigen Sie die Zeile (n) wo der Schwelle(14) wurde überschritten.

      [,1] [,2] [,3] [,4] [,5] [,6]
[1,]   NA   95   93   92   91   91
[2,]   NA   NA   93   93   91   91
[3,]   NA   NA   94   94   91   91
[4,]   NA   NA   95   95   91   91
[5,]   NA   NA   NA   95   91   91

Kann mir jemand helfen mit diesem.

Vielen Dank im Voraus. Ist auch einfacher, die Ideen als geschachtelte Schleifen sind geschätzt.

Nicht sicher, ob ich Verstand Ihr problem, aber in der 4. Zeile in der Schleife sollte nur i:n oder seq(i,n) und nicht seq(i:n)

InformationsquelleAutor user3472591 | 2015-01-06

loops r

1

Hier ist der Ansatz ähnlich wie das von @josiber aber mit dplyr (vorausgesetzt, es ist ein Fehler in die erwartete Ausgabe). Der erste block von code, den ich.e bis zu der Erstellung von "df1" ist im Grunde erstellen von Replikationen von "0" - Zeilen für eine Gruppierung der variable "ind". Wir erstellen ein n Wiederholungen von "a" (rep(list(a), n)) dann rbind es um die form der matrix "a1", ein index matrix "i1" von "1s" die gleichen Abmessungen wie "a1". Das "i1" matrix aufgefüllt mit "0" - Zeilen basierend auf dem "m1" - index, und erstellen Sie eine Daten.Rahmen("df1") mit grouping variable "ind" und die anderen Spalten wird das Ergebnis von "a1*i1", so dass Elemente in "i1", das sind "0" "0" ausgegeben.
```
library(dplyr)
a1 <- do.call(rbind,rep(list(a),n))
i1 <- matrix(1,ncol=n,nrow=nrow(a1))
m1 <- matrix(seq(6,nrow(a1), by=n)+rep(0:3,each=4),ncol=4,byrow=TRUE)
i1[m1[upper.tri(m1,diag=TRUE)],] <- 0
df1 <- data.frame(ind=rep(1:n,each=n),a1*i1)
```
Verwendung von "df1", verwenden wir "ind", wie die Gruppierung der Variablen (group_by) und fasse die Spalten ("X1:X5") mit summarise_each. Erhalten wir den index der Elemente in jeder Spalte, die eine kumulative Summe größer als 14, und verwenden Sie diese als index, um die Elemente von "a" (o[which(cumsum(.)>14)..)
```
df1 %>%
    group_by(ind)%>%
    summarise_each(funs(o[which(cumsum(.)>14)[1]]))
#  ind X1 X2 X3 X4 X5
#1   1 95 93 92 91 91
#2   2 NA 93 93 92 92
#3   3 NA 94 94 93 93
#4   4 NA 95 95 94 94
#5   5 NA NA 95 95 95
```
Benötigen oder wenn Sie eine matrix-Ausgang, anstelle des "dplyr" erstellen Sie die "indx" mit Map und verwenden tapply um das Ergebnis zu erhalten.
```
indx <- do.call(rbind,Map(`+`,list(matrix(1:5,ncol=5,
                       nrow=5,byrow=TRUE)), seq(0,20, by=5)) )
res <- matrix(tapply(a1*i1, indx, FUN=function(x)
              o[which(cumsum(x) >14)[1]]),ncol=5,byrow=TRUE)
 res
 #     [,1] [,2] [,3] [,4] [,5]
 #[1,]   95   93   92   91   91
 #[2,]   NA   93   93   92   92
 #[3,]   NA   94   94   93   93
 #[4,]   NA   95   95   94   94
 #[5,]   NA   NA   95   95   95
```
Daten
```
a <- matrix(1:25, nrow=5)
o <- 91:95
n <- 5
```
- Dank josilber und akrun für Eure Antworten. Ja, Sie hatten Recht, es war ein Fehler in meinem Beispiel. Ich sehe deine Antworten ganz gut funktionieren. Nichtsdestotrotz bin ich aber immer noch nur versuchen, durch zu bekommen Sie Ihre Antwort, und herauszufinden, wie Sie Ihre Funktion(x) und die Funktion(p) genau funktionieren und was Sie berechnen.
- Ich aktualisiert mit einigen Erklärungen. Es ist erstaunlich, dass Sie kommentiert, nach so vielen Tagen. Ich habe vergessen, die meisten der Schritte 🙂
- Hey akrun. Vielen Dank. Ich werde in Sie Graben. Ja, ich war damit beschäftigt, mit einigen anderen arbeiten. Grüßt.
InformationsquelleAutor akrun
2

Als ich Ihre Frage verstehe, ist das element (i, j) von Ihrem Ergebnis enthalten soll die Zeile identifier für die erste Zeile, wo die kumulative Summe (i, j) + (i+1, j) + (i+2, j) + ... übersteigt 14, und es sollte NA wenn das kumulative Summe nie mehr als 14 vor dem erreichen der unteren Rand des input-matrix. In diesem Fall ist die erwartete Ausgabe enthält einen Fehler -- die beiden rechten Spalten werden sollte, 91, 92, 93, 94, 95 statt 91, 91, 91, 91, 91.

Für jede Spalte der input-matrix (wir arbeiten mit apply mit Marge 2) Sie können eine for-Schleife (ich benutze sapply mehr als eine komfortable alternative unten), Berechnung der kumulierten Summe der Spalte ignorieren-element über die element-Sie kalkulieren:
```
a <- matrix(1:25, nrow=5)
o <- 91:95
apply(a, 2, function(x) sapply(1:length(x), function(p) {
    cs <- cumsum(x * c(rep(0, p-1), rep(1, length(x)-p+1)))
    if (max(cs > 14)) {
        return(o[min(which(cs > 14))])
    } else {
        return(NA)
    }
}))
#      [,1] [,2] [,3] [,4] [,5]
# [1,]   95   93   92   91   91
# [2,]   NA   93   93   92   92
# [3,]   NA   94   94   93   93
# [4,]   NA   95   95   94   94
# [5,]   NA   NA   95   95   95
```
InformationsquelleAutor josliber

Schreibe einen Kommentar

Du musst angemeldet sein, um einen Kommentar abzugeben.

Daten