Hi everyone!
I have 170 samples. For each sample I have a vcf file with the SNPs and an other file with values corresponding to Drug sensitivity values for 500 drugs.
The final files are:
File2:
FID IID/Drug Drug1 Drug2 Drug3. ..... Drug500 0 Sample1 0.8 0.4 ........ ...... 0.9 0 Sample2 0.3 -9 ........... ........... 0.3 . . . 0 Sample170
and the merged vcf file for all samples
FIle1 (vcf):
Chr Pos ... Sample1 Sample2 ... Sample170 1 1 1 . . . . . Y
Can you help me on how to make association between specific SNPs and the abberations (changes) in drug sensitivity of a particular drug?
Do I have to use a command like this? Is this right? Have you other suggestions?
plink --vcf-half-call r --const-fid 0 --vcf MERGED.vcf.gz --pheno DrugSensitivityAnalysisFinal.tsv --allow-no-sex --all-pheno --linear --pfilter 5e-4 --ouT results