[R] Cross-over Data with Kenward-Roger correction

2015-06-07 Thread knouri
Dear all:for the folowing data, a two-period, two treatment (A=1 vs. B=2) cross-over is fitted using the folowing SAS code.  data one; input  Sbj  Seq  Per  Trt  PEF; cards; 1  1   1   1   310 1  1   2   2   270 4  1   1   1   310 4  1   2   2   260 6  1   1   1   370 6  1

[R] MCMC with cumulative link models

2013-07-30 Thread knouri
Hi: Could anyone please let me know where I can find [R] code for implementing MCMC with cumulative link models (i.e. for analysis of ordinal probit or ordinal logistic models).   Thanks in advance, Regards, Keramat Nourijelyani, PhD Associate Professorof Biosta

[R] Bivariate skew normal cdf; very slow

2013-08-26 Thread knouri
Dear all, I am calculating the bivariate skew normal cdf in "sn" package using "pmsn" function. Although it is quite convenient ( thanks to prof. Azzalini) but it seems to be slow. For example, it takes about 1 minute in calculation of 100k of such cdf values. I am thinking to write a c++ code fo

[R] MLE for probit regression. How to avoid p=1 or p=0

2013-05-27 Thread knouri
Dear all: I am writing the following small function for a probit likelihood. As indicated, in order to avoid p=1 or p=0, I defined some precisions. I feel however, that there might be a better way to do this. Any help is greatly appreciated.

[R] Bayesian Haplotype analysisin R

2012-12-11 Thread knouri
Dear members: Could you please tell me the package for Bayesian haplotype estimation  in R. Currently PHASE is available for such purpose but I need a similar R package. Best regards, Keramat Nourijelyani, PhD Associate Professorof Biostatistics Tehran University

[R] problem with the MLE of the skew-normal shape parameter

2012-03-29 Thread knouri
Hi everyone: my likelihood function involves the shape parameter of the skew-normal in addition to other parameters. I used both optim and nlm function to find the MLS's. However, every time that I use a different initial value for the shape parameter, I get a different estimate for it. No such p

[R] clogit & weights

2015-12-21 Thread knouri via R-help
Merry Christmas everyone: I have the following data(mydat) and would like to fit a conditional logistic regression model considering "weights". id  case exposure   weights 1 1 1  2 1 0 0  2 2 1 1  2 2 0 0  2 3 1