An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120205/00abdf9f/attachment.pl>
R- Fisher Information
2 messages · Christopher Kelvin, David Winsemius
On Feb 5, 2012, at 1:20 PM, Christopher Kelvin wrote:
Dear All, Can you help me, with the code below how do I obtain the fisher information from it.
A modest amount of searching the archives (which is really your responsibility): search.r-project.org/cgi-bin/namazu.cgi?query="fisher + information "+ optim &max = 100 &result = normal &sort =score&idxname=functions&idxname=Rhelp08&idxname=Rhelp10&idxname=Rhelp02 ... produces: http://finzi.psych.upenn.edu/R/Rhelp02/archive/72856.html You should go through the rest of that thread because there is other usefull inofrmation in it.
Is my q<-replicate(1000,x) the right way to do simulation. thank you. x<-rweibull(100,0.8,1.5) q<-replicate(1000,x)
I doubt that is what you want. Did you even look at what it produced? > x<-rweibull(100,0.8,1.5) > q<-replicate(2,x) > str(q) num [1:100, 1:2] 2.1528 1.8345 5.3473 0.3728 0.0434 ... > q[1,] [1] 2.152798 2.152798
David
> z<-function(p){
> beta<-p[1]
> eta<-p[2]
> log1<-(n*log(beta)-n*beta*log(eta)+(beta-1)*sum(log(x))-sum((x/
> eta)^beta))
> return(-log1)
> }
> zz<-optim(c(0.5,0.5),z)
> zz
>
>
>
> Chris Guure
> postgraduate researcher/tutor
> Institute for Mathematical Research
> Universiti Putra Malaysia
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
David Winsemius, MD
West Hartford, CT