Kembali mentransformasikan hasil regresi saat memodelkan log (y)

Saya menyesuaikan regresi pada . Apakah valid untuk mendukung estimasi titik transformasi (dan interval kepercayaan / prediksi) dengan eksponensial? Saya tidak percaya begitu, karena tetapi menginginkan pendapat orang lain. $\log(y)$ $E[f(X)] \ne f(E[X])$

Contoh saya di bawah ini menunjukkan konflik dengan transformasi balik (.239 vs .219).

set.seed(123)

a=-5
b=2

x=runif(100,0,1)
y=exp(a*x+b+rnorm(100,0,.2))
# plot(x,y)

### NLS Fit
f <- function(x,a,b) {exp(a*x+b)} 
fit <- nls(y ~ exp(a*x+b),  start = c(a=-10, b=15)) 
co=coef(fit)
# curve(f(x=x, a=co[1], b=co[2]), add = TRUE,col=2,lwd=1.2) 
predict(fit,newdata=data.frame(x=.7))
[1] 0.2393773

### LM Fit
# plot(x,log(y))
# abline(lm(log(y)~x),col=2)
fit=lm(log(y)~x)
temp=predict(fit,newdata=data.frame(x=.7),interval='prediction')
exp(temp)
        fit       lwr       upr
1 0.2199471 0.1492762 0.3240752

regression back-transformation Lembah kecil
sumber

Apakah ini bukan salah satu masalah yang dipecahkan oleh GLM gaussian log-linked?

generic_user

@ARM Ya saya percaya begitu. Terima kasih telah menunjukkannya. Namun menggunakan GLM lebih sulit untuk mendapatkan interval prediksi tapi saya pikir saya bisa menyelesaikannya.

Glen

@Glen Lakukan pencarian untuk mengolesi Duan di situs ini.

Dimitriy V. Masterov

Itu tergantung pada apa yang ingin Anda peroleh di ujung yang lain.

Interval kepercayaan untuk parameter yang ditransformasikan mentransformasikannya dengan baik. Jika memiliki cakupan nominal pada skala log, maka akan memiliki cakupan yang sama kembali pada skala asli, karena monotonitas transformasi.

Interval prediksi untuk pengamatan di masa depan juga mentransformasikannya dengan baik.

Interval untuk rata-rata pada skala log umumnya tidak akan menjadi interval yang cocok untuk rata-rata pada skala asli.

Namun, kadang-kadang Anda dapat secara tepat atau sekitar menghasilkan estimasi yang masuk akal untuk rata-rata pada skala asli dari model pada skala log.

Namun, perawatan diperlukan atau Anda mungkin menghasilkan perkiraan yang memiliki sifat yang agak mengejutkan (mungkin untuk menghasilkan perkiraan yang tidak memiliki populasi rata-rata misalnya; ini bukan ide semua orang tentang hal yang baik).

Jadi misalnya, dalam kasus lognormal, ketika Anda berekspresi kembali, Anda memiliki perkiraan , dan Anda mungkin mencatat bahwa rata-rata populasi adalah , jadi Anda mungkin berpikir untuk meningkatkan dengan menskalakannya dengan beberapa perkiraan . $\exp(\mu_i)$ $\exp(\mu_i+\frac{1}{2}\sigma^2)$ $\exp(\hat{\mu_i})$ $\exp(\frac{1}{2}\sigma^2)$

Orang setidaknya harus bisa mendapatkan estimasi yang konsisten dan memang beberapa asimtotik distribusi melalui teorema Slutsky (khususnya bentuk produk) selama seseorang dapat secara konsisten memperkirakan penyesuaian. Teorema pemetaan berkelanjutan mengatakan bahwa Anda dapat jika Anda dapat memperkirakan secara konsisten ... yang merupakan kasusnya. $\sigma^2$

Jadi selama adalah penaksir yang konsisten dari , maka konvergen dalam distribusi ke distribusi (yang dengan inspeksi kemudian akan didistribusikan secara asimptotik secara lognormal. ). Karena akan konsisten untuk , bu teorema pemetaan berkelanjutan, akan konsisten untuk , dan karenanya kami memiliki penduga yang konsisten dari berarti pada skala aslinya. $\hat{\sigma}^2$ $\sigma^2$ $\exp(\hat{\mu_i})\cdot \exp(\frac{1}{2}\hat{\sigma}^2)$ $\exp(\hat{\mu_i})\cdot \exp(\frac{1}{2}\sigma^2)$ $\hat{\mu_i}$ $\mu_i$ $\exp(\hat{\mu_i})$ $\exp(\mu_i)$

Lihat di sini .

Beberapa pos terkait:

Kembali transformasi model MLR

Kembali Transformasi

Interval kepercayaan kembali ditransformasikan

Glen_b -Reinstate Monica
sumber

Terima kasih, saya melihat posting sebelumnya dan, sementara mencerahkan, masih agak bingung, maka pertanyaan saya.

Glen

+1 Jawaban bagus! Hanya klarifikasi cepat: Dari mana berasal sebagai scaler untuk ? Saya melihatnya dalam definisi lognormal di Wikipedia tetapi tidak dijelaskan di sana juga, apakah itu hanya mengintegrasikan maksud dari PDF?

\frac{1}{2}

$\frac{1}{2}$

\hat{σ^{2}}

$\hat{\sigma^2}$

usεr11852

Anda harus bisa mendapatkannya hanya dengan langsung mengintegrasikan: mana adalah kepadatan untuk lognormal, tetapi mungkin lebih mudah dilakukan dengan menghitung untuk normal (di mana ), tetapi mungkin lebih baik untuk menemukan MGF untuk - yang tidak lebih sulit - dan dari saat-saat untuk sangat mudah diperoleh (dengan mengganti oleh pada gilirannya), pada dasarnya mendapatkan momen yang lebih tinggi secara gratis.

E (Y) = \int_{0}^{\infty} y f (y) d y

$E(Y) = \int_0^\infty y\, f(y)\, dy$

f

$f$

E (e^{X})

$E(e^X)$

X = \log Y

$X=\log Y$

X

$X$

Y

$Y$

t

$t$

1, 2, . . .

$1,2,...$

Glen_b -Reinstate Monica

e^{x}

$e^x$

e^{t x}

$e^{tx}$

e^{. . .}

$e^{...}$

x

$x$

\frac{1}{2}

$\frac12$

t

$t$

e^{μ t + \frac{1}{2} σ^{2} t^{2}}

$e^{\mu t + \frac12 \sigma^2t^2}$

Kembali mentransformasikan hasil regresi saat memodelkan log (y)

Jawaban: