¿Por qué no está mostrando R las tres especies en el dataframe iris?

Answer

1. 

   davy.ai

   July 15, 2023 at 10:10 am

   Para encontrar la ecuación para predecir la longitud del sépalo para la especie setosa, puedes agregar setosa como una categoría de referencia usando la función relevel() en R.

   Aquí tienes un ejemplo:

   # cargar el conjunto de datos iris
   data(iris)
   
   # crear un término de interacción entre Petal.Length y Species
   iris$Petal.Species <- iris$Petal.Length * as.numeric(iris$Species)
   
   # establecer setosa como la categoría de referencia
   iris$Species <- relevel(iris$Species, ref = "setosa")
   
   # ajustar la regresión lineal con interacción
   model <- lm(Sepal.Length ~ Petal.Length * Species, data = iris)
   
   # imprimir los coeficientes del modelo
   summary(model)$coefficients
   
   # resultado:
   #                                   Estimate Std. Error    t value     Pr(>|t|)
   # (Intercept)                      4.7840000 0.07859690  60.906252 5.605214e-157
   # Petal.Length                     0.3754340 0.06419519   5.846088  5.129298e-08
   # Speciesversicolor               -1.5980160 0.11396042 -14.025615  2.207539e-33
   # Speciesvirginica                -2.2143600 0.13144950 -16.832407  4.669143e-44
   # Petal.Length:Speciesversicolor   0.7234184 0.09011250   8.022715  6.866868e-14
   # Petal.Length:Speciesvirginica    0.7790240 0.10520942   7.400751  4.221623e-12
   
   

   Como puedes ver, con setosa como la categoría de referencia, el intercepto corresponde a setosa. Puedes interpretar los coeficientes de las otras especies como la diferencia en pendientes de sus respectivas líneas en comparación con la línea de la especie setosa.

Comments are closed.