- 問題: 當你想要問問題時，請使用這個類別
[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
新手(沒寫過程式，R 是我的第一次)
[問題敘述]:
不知道如何從ols法 改成mle法
[程式範例]:
用OLS法去估 alpha beta時用這樣
F1=function(n)
{
A=c()
B=c()
for(i in 1:10000)
{
X=rnorm(n,0,1)
Y=rt(n,10)
L=lm(Y~X)
A[i]=L$coefficients[1]
B[i]=L$coefficients[2]
}
C=c(mean(A),var(A),mean(B),var(B))
return(C)
}
同樣的一件是若要使用MLE法去估
不知該如何輸入
[環境敘述]:
請提供 sessionInfo() 的輸出結果，
裡面含有所有你使用的作業系統、R 的版本和套件版本資訊，
讓版友更容易找出錯誤
[關鍵字]:
選擇性，也許未來有用

二者的結果在簡單迴歸不是一樣的嗎？

標準差不一樣

stats4::mle

不好意思 請問是改在哪邊?

先試著自訂一個function，這裡有例子 goo.gl/UmtsHx做好function後再拿它來套入你的code，可能比較容易

我記得MLE跟最小平方法出來的標準差只是一個scale差異而已可以直接用lm做出來之後，自己調整吧mle是1/n, least-square是用1/(n-2)mle的標準差不是不偏估計輛

有點好奇 我記得mle/ols beta^hat 的長相是一樣的http://slideplayer.com/slide/6184052/ page20剛用 lm, glm(gaussian) 取係數 跟(X'X)^X'Y 數值都一致上面少惹inverse

係數一樣沒錯，是標準差的分母差1而已。