Czy jest jakieś polecenie, aby znaleźć błąd standardowy średniej w R?W R, jak znaleźć standardowy błąd średniej?
64
A
Odpowiedz
25
113
Standardowy błąd jest tylko odchylenie standardowe podzielone przez pierwiastek kwadratowy z wielkością próby. Dzięki czemu można łatwo zrobić własną funkcję:
> std <- function(x) sd(x)/sqrt(length(x))
> std(c(1,2,3,4))
[1] 0.6454972
74
Prawdopodobnie jest bardziej wydajny w użyciu var ... skoro faktycznie sqrt dwa razy w kodzie, raz aby uzyskać sd (kod do SD w badania i ujawnił poprzez wpisanie „SD”) ...
se <- function(x) sqrt(var(x)/length(x))
0
bardziej ogólnie, do standardowych błędów na każdym innym parametrem, można skorzystać z pakietu startowego dla symulacji bootstrap (lub zapisać je na własną rękę)
51
wersja powyższej odpowiedzi Johna, która usuwa brzydkie NA:
stderr <- function(x) sqrt(var(x,na.rm=TRUE)/length(na.omit(x)))
3
Pakiet sciplot ma wbudowaną funkcję se (x)
-10
y <- mean(x, na.rm=TRUE)
sd(y)
dla odchylenia standardowego var(y)
dla wariancji.
Obie pochodne używają w mianowniku n-1
, więc bazują na danych przykładowych.
Powiązane problemy
- 1. R obliczyć błąd standardowy użyciem bootstrap
- 2. Błąd: nie można znaleźć funkcji w R ...
- 3. Jak znaleźć standardowy katalog site_perl dla Perla?
- 4. Jaki algorytm R używa do obliczania średniej?
- 5. Jak obliczyć błąd standardowy z wyników ODR?
- 6. obliczenie średniej wartości w tabelach za pomocą wzoru [R]
- 7. Sortowanie boxplot od średniej (a nie średnia) w R
- 8. R: Obliczanie średniej z 5 lat w danych panelu
- 9. Jak znaleźć Waldo z R?
- 10. Jak obliczyć wartość P i standardowy błąd z analizy korelacji R cor()
- 11. Standardowy mechanizm ma opcje specyficzne dla pakietu w R?
- 12. Błąd kompilacji szablonu - standardowy czy nie?
- 13. Co oznacza "standardowy interfejs formuły do danych.frame" w języku R?
- 14. Wysoka wydajność macierzy średniej
- 15. Jak znaleźć maksimum avg
- 16. Średniej wielkości czcionki w css
- 17. Błąd: nieoczekiwane "}" w "}" w R
- 18. Standardowy błąd ignorowania NaN w grupie pandy według grup
- 19. Obliczanie średniej kroczącej w F #
- 20. Obliczanie średniej ruchomej w C++
- 21. Gdzie można znaleźć najnowszy standardowy dokument ANSI C
- 22. znaleźć i zastąpić sekwencję liczbową w R
- 23. Średniej wielkości aplikacja Clojure?
- 24. Ogólne - Czy istnieje sposób obliczania średniej na podstawie istniejącej średniej i nowych parametrów dla średniej wynikowej?
- 25. "Standardowy czas wschodni USA" a "Czas wschodnio-standardowy" w .NET
- 26. Standardowy wywoływacz SQL UPSERT
- 27. najbliższa wartości średniej SQL
- 28. R rvest: nie mógł znaleźć funkcję „xpath_element”
- 29. R - Jak znaleźć punkty w obrębie określonego konturu
- 30. Jak znaleźć największy N elementów na liście w R?
Co ciekawe, twoja funkcja i Ian są prawie identyczne szybko. Testowałem je zarówno 1000 razy w porównaniu do 10^6 milionów losowań rnorm (niewystarczająca moc, aby popchnąć je mocniej niż to). Odwrotnie, funkcja plotrix była zawsze wolniejsza niż nawet najwolniejsze przebiegi tych dwóch funkcji - ale ma też o wiele więcej działań pod maską. –
Jest tak prawdopodobnie dlatego, że przewaga funkcji Johna pozwala uniknąć jednego wywołania funkcji i obliczenia pierwiastka kwadratowego z jednej liczby. Jeśli ma jakąkolwiek przewagę, jest bardzo mały. –
Nie tylko prawdopodobne, ale prawdziwe ... Robiłem taktowanie rnorm (5), ale wykonałem 5e4 egzekucji polecenia se. Moja wersja zajmuje około 2,5 sekundy, podczas gdy ta z dwoma sqrts zajmuje około 2.8. Niewielkie przyspieszenie. (Wykonałem kilka uruchomień i raportuję średnią zaokrągloną) – John