Is there a way to get `cor` to work?

Offen

#90 erstellt am 15.4.2021 0 Kommentare 0 Reaktionen 0 zugewiesene Personen Auf GitHub ansehen

Dieses Issue hat noch niemand übernommen.

Bewertung

Schwierigkeit: 4/5
Geschätzter Aufwand: 3-5 Tage
Anfängerfreundlichkeit: 42/100

Issue-Typ: Feature
Klarheit: Größtenteils klar
Aktivitätsstatus: Veraltet
Tech-Stack: r
Bereich: data

Rechercherichtung

Beginnen Sie damit, den gemeldeten Fehler cor(DelayedArray(x), DelayedArray(y)) zu reproduzieren, und vergleichen Sie ihn mit den bereitgestellten delayedCov-Beispielen. Verwenden Sie das im Issue beschriebene Verhalten der Pearson-Korrelation als Ziel und überprüfen Sie, dass äquivalente Matrix- und DelayedArray-Eingaben übereinstimmende Ergebnisse ohne den numerischen Typfehler liefern.

Vom Indexierungsmodell aus dem Issue-Text verfasst.

Beschreibung

When I try cor, I get errors like Error in cor(DelayedArray(x), DelayedArray(y)) : 'x' must be numeric. Am I missing something simple?

If it's not implemented, I have a little function for this and I could send a pull request.

#' Drop-in replacement for cov or cor; tolerates DelayedArray input. Pearson correlation only. 
#'
delayedCov = function(x, y = x, do_cor = F){
  mux = colMeans(x)
  muy = colMeans(y)
  zx = sweep(x, 2, mux, "-")
  zy = sweep(y, 2, muy, "-")
  covxy = t(zx) %*% zy
  covxy = covxy / nrow(x)
  if(do_cor){
    sx = sqrt(colMeans(zx^2))
    sy = sqrt(colMeans(zy^2))
    covxy = Matrix::Diagonal(x = 1/sx) %*% covxy %*% Matrix::Diagonal(x = 1/sy)
  }
  covxy
}
# tests
x = matrix(rnorm(1000), ncol = 10)
y = matrix(rnorm(1000), ncol = 10)
plot(delayedCov(x, y),  cov(x, y)); abline(a = 0, b = 1)
plot(delayedCov(x, y, do_cor = T),  cor(x, y)); abline(a = 0, b = 1)
plot(delayedCov(DelayedArray(x), DelayedArray(y), do_cor = T),  cor(x, y)); abline(a = 0, b = 1)

Vorherrschende Sprache: R
Sterne: 29
Forks: 12
PR-Merge-Kennzahlen: Keine gemergten PRs in 30 T.

Beitragsleitfaden

Für dieses Repository ist kein Beitragsleitfaden indexiert

Erste Schritte

Lesen Sie das ganze Issue und danach den Beitragsleitfaden des Projekts.
Schreiben Sie ins Issue, dass Sie es übernehmen — das erspart doppelte Arbeit.
Forken Sie das Repository und arbeiten Sie in einem Branch.
Öffnen Sie einen Pull Request, der die Issue-Nummer nennt.

Mehr aus Bioconductor/DelayedArray

Conversion to dgCMatrix Causes Segmentation Fault Offen

Schwierigkeit 4/5 3-5 Tage Anfängerfreundlichkeit 38/100

Bioconductor/DelayedArray#129 · 11 Kommentare · Erstellt am 10.9.2026
Custom delayed operations Offen

Schwierigkeit 4/5 3-5 Tage Anfängerfreundlichkeit 45/100

Bioconductor/DelayedArray#127 · 1 Kommentar · Erstellt am 22.5.2026
the supplied seed must support extract_array() Offen

Schwierigkeit 4/5 3-5 Tage Anfängerfreundlichkeit 25/100

Bioconductor/DelayedArray#125 · 1 Kommentar · Erstellt am 13.10.2025
Summation of ConstantArray of 0 with vector of 0's fails Offen

Schwierigkeit 3/5 1-2 Tage Anfängerfreundlichkeit 45/100

Bioconductor/DelayedArray#123 · Erstellt am 05.3.2025
Proposal to migrate `DelayedMask` to the DelayedArray package Offen

Schwierigkeit 5/5 Über eine Woche Anfängerfreundlichkeit 20/100

Bioconductor/DelayedArray#122 · Erstellt am 23.1.2025

Alle Issues in Bioconductor/DelayedArray