推 fairwind:2. D <- ifelse(y > 100, 1, 0) 01/22 00:49
推 acecc:nlm;nlminb;optim 01/22 01:05
請問一下
1.我要跑一個常態的MLE，要估計平拘數和變異數
如果我有一組x的資料，我的概似函數可以寫成
fn=function(sigma,mu)
(n/2)*log(2*pi*sigma*sigma)+(1/(2*sigma*sigma))*t(x-mu)%*%(x-mu)
sigma是要估的變異數，mu是要估計的平均數
那有了mle的函數之後，我要用什麼指令去求出最大概似估計法的變異數估計值和
平均數估計值呢??
2.假設y是一組資料，假設是n*1的矩陣，如果我想要設dummy variable
例如 D=1 y>100
D=0 y<100
那我的語法要怎麼寫才能讓y的資料變成我要的dummy資料呢??
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.214.84