Pokouším se vytvořit predikční intervaly pomocí funkce predict() pro novou sadu dat, ale na více než jeden model, který jsem vytvořili pro dataset. Jsem relativně nezkušený v používání lapply, ale postava to by mělo být užitečné v tomto procesu:
#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
)
Nicméně, toto čerpá následující chyba:
Chyba v(jako.charakteru(ZÁBAVA), mode = "funkce", envir = žp) : objekt 'model_dna' režim 'funkce' nebyl nalezen
Nejsem si jistý, jak nejlépe strukturovat tento lapply funkce, zejména při použití po více než jeden model. Je obecně čistší způsob, jak na to?