Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: permute xtmixed


From   Fotios Drenos <rmhafdr@ucl.ac.uk>
To   statalist@hsphsun2.harvard.edu
Subject   st: permute xtmixed
Date   Tue, 23 Oct 2007 17:30:59 +0100

Dear All,
I am trying to run a permutation tests on an xtmixed command and I get a number of errors.

I am using the "use http://www.stata-press.com/data/r9/pig"; dataset:

The program for the xtmixed command is:

program test7, rclass
version 9
args dependent fac_week
xi: xtmixed `dependent' i.`fac_week' ||id: `fac_week', cov(unstruct)
assert e(converged)==1
testparm _I`fac_week'*
return scalar pgene=r(p)
end

and then using 0 , 1 , 2 or 3 as seeds I run the command:
permute weight t=r(pgene), reps(2) noi: test7 weight week

For seeds 0 1 and 2 problems arise during the permutation, while seed 3 seems to work.

I have used the assert e(converged)==1 after the xtmixed and fails for the first three seeds.

Is there something that I fail to see? Is there a specific reason while the xtmixed command is not converging?

All the best

*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/




© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index