I am using R version 1.8.1. I am doing following tast:
First generate 100 Gaussion(3,1) numbers, then write the likelihood function to estimate the parameters of Gaussian distribution by direct maximizing the likelihood function.
My likelihood function is:
fn<-function(x) (-50*log((sd(x))^2))-50*log(sqrt(2*pi))-(1/2*((mean(x))^2))*(sum((x-(mean(x))^2))
After I input above function, the '' >'' turned to be '' + ''. and it did not work.
I am looking for the help to solve this tast by writting a likelihood function.
Thanks a lot.
Best regards. edward sun p.s
x=rnorm(100, mean=3, sd=1)
library(MASS) fitdistr(x, "normal")
and
fitdistr(x, "normal", list(mean=0, sd=1))
do not work.
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
