Statalist


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

Re: st: genhw error


From   David Airey <[email protected]>
To   [email protected]
Subject   Re: st: genhw error
Date   Tue, 23 Oct 2007 14:17:48 -0500

.

hwsnp is also by the Mario Cleves' and allows testing of multiple snps at once.

There is also David Clayton's fantastic set of programs for genetic analysis in Stata. HWE is tested by one of those as well, and there is helpful documentation on manipulating data so that each row is a person (genotype models) or each row is a chromosome (allele counting models); the latter assume HWE.

http://www-gene.cimr.cam.ac.uk/clayton/software/stata/

On Oct 23, 2007, at 2:08 PM, n j cox wrote:


Tony and David's experiences can be reconciled. Tony is
using an old -genhw-. It is expecting three numbers,
the frequencies of different phenotypes [? I'm a geographer],
and is only given two, and it chokes on a comma where
it expects to see a number. This is because Tony's data
is 4 AA, 2 AG -- and, importantly, 0 GG, except that last
frequency is not being passed to -genhwi-. Evidently this was
an old bug in -genhw- that was caught and fixed. Thus
Tony should update by using -findit- to find the latest
version on Mario Cleves' website.

Nick
[email protected]

David Airey

I am not seeing the same problem with genhw. I'm Stata MP OS X.

. which genhw
/Users/dairey/Library/Application Support/Stata/ado/plus/g/genhw.ado
*! version 7.0.0 14mar2005

Tony Black

I am having difficulty with the "genhw" command
With this simple dataset:

allele1 allele2
A A
A G
A A
A G
A A
A A

I issue the following command:

. genhw allele1 allele2

and get this result (in red)

',' found where integer expected

I set trace on and it seems to be blowing up in "genhwi" -- here is a snip
from the trace file

- genhwi `vars', `labopt' `binvar'
= genhwi 4 2, label(AA AG)

---------------------------------------------------------------------- ------
------------------- begin genhwi ---
- version 6
- qui {
- gettoken AA 0 : 0, parse(" ,")
- gettoken Aa 0 : 0, parse(" ,")
- gettoken aa 0 : 0, parse(" ,")
- confirm integer number `AA'
= confirm integer number 4
- confirm integer number `Aa'
= confirm integer number 2
- confirm integer number `aa'
= confirm integer number ,
',' found where integer expected


*
* 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/
--
David C. Airey, Ph.D.
Pharmacology Research Assistant Professor
Center for Human Genetics Research Member

Department of Pharmacology
School of Medicine
Vanderbilt University
Rm 8158A Bldg MR3
465 21st Avenue South
Nashville, TN 37232-8548

TEL   (615) 936-1510
FAX   (615) 936-3747
EMAIL [email protected]
URL   http://people.vanderbilt.edu/~david.c.airey/dca_cv.pdf
URL   http://www.vanderbilt.edu/pharmacology


*
*   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–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index