From: Hathaikan Chootrakool <Hathaikan.Chootrakool_at_newcastle.ac.uk>

Date: Wed 03 Aug 2005 - 05:36:23 EST

(maximum likelyhood function)

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Aug 03 05:44:08 2005

Date: Wed 03 Aug 2005 - 05:36:23 EST

Dear everyone

I am a new user,would like to combine these code together by using a loop,each function has three value as Tr = 1 - 3,how can i combine together?

logitTr1 <-logit[logit[,"Study"]&logit[,"Tr"]==1,]

(number of row in each group (1-3) is difference but equal in colume)

fnTr1 <- function (p) sum(

n/2*log(2*pi)+log(1/logitTr1$sd)+1/2*(logitTr1$logitp*logitTr1$logitp-2*logitTr1$logitp*p+p^2)

*1/logitTr1$sd*logitTr1$sd )

(maximum likelyhood function)

outTr1<- nlm (fnTr1,p=c(10),hessian=TRUE)
minimumTr1 <- outTr1$minimum

valueTr1 <- outTr1$estimate

(estimate the value)

The problem is the program couldn't work by using logitTr[i],fnTr[i] outTr[i],minimumTr[i],value[i] in a loop. The function logitTr[i] is the data matrix which is not equal each group, fnTr[i] is the maximumlikelyhood function for estimation in next step; outTr[i],minimumTr[i],valueTr[i].

Has anyone got any idea to help me?,thank you very much.

Hathaikan

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Aug 03 05:44:08 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:39:39 EST
*