Sponsored Content
Top Forums Shell Programming and Scripting Problem with output file plink GWAS Post 302905631 by boro82 on Thursday 12th of June 2014 02:10:10 PM
Old 06-12-2014
Problem with output file plink GWAS

hi,

Can anybody help me about the odd output file yielded by the following command:

Code:
plink --noweb --allow-no-sex --bfile input_file --assoc --adjust --out output_file

In fact, I'd like to perform a simple Bonferroni basic association test, and I obtains the following output file (extract for first ten rows):

Code:
 CHR         SNP      UNADJ         GC       BONF       HOLM   SIDAK_SS   SIDAK_SD     FDR_BH     FDR_BY
   1   rs7540001        INF        INF        INF        INF        INF        INF        INF        INF 
   1   rs2743476        INF        INF        INF        INF        INF        INF        INF        INF 
   1   rs9435653        INF        INF        INF        INF        INF        INF        INF        INF 
   1  rs11810307        INF        INF        INF        INF        INF        INF        INF        INF 
   1  rs12410470        INF        INF        INF        INF        INF        INF        INF        INF 
   1   rs2873038        INF        INF        INF        INF        INF        INF        INF        INF 
   1   rs1553451        INF        INF        INF        INF        INF        INF        INF        INF 
   1   rs2580567        INF        INF        INF        INF        INF        INF        INF        INF 
   1  rs11578605        INF        INF        INF        INF        INF        INF        INF        INF 
   ....
   ....

This file contains (528000 rows (SNPs)), in particular: 328 rows with BONF = "INF", and 1300 others rows with BONF < 0.05, the remaining SNPs present a higher BONF-value than 0.05.
My question is: 1300+328 = 1628 (huge) is really the number of SNPs declared significants (<0.05) ? is it normal to get this kind of file ?

Thanks for you
Moderator's Comments:
Mod Comment Please use CODE tags when displaying sample input, output, and code.

Last edited by Don Cragun; 06-12-2014 at 03:18 PM.. Reason: Change BOLD tags to CODE tags and add CODE tags for sample output.
 

10 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

Plink problem....only works if passwd is in the script

I have a one line bat script run off a XP machine that tar's and compresses some files from a Sol 8 box. It goes something like this (a bit simplified).... plink -pw <passwd> user@host "tar -cvf - -C / tmp/a_file | compress " > a_file.tar.Z So this works....and it's worked many times. But now... (3 Replies)
Discussion started by: Yinzer955i
3 Replies

2. Shell Programming and Scripting

problem parsing output file

Hi Gurus, I am using the following code to parse the output of a file. This code basically parses the file and adds | at the end of each field. But I am getting it wrong in some cases. I have explained it below #----- parse output file, get each field position, ----- #----- and construct... (13 Replies)
Discussion started by: ragha81
13 Replies

3. Shell Programming and Scripting

How we can use plink?

Hi, How we can use use plink to access unix system using Dos. Could someone send me the commands that can be use in Batch file to call unix system using plink utility. Thanks in advance (0 Replies)
Discussion started by: shekhar_ssm
0 Replies

4. Shell Programming and Scripting

Error with Plink

Hello. I have a TCL script that logs in to a server using SSH. As SSH isn't available in windows,I used Plink to do the job.The script works fine on my PC and 2 of my friend's PC. However, on one PC, I get the following error message: "'D:\scripts\plink.exe' is not a Win32 Console... (1 Reply)
Discussion started by: plasmalightwave
1 Replies

5. Shell Programming and Scripting

problem with suppressed output to file using echo and tee command

Hi, When I run the following command in terminal it works. The string TEST is appended to a file silently. echo TEST | tee -a file.txt &>/dev/null However, when I paste this same line to a file, say shell1.sh, and use bourne shell . I run this file in terminal, ./shell1.sh. However I... (1 Reply)
Discussion started by: shahanali
1 Replies

6. Shell Programming and Scripting

problem with print append to output file syntax

I'm trying to output the contents of the infile to the outfile using Append. I will want to use append but the syntax doesn't seem to be working ! Input file (called a.txt) contains this: a a a b b b I'm running shell script (called k.sh) from Unix command-line like this: ./k.sh .... (1 Reply)
Discussion started by: script_op2a
1 Replies

7. Shell Programming and Scripting

problem with output re director to a file

I have a file named XXX_TEF_20120530234576.txt.pgp I am creating a new file with the below syntax but looks it doesnt like the redirection for file in XXX_TEF_20120530234576.txt.pgp ;do; echo $file cat $file | /usr/bin/gpg/ -d > echo $file | cut -f1,2 -d'.' done; What I want actually is... (6 Replies)
Discussion started by: dsravan
6 Replies

8. Shell Programming and Scripting

PLINK help

Hi Plink users, I am very new in GWAS and decided to use PLINK for this. But I am struggling with file formats. I have one exome data with 200000 SNPs in it. But it is in .txt format. But all the necessary fields are there in different columns. Can you please suggest me how I can create the... (5 Replies)
Discussion started by: smitra
5 Replies

9. Shell Programming and Scripting

Need help on Plink

Hi All, Iam a newbie to the plink and need your assistance. I have referred some posts but it doesn't helps me much. I have two steps to do. 1. I have a config file which has a list of servers,username and password. 2. I have a shell script in windows which accepts arguments and need to... (0 Replies)
Discussion started by: cskumar
0 Replies

10. Windows & DOS: Issues & Discussions

Plink wait problem

Hi, I have run into a problem to which i can't seem to find any solution, posting here is my last resort. Problem: I am using plink to access my router and run a few configuration commands. When in enter configurations mode, instead of sending next command plink keeps on waiting for manual... (7 Replies)
Discussion started by: zaainabbas
7 Replies
ppmtosixel(1)                                                 General Commands Manual                                                ppmtosixel(1)

NAME
ppmtosixel - convert a portable pixmap into DEC sixel format SYNOPSIS
ppmtosixel [-raw] [-margin] [ppmfile] DESCRIPTION
Reads a portable pixmap as input. Produces sixel commands (SIX) as output. The output is formatted for color printing, e.g. for a DEC LJ250 color inkjet printer. If RGB values from the PPM file do not have maxval=100, the RGB values are rescaled. A printer control header and a color assignment table begin the SIX file. Image data is written in a compressed format by default. A printer control footer ends the image file. OPTIONS
-raw If specified, each pixel will be explicitly described in the image file. If -raw is not specified, output will default to com- pressed format in which identical adjacent pixels are replaced by "repeat pixel" commands. A raw file is often an order of magni- tude larger than a compressed file and prints much slower. -margin If -margin is not specified, the image will be start at the left margin (of the window, paper, or whatever). If -margin is speci- fied, a 1.5 inch left margin will offset the image. PRINTING
Generally, sixel files must reach the printer unfiltered. Use the lpr -x option or cat filename > /dev/tty0?. BUGS
Upon rescaling, truncation of the least significant bits of RGB values may result in poor color conversion. If the original PPM maxval was greater than 100, rescaling also reduces the image depth. While the actual RGB values from the ppm file are more or less retained, the color palette of the LJ250 may not match the colors on your screen. This seems to be a printer limitation. SEE ALSO
ppm(5) AUTHOR
Copyright (C) 1991 by Rick Vinci. 26 April 1991 ppmtosixel(1)
All times are GMT -4. The time now is 08:15 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy