Tässä artikkelissa aiomme nähdä kuinka laskea standardivirhe R-ohjelmointikielessä.
Matemaattisesti voimme laskea keskivirheen käyttämällä kaavaa:
keskihajonta/neliöjuuri(n)
R-kielessä voimme laskea seuraavilla tavoilla:
- sd()-funktion käyttö pituusfunktion kanssa
- Käyttämällä standardivirhekaavaa.
- Plotrix-paketin käyttö.
Tapa 1: sd()-funktion käyttö pituusfunktion kanssa
Tässä aiomme käyttää sd()-funktiota, joka laskee keskihajonnan ja sitten pituus()-funktion löytääkseen havaintojen kokonaismäärän.
Syntaksi: sd(data)/sqrt(length((data)))
Esimerkki: R-ohjelma, joka laskee standardivirheen 10 arvon joukosta vektorissa
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))))> |
valitse nimellä
>
>
Lähtö:
[1] 26.20274>
Tapa 2: Käyttämällä standardivirhekaavaa
Tässä käytetään standardivirhekaavaa havaintojen saamiseen.
Syntaksi: sqrt(sum((a-mean(a))^2/(pituus(a)-1)))/sqrt(pituus(a))
missä
- data on syöttötieto
- sqrt-funktio on neliöjuuren etsiminen
- summaa käytetään datan elementtien summan selvittämiseen
- keskiarvo on funktio, jota käytetään tietojen keskiarvon löytämiseen
- pituus on funktio, jota käytetään palauttamaan datan pituus
Esimerkki: R-ohjelma, joka laskee standardivirheen kaavan avulla
R
merkkijono ti int
# 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)))> |
>
>
Lähtö:
java-tyyppinen muuttuja
[1] 26.20274>
Tapa 3: std.error()-funktion käyttäminen (plotrix-paketti)
Tämä on sisäänrakennettu toiminto, joka laskee suoraan vakiovirheen. Se on saatavana plotrix-paketissa
Syntaksi: std.error(data)
Esimerkki: R-ohjelma standardivirheen laskemiseen käyttämällä 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))> |
>
>
Lähtö:
[1] 26.20274>