Very useful, but I had two questions. (1) LogRank gives an error at line 116, at [table1 table12 t1 T1 xcg1 ycg1 lambda1]=kmplot(x1,0.05,cflag,0);
This is fixed by deleting the last arg to kmplot, giving kmplot(x1,0.05,flag);
2) The p-value for LogRank is twice what it should be? 0.03058 vs. 0.01529 in the file header. Perhaps, mismatch between old documentation and changed function, or extra factor of 2 in p = 2*(1-0.5*erfc(-z/realsqrt(2))); %p-value
Removing Yates' gives 0.019250, which is correct. I tend to think 0.03058 is correct, with the Yates' correction. Thanks for the super-useful code!
Great Job ! It;s been really helpful ! thanks !
I am trying to modify the code for use in 3D . I have a problem understanding the function fastHessian_getResponse . Could you please give a little description on it so i would know what should be expected from this part so I would know how to change it for 3D ?