Calcular error estándar en R

En este artículo vamos a ver cómo calcular el error estándar en el lenguaje de programación R. 

Matemáticamente podemos calcular el error estándar usando la fórmula:

desviación estándar/raíz cuadrada(n)

En lenguaje R, podemos calcular de estas formas:

• Usando la función sd() con la función de longitud
• Usando la fórmula del error estándar.
• Usando el paquete plotrix.

Método 1: usar la función sd() con la función de longitud

Aquí vamos a utilizar la función sd() que calculará la desviación estándar y luego la función length() para encontrar el número total de observaciones.

Sintaxis: sd(datos)/sqrt(longitud((datos)))

Ejemplo: programa R para calcular un error estándar a partir de un conjunto de 10 valores en un vector

# 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))))

Producción:

[1] 26.20274

Método 2: mediante el uso de la fórmula de error estándar

Aquí usaremos la fórmula del error estándar para obtener las observaciones.

Sintaxis: sqrt(suma((a-media(a))^2/(longitud(a)-1)))/sqrt(longitud(a))

dónde

• los datos son los datos de entrada
• la funcion sqrt es encontrar la raiz cuadrada
• sum se utiliza para encontrar la suma de los elementos en los datos
• mean es la función utilizada para encontrar la media de los datos
• length es la función utilizada para devolver la longitud de los datos

Ejemplo: programa R para calcular el error estándar usando la fórmula

# 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)))

Producción:

[1] 26.20274

Método 3: Uso de la función std.error() (paquete plotrix)

Esta es la función incorporada que calculó directamente el error estándar. Está disponible en paquete plotrix 

Sintaxis: std.error(datos)

Ejemplo: programa R para calcular el error estándar usando std.error()

# 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))

Producción:

[1] 26.20274