lincomest and STATA 15.1

[laure rouch]

Dear STATA users,

I had written a few months ago a program that I was using with a previous version of STATA.

The program is as follows:

************************************************** ****************
************************************************** ****************

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

tempfile tf1 tf2 tf3 tf4 tf5 tf6

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

** 1 vs 0
lincomest 1.group_96+1.group_96#2006.year
parmest,label saving(tf1,replace)

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

** 2 vs 0
lincomest 2.group_96+2.group_96#2006.year
parmest,label saving(tf2,replace)

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

** 3 vs 0
lincomest 3.group_96+3.group_96#2006.year
parmest,label saving(tf3,replace)

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

** 3 vs 1
lincomest 3.group_96+3.group_96#2006.year-1.group_96-1.group_96#2006.year
parmest,label saving(tf4,replace)

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

** 3 vs 2
lincomest 3.group_96+3.group_96#2006.year-2.group_96-2.group_96#2006.year
parmest,label saving(tf5,replace)

xi:xtmixed perfgen group_96##year i.age96 sexe etude3 sedentarite obesite tabagisme diabete dyslipidemie alcoolisme pathocv || numsujet : i.year

** 2 vs 1
lincomest 2.group_96+2.group_96#2006.year-1.group_96-1.group_96#2006.year
parmest,label saving(tf6,replace)

dsconcat tf1 tf2 tf3 tf4 tf5 tf6

qqvalue p, method(hochberg) qvalue(myqval)
list

************************************************** ****************
************************************************** ****************

When I want to run it with STATA 15.1 it does not work, especially the command lincomest.

Could you please help with to solve this problem? How should I write the program to make able to be used with STATA 15.1 ?

Thank you so much in advance for your help +++

Best,

Laure

[Richard Williams]

What exactly do you mean by "does not work?" Are you getting error messages and if so what are they? Are you getting results, but you think they are wrong? Or what?

If you were using, say Stata 12 before, try making your first command

version 12

lincomest and parmest are user-written commands. Do you have them installed?

Finally, note that you are using some outdated commands here, e.g. xi:, xtmixed. If they still do what you want, fine, but you may want to use newer commands.

[laure rouch]

Actually, lincomest was unrecognized but indeed I think I had not installed the package.

Thank you very much Richard for your help,

Best,

Laure