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
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
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
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.
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
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
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
7 matches
Mail list logo