This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GCTA make GRM error contains illegal chr number

Hi

I created bim, fam, and bed files by Plink with the use of these codes.

> shell("plink --file r --maf 0.05 --geno 0.01 --mind 0.1 --hwe 1e-6
> --nonfounders --horse --make-bed --out r") 
> 
> shell("plink --file r --maf 0.05 --geno 0.01 --mind 0.1 --hwe 1e-6
> --nonfounders --horse --recode --out r") 
> 
> 
> shell("plink --file r --maf 0.05 --geno 0.01 --mind 0.1 --hwe 1e-6
> --nonfounders --horse --recodeA --out rA") 
> 
> shell("plink --file r --maf 0.05 --geno 0.01 --mind 0.1 --hwe 1e-6
> --nonfounders --horse --recodeAD --out rAD")

Now, I want to create GRM by GCTA software.

> gcta64  --bfile r  --autosome  --make-grm  --out r

But I am faced with an error.

Error: Line 1015 of [file.bim] contains illegal chr number, please check
Line of 1015:   27  M2609   0   913820  1   2

I removed this line but I see again this error in another line.

How can I solve this problem?

Thank you for your time

plink gcta
• 4,627 views • link
written by Star:) • Star:)

Maybe find all the rows with "illegal chr number" and fix them.

• 0 views • link
written by zx8754 zx8754

2 answers

By default GCTA assumes human genome (22 pairs of autosomal chromosomes). For your horse genome, set the autosome-num flag to the appropriate number.

From https://cnsgenomics.com/software/gcta/#Datamanagement :

--autosome-num 22
Specify the number of autosomes for a species other than human. For example, if you specify the number of autosomes to be 19, then chromosomes 1 to 19 will be recognized as autosomes and chromosome 20 will be recognized as the X chromosome. The default number is 22 if this option not specified.

• 0 views • link
written by Lemire Lemire

Thank you for your help. I did that.

I put --autosome-num 29 for making GRM and MLMA.

• 0 views • link
written by Star:) • Star:)

Log in to answer this question.

More posts like this