Hi,
I've written this function:
g = function(test,p1,p2)
{
test=sort(test)
merke=0
for (z in 1:length(test))
{
F1=((2*z-1)/length(test))
F21=log(plnorm(test[z],p1,p2))
F22=log(1-plnorm(test[length(test)+1-z],p1,p2))
F2= F21+F22
merke=merke+F1*F2
}
return(-length(test)*-merke)
}
When I call it for e.g. test = 1:10 it only returns ONE value. Instead what
I want to have is 10 values (like it makes the dlnorm for instance).
I could make that with a loop, but maybe there is there a more simple trick
?
Carsten
[[alternative HTML version deleted]]