Esiste una funzione singola, simile a "runif", "rnorm" e simili che produce previsioni simulate per un modello lineare? Posso codificarlo da solo, ma il codice è brutto e presumo che questo sia qualcosa che qualcuno ha già fatto.Esiste una funzione o un pacchetto che simulerà le previsioni per un oggetto restituito da lm()?
slope = 1.5
intercept = 0
x = as.numeric(1:10)
e = rnorm(10, mean=0, sd = 1)
y = slope * x + intercept + e
fit = lm(y ~ x, data = df)
newX = data.frame(x = as.numeric(11:15))
Quello che mi interessa è una funzione che sembra la linea di seguito:
sims = rlm(1000, fit, newX)
Tale funzione sarebbe tornato 1000 simulazioni di valori y, in base alle nuove variabili x.
L'ultima riga nel Q mi ha confuso. 'x' è fisso; intendi simulare 'y' (la risposta) per i nuovi dati' x'? –
Scusa, Gavin, hai ragione. Intendevo dire che le risposte sarebbero state simulate. Questo è stato modificato. – PirateGrunt
OK, quindi si potrebbe guardare a '? Simulate' ma funziona solo con' x' corrente. Ma potreste alterarlo ('simulate.lm()') per chiamare 'predicict()' sull'oggetto del modello con 'newdata = newX' invece della chiamata corrente a' fitted() 'e quindi permettergli di procedere come con il codice normale. Supponendo che i "pesi" non siano stati usati in quanto complicherebbe le cose ... –