2016-06-09 8 views
5

Jak mogę obliczyć wariancję populacji moich danych za pomocą R?Wariancja populacji w r

czytam tam jest pakiet o nazwie popvar ale mam w wersji 0.99.892 i nie znajdę pakiet

+1

Jakie oprogramowanie posiada wersja 0.99.892? Aktualna wersja R to 3.3.0. Może RStudio? Nie mylić R i RStudio - RStudio to tylko narzędzie do pisania kodu R. – Gregor

+2

Funkcja 'var' R' dzieli domyślnie przez n-1. Pomnożenie wyniku przez var (n-1)/n da ci to, czego chcesz. – Dason

Odpowiedz

11

Funkcjaw podstawie R oblicza wariancję próbki, a wariancja populacji różni się w zależności od wariancji o współczynnik n/n - 1. Więc alternatywą obliczyć wariancję populacji będzie var(myVector) * (n - 1)/n gdzie n to długość wektora, oto przykład:

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

Z definicji wariancji populacji:

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

można znaleźć szczegółowe informacje na temat pakietu popvar tutaj: https://cran.r-project.org/web/packages/PopVar/index.html - Można ją zainstalować za pomocą polecenie install.packages("PopVar"); Należy zwrócić uwagę, że w nazwie rozróżniana jest wielkość liter (kapitał P, kapitał V).

+0

Pomocne będzie wyjaśnienie, jak z niego korzystać. Próbowałem ** PopVar (c (1,2,3)) ** i otrzymałem komunikat o błędzie: * Błąd: nie można znaleźć funkcji "PopVar" *. –

0

masz już wielki odpowiedź, ale chciałbym pokazać, że możesz łatwo tworzyć własne funkcje wygody. Zaskakujące jest to, że funkcja odchylenia populacji/odchylenia standardowego nie jest dostępna w bazie R. Jest dostępny w Excel/Calc i innym oprogramowaniu. Nie byłoby trudno mieć taką funkcję. To może być nazwany sdp lub sd.p lub być wywołany z sd(x, pop = TRUE)

Oto podstawowej wersji wariancji populacji bez typu kontroli:

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

skalować w górę, jeśli prędkość jest problem, colSums i/lub colMeans mogą być używane (patrz: https://rdrr.io/r/base/colSums.html)

+0

Również zaskakujące, że nie ma skośności (z/bez korekty) lub kurtozy lub surowych momentów w bazie '' R'' ... – PatrickT