Displaying 1 result from an estimated 1 matches for "fntr1".
Did you mean:
intr1
2005 Aug 02
1
Loop problem
...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$e...