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@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
______________________________________________
R-help@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.
