> x <- readLines(textConnection("Testing marker:
s_1
+ ---------------------------------------------
+
+ Allele   df(0) -LnLk(0)   df(T) -LnLk(T)   ChiSq       p
+     3 7995 29320.30 7994 29311.85 16.90 4e-05 (2229/8000 probands)
+
+ Testing marker:                           s_2
+ ---------------------------------------------
+
+ Allele   df(0) -LnLk(0)   df(T) -LnLk(T)   ChiSq       p
+     3 7995 29339.25 7994 29338.31 1.88 0.1702 ( 121/8000 probands)
+
+ Testing marker:                           s_3
+ ---------------------------------------------
+
+ Allele   df(0) -LnLk(0)   df(T) -LnLk(T)   ChiSq       p
+     3 7995 29340.83 7994 29338.97 3.71 0.0540 (2429/8000 probands)"))
> closeAllConnections()
>
> # find the markers
> markers <- grep("^Testing marker", x)
> # get the marker number
> nums <- sub(".*:\\s+(\\S <file://s/>+).*", "\\1 <file://0.0.0.1/>",
x[markers])
> # get the 'p' (7th value on +4 line)
> pvalues_list <- strsplit(x[markers + 4], "\\s <file://s/>+")
> # separate out the 7th value (really 8th since first is blank)
> pvalues <- sapply(pvalues_list, '[', 8)
> # assume data structured well
> names(pvalues) <- nums
> pvalues
     s_1      s_2      s_3
 "4e-05" "0.1702" "0.0540"
>


On Fri, Mar 26, 2010 at 3:14 PM, <lixxx554@umn.edu> wrote:

> Dear R listers,
>
> I have a data file looks like the following:
>
> Testing marker:                           s_1
> ---------------------------------------------
>
> Allele   df(0) -LnLk(0)   df(T) -LnLk(T)   ChiSq       p
>     3 7995 29320.30 7994 29311.85 16.90 4e-05 (2229/8000 probands)
>
> Testing marker:                           s_2
> ---------------------------------------------
>
> Allele   df(0) -LnLk(0)   df(T) -LnLk(T)   ChiSq       p
>     3 7995 29339.25 7994 29338.31 1.88 0.1702 ( 121/8000 probands)
>
> Testing marker:                           s_3
> ---------------------------------------------
>
> Allele   df(0) -LnLk(0)   df(T) -LnLk(T)   ChiSq       p
>     3 7995 29340.83 7994 29338.97 3.71 0.0540 (2429/8000 probands)
>
> There are 10,000 testing markers, s_1 to s_10000. I want to read in all the
> 'p' values of the testing markers from this file. How can I do this?
>
> Thank you!
>
> xiang
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
> and provide commented, minimal, self-contained, reproducible code.
>



-- 
Jim Holtman
Cincinnati, OH
+1 513 646 9390

What is the problem that you are trying to solve?

	[[alternative HTML version deleted]]

