Displaying 1 result from an estimated 1 matches for "tp1559382p1560118".
2010 Feb 17
2
Trouble with optim function
Hi all,
I'm trying to make a little script to determine an "unknown" rate for a
number of known exponential trials.
My Code:
#Set Trials and generate number
trials=100
rand<-runif(1,0,1)
vector=0
#Generate vector of 100 random exponentials and sum them
for (i in 1:100) {
vector<-rexp(trials,rate=rand)
}
sumvect=sum(vector)
#Create the log likelihood function