Šajā rakstā mēs redzēsim, kā aprēķināt standarta kļūdu R programmēšanas valodā.
Matemātiski mēs varam aprēķināt standarta kļūdu, izmantojot formulu:
standarta novirze/kvadrātsakne(n)
R Valodā mēs varam aprēķināt šādos veidos:
- Funkcijas sd() izmantošana ar garuma funkciju
- Izmantojot standarta kļūdu formulu.
- Izmantojot plotrix pakotni.
1. metode: funkcijas sd() izmantošana ar garuma funkciju
Šeit mēs izmantosim funkciju sd (), kas aprēķinās standarta novirzi un pēc tam funkciju garums (), lai atrastu kopējo novērojumu skaitu.
Sintakse: sd(dati)/sqrt(garums((dati)))
Piemērs: R programma, lai aprēķinātu standarta kļūdu no 10 vērtību kopas vektorā
R
# consider a vector with 10 elements> a <->c>(179, 160, 136, 227, 123, 23,> >45, 67, 1, 234)> # calculate standard error> print>(>sd>(a)/>sqrt>(>length>((a))))> |
atgriežot masīvu java
>
>
stīgu maiņa c
Izvade:
[1] 26.20274>
2. metode: izmantojot standarta kļūdu formulu
Šeit novērojumu iegūšanai izmantosim standarta kļūdu formulu.
Sintakse: kvadrāts(summa((a-mean(a))^2/(garums(a)-1)))/sqrt(garums(a))
kur
- dati ir ievades dati
- sqrt funkcija ir atrast kvadrātsakni
- summa tiek izmantota, lai atrastu elementu summu datos
- vidējais ir funkcija, ko izmanto, lai atrastu datu vidējo vērtību
- garums ir funkcija, ko izmanto, lai atgrieztu datu garumu
Piemērs: R programma standarta kļūdas aprēķināšanai, izmantojot formulu
R
# consider a vector with 10 elements> a <->c>(179, 160, 136, 227, 123, 23,> >45, 67, 1, 234)> # calculate standard error> print>(>sqrt>(>sum>((a ->mean>(a)) ^ 2/(>length>(a) - 1)))> >/>sqrt>(>length>(a)))> |
>
>
Izvade:
[1] 26.20274>
3. metode: funkcijas std.error() izmantošana (plotrix pakotne)
Šī ir iebūvētā funkcija, kas tieši aprēķināja standarta kļūdu. Tas ir pieejams plotrix iepakojumā
Sintakse: std.error(data)
Piemērs: R programma standarta kļūdas aprēķināšanai, izmantojot std.error()
R
# import plotrix package> library>(>'plotrix'>)> # consider a vector with 10 elements> a <->c>(179,160,136,227,123,> >23,45,67,1,234)> # calculate standard error using in built> # function> print>(>std.error>(a))> |
supw
>
>
Izvade:
[1] 26.20274>