Estoy tratando de generar intervalos de predicción utilizando la función de predecir() para un conjunto de datos nuevos, pero a través de más de un modelo que he generado para un conjunto de datos. Yo soy relativamente inexpertos en el uso de lapply, pero la figura que debe ser útil en este proceso:
#Calling in my libraries:
library(dplyr)
#Creating dataset:
DNase <- DNase
#Generating models, one for each "Run" in DNAse:
model_dna <- DNase %>%
group_by(Run) %>%
do(model_dna_group = lm(log(density) ~ log(conc), data = .)) %>% ungroup()
#Creating a new data set to be used to generate predictions:
new_dna <- as.data.frame(DNase$conc) %>%
mutate(conc = DNase$conc * 2) %>% select(conc)
#Attempting to apply predict to these models for a new data frame:
new_dna_w_predictions <- lapply(
X = model_dna,
FUN = predict,
newdata = new_dna,
interval = "prediction",
level = 0.9
)
Sin embargo, este saca el siguiente error:
Error en llegar(como.carácter(DIVERSIÓN), mode = "función", no = no) : el objeto 'model_dna' de modo de 'función' no se ha encontrado
No estoy seguro de cómo mejor la estructura de esta lapply función, especialmente cuando se utilizan en más de un modelo. Hay generalmente un limpiador de manera acerca de esto?