Print an array without the duplicate

2004-11-18 Thread Li, Aiguo (NIH/NCI)
Hi, all. I need to print a header line without duplicated item and am struglling with the elimination of replicated items. The following codes are nor working yet. Could anyone help me to make this work? #!usr/bin/perl use strict; use warnings; my @header; my $i ; my $size = @header; while

Where is it wrong with my code

2004-12-02 Thread Li, Aiguo (NIH/NCI)
Hello, all. I am trying to assign a "P" for any values greater than 1.0 and assign a "A" otherwise. However, I need to skip the header line and the first column. Something is wrong with my code and it does not skip the first column well. Please help me to detect the bug. Thanks, Aiguo #!/usr/b

RE: Where is it wrong with my code

2004-12-02 Thread Li, Aiguo (NIH/NCI)
Never mind. It is working now. Thanks, Aiguo -Original Message- From: Li, Aiguo (NIH/NCI) Sent: Thursday, December 02, 2004 10:11 AM To: Perl Beginners Subject: Where is it wrong with my code Hello, all. I am trying to assign a "P" for any values greater than 1.0 and a

Merge two files

2004-12-07 Thread Li, Aiguo (NIH/NCI)
Hello, all. I need to merge two files according a id list in each file. What I am trying to do is to have a sub routine that get the id from data file 1 as key for a hash, a second sub routine to put the value from data file 1 into the hash, then while reading in the data file 2 using the while

Why the if loop does not work??

2004-12-09 Thread Li, Aiguo (NIH/NCI)
I am trying to merge two file based on a SNP_A-## list in each file. For some reason, the regular expression in the if loop does not work and I can not match the key values generated from hash to the string from the input file. Could anybody help me detect the problem? Thanks, AG #!usr/lo

How to go back to a line and do something with it in an input fil e?

2005-01-24 Thread Li, Aiguo (NIH/NCI)
hello, all. I need to transform the column3 data from an input file as indicated below into three values: 0, 0.4, 1. if there are more than 3 NoCall continuously in column3, change NoCall into "1". this option has the highest precedence. for the rows where column2 is between 52 and 105 (range from

Number or string?

2005-02-22 Thread Li, Aiguo (NIH/NCI)
Hi, all. I have the following data from a file. __data from file IN__ SNP_A-1512608 23 148840899 0.8281090.823391 11128 SNP_A-1512550 23 148841154 1.7173971.750476 11129 SNP_A-1518843 23 149078514 0.8322850.99744 11130 SNP_A-1507809

question related to readdir function

2005-07-15 Thread Li, Aiguo (NIH/NCI)
Dear all. I wrote a piece of code to read and open all files from a directory to do something with it. But I found that the variable assigned to @file_array split the file names into pieces. For example ". . test.rpt" file will be separated into three elements in an array, which are "." "." "tes

A question related to column subsetting

2005-08-19 Thread Li, Aiguo (NIH/NCI)
Hi, all. I need to write a perl script to extract a subset of columns based on the column header. For example, I have a dataset containing 200 columns and a list of 10 column names. I like to extract a dataset that contains the 10 columns of data only from the whole dataset. What in my mi

Questions related to draw graphy using perl

2004-10-15 Thread Li, Aiguo (NIH/NCI)
Dear all. I am a new user of Perl with some experience in perl data extraction with pattern matching, but never did anything with making graphy using perl. I need to draw a chromosome copy number graphy and p-values. The data looks are as follow: SNP id physical location copy number

Generate a perl list file

2004-10-19 Thread Li, Aiguo (NIH/NCI)
Dear all. I have a tab delimited file as follow: V namep 1.0 AAA 0.001 0.9 BBB 0.003 0.8 CCC 0.004 . I need to convert the file into following format: { labels = ( {v="1.0"; name = "AAA"; p = "0.001"; }, {v="0.9

P/A ratio

2004-11-09 Thread Li, Aiguo (NIH/NCI)
Hi, all. I have the following dataset and need to calculate a P/A ratio for each replicates in the dataset. In this case, treatment 1 has 4 replicats and treatment2 has 3 replicates. The P/A = [((#of P)*2) + (# of M)]/# of replicates. The output should be two columns of P/A ratios for two treatme

Why doesn't this work?

2004-11-09 Thread Li, Aiguo (NIH/NCI)
Hello, I have the following dataset and want to calculate a P/A ratio for each replicates in the dataset. In this case, treatment 1 has 4 replicats and treatment2 has 3 replicates. The P/A = [((#of P)*2) + (# of M)]/# of replicates. The output should be two columns of P/A ratios for two treatment

RE: Why doesn't this work?

2004-11-11 Thread Li, Aiguo (NIH/NCI)
Hi, Zeus. Thanks for your comment. "Replicates" line is the header line and there are two treatments in this case. Probe idTreat1 Treat1 Treat1 treat1 Treat2 treat2 treat2 AFFX-BioB-5_at P P P P P P P AFFX-BioB-M_at P P M P P

a question related to working with array

2008-08-15 Thread Li, Aiguo (NIH/NCI) [E]
Hello, I have a programming questions to you all, I have three arrays: arrayA1, arrayA2, (arrayA1 and arrayA2 are same length) and array B. Using the following code I created the fourth arrayC that contain intersection elements between arrayA1 and arrayB. Now I need to get the corresponding

A questions about how to split files

2009-01-15 Thread Li, Aiguo (NIH/NCI) [E]
H all, I need to split a file containing three columns of data as shown below into three separate files. Each split file should contain row names and one column of data and the column name should be the file names. Is there any perl advanced function that allow me to do this? probe set E

a question related to file selection

2009-01-23 Thread Li, Aiguo (NIH/NCI) [E]
Hi all, I need copy files from a directory daily to a folder. How can I select files based on dates that were created? Thanks, AG Lee