Displaying 1 result from an estimated 1 matches for "minimumtr1".
Did you mean:
minimumtr
2005 Aug 02
1
Loop problem
...h 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 fun...