2016-06-09 17 views
5

Come posso calcolare la variazione della popolazione dei miei dati utilizzando R?Varianza popolazione in r

ho letto c'è un pacchetto chiamato popvar ma ho la versione 0.99.892 e non trovo il pacchetto

+1

Che software avete della versione 0.99.892 di? L'attuale versione R è 3.3.0. Forse RStudio? Non confondere R e RStudio - RStudio è solo uno strumento per scrivere il codice R. – Gregor

+2

La funzione 'var' di R si divide per n-1 per impostazione predefinita. Moltiplicare l'output di var per (n-1)/n ti darà quello che vuoi. – Dason

risposta

11

La funzione var() nella R base calcola la varianza campionaria e la varianza della popolazione differisce dalla varianza campionaria di un fattore di n/n - 1. Così un'alternativa per calcolare varianza della popolazione sarà var(myVector) * (n - 1)/n dove n è la lunghezza del vettore, ecco un esempio:

x <- 1:10 
var(x) * 9 /10 
[1] 8.25 

Dalla definizione di varianza della popolazione:

sum((x - mean(x))^2)/10 
[1] 8.25 
1

È possibile trovare i dettagli sul pacchetto popvar qui: https://cran.r-project.org/web/packages/PopVar/index.html - È possibile installarlo utilizzando il comando install.packages("PopVar"); Si noti che il nome è case sensitive (maiuscola P, maiuscola V).

+0

Sarà utile se spiegherai come usarlo. Ho provato ** PopVar (c (1,2,3)) ** e ho ottenuto un errore: * Errore: impossibile trovare la funzione "PopVar" *. –

0

si dispone già di un grande rispondi, ma mi piacerebbe mostrare che puoi facilmente fare le tue funzioni di convenienza. È sorprendente che una variazione di popolazione/funzione di deviazione standard non sia disponibile nella base R. È disponibile in Excel/Calc e in altri software. Non sarebbe difficile avere una tale funzione. Potrebbe essere nominato sdp o sd.p o essere richiamato con sd(x, pop = TRUE)

Ecco una versione di base di varianza della popolazione senza alcun tipo di controllo:

x <- 1:10 
    varp <- function(x) mean((x-mean(x))^2) 
    varp(x) 
    ## [1] 8.25 

Per scalare, se la velocità è un problema, colSums e/o si possono usare i moli (vedi: https://rdrr.io/r/base/colSums.html)

+0

Sorprendente inoltre che non ci siano asimmetria (con/senza regolazione) o kurtosi o momenti grezzi in base '' R'' ... – PatrickT