作者k84526991 (哇哈哈)
看板Statistics
标题[问题] R MLE问题和设虚拟变数问题
时间Thu Jan 21 23:51:38 2010
请问一下
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
1F：推 fairwind:2. D <- ifelse(y > 100, 1, 0) 01/22 00:49
2F：推 acecc:nlm;nlminb;optim 01/22 01:05