Hi,
I have got problems integrating the following function using "integrate":
lambdat<-function(t){
tempT<-T[k,][!is.na(T[k,])]#available values from k-th row of matrix T
tempJ<-J[k,][!is.na(J[k,])]
hg<-length(tempT[tempT<=t & tempJ==0])#counts observations satisfing the
conditions
ag<-length(tempT[tempT<=t & tempJ==1])
lambdaXY[hg+1,ag+1]#takes values from a 10x10 matrix
}
I keep receiving this message:
1: longer object length
is not a multiple of shorter object length in: tempT <= t
2: longer object length
is not a multiple of shorter object length in: tempT <= t & tempJ == 0
What I suspect is that the "integrate" function submits the whole vector of
points at which the integral is to be evaluated at once. For my function to be
integrated it would rather have to be evaluated at each point after another in
a loop of some kind.
Can you think of a way to solve this problem without me having to write the
integrating procedure from scratch (I have no idea about FORTRAN and this is
what the "integrate" description refers to)?
Thank you.
Kind Regards,
Lukasz Szczepanski
Student
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
