[R] Help using Cast (Text) Version

[Steve Sidney]

David

Thanks, I'll try that......but no what I need is the total (1's) for each of 
the rows, labelled 1-6 at the top of each col in the table provided.

What I guess I am not sure of is how to identify the col after the melt and 
cast.

Steve

----- Original Message ----- 
From: "David Winsemius" <dwinsemius at comcast.net>
To: "Steve Sidney" <sbsidney at mweb.co.za>
Cc: <r-help at r-project.org>
Sent: Sunday, January 17, 2010 4:39 PM
Subject: Re: [R] Help using Cast (Text) Version


>
> On Jan 17, 2010, at 5:31 AM, Steve Sidney wrote:
>
>> Sorry to repeat the meassage, not sure if the HTML version has been 
>> received - Apologies for duplication
>>
>> Dear list
>>
>> I am trying to count the no of occurances in a column of a data  frame 
>> and there is missing data identifed by NA.
>>
>> I am able to melt and cast the data correctly as well as sum the 
>> occurances using margins and sum.
>>
>> Here are the melt and cast commands
>>
>> bw = melt(res, id=c("lab","r"), "pf_zbw")
>> b = cast(bw, lab ~ r, sum, margins = T)
>>
>> Sample Data (before using sum and margins)
>>
>>    lab  1  2  3  4  5  6
>> 1  4er66  1 NA  1  0 NA  0
>> 2  4gcyi  0  0  1  0  0  0
>> 3  5d3hh  0  0  0 NA  0  0
>> 4  5d3wt  0  0  0  0  0  0
>> .
>> . lines deleted to save space
>> .
>> 69 v3st5 NA NA  1 NA NA NA
>> 70 a22g5 NA  0 NA NA NA NA
>> 71 b5dd3 NA  0 NA NA NA NA
>> 72 g44d2 NA  0 NA NA NA NA
>>
>> Data after using sum and margins
>>
>>    lab 1 2 3 4 5 6 (all)
>> 1  4er66 1 0 1 0 0 0     2
>> 2  4gcyi 0 0 1 0 0 0     1
>> 3  5d3hh 0 0 0 0 0 0     0
>> 4  5d3wt 0 0 0 0 0 0     0
>> 5  6n44r 0 0 0 0 0 0     0
>> .
>> .lines deleted to save space
>> .
>> 70 a22g5 0 0 0 0 0 0     0
>> 71 b5dd3 0 0 0 0 0 0     0
>> 72 g44d2 0 0 0 0 0 0     0
>> 73 (all) 5 2 4 3 5 7    26
>>
>> Uisng length just tells me how many total rows there are.
>
>
>> What I need to do is count how many rows there is valid data, in  this 
>> case either a one (1) or a zero (0) in b
>
> I'm guessing that you mean to apply that test to the column in b  labeled 
> "(all)" . If that's the case, then something like (obviously  untested):
>
> sum( b$'(all)' == 1 | b$'(all)' == 0)
>
>
>
>>
>> I have a report to construct for tomorrow Mon so any help would be 
>> appreciated
>>
>> Regards
>> Steve
>
> David Winsemius, MD
> Heritage Laboratories
> West Hartford, CT
>
>