Eu sunt încercarea de a genera intervale de predicție folosind funcția prezice() pentru un nou set de date, dar peste mai mult de un model pe care l-au generat pentru un set de date. Eu sunt relativ lipsit de experiență în utilizarea lapply, dar cifra ar trebui să fie de ajutor în acest proces:
#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
)
Cu toate acestea, acest lucru atrage următoarea eroare:
Eroare în ia(ca.caracter(DISTRACȚIE), modul = "funcția", mediu = mediu) : obiect 'model_dna' de modul de "funcție" nu a fost găsit
Eu nu sunt sigur cum cel mai bine pentru a structura acest lapply funcție, mai ales atunci când sunt utilizate în mai mult de un model. Există, în general, o mai curat mod de a aborda acest lucru?