[R] compare and count data

arun smartpink111 at yahoo.com
Mon Jan 21 15:34:59 CET 2013


Hi Roslina,
No problem.
If you take out the i from the code,
1-(res/1000)
#   V1    V3    V2    V4 
#1.000 1.000 0.981 1.000 
Is it what you were looking for?

Also, in addition, you could also try:
  #if the total column number is even
critical<-rep(c(1.61,75.89),ncol(mydata)/2)
#if it is odd,
critical1<-c(rep(c(1.61,75.89),(ncol(mydata)-1)/2),1.61)

colSums(mapply("<",mydata,critical))
#V1 V2 V3 V4 
# 0 19  0  0 


A.K.







________________________________
From: Roslina Zakaria <zroslina at yahoo.com>
To: arun <smartpink111 at yahoo.com> 
Sent: Monday, January 21, 2013 2:45 AM
Subject: Re: [R] compare and count data


Hi Arun,

Thank you so much for your fast respond.  

Now I would like to use the 'res' value and calculate using the formula below:

power <- 1 - (res(i) /1000);  

power

I tried the above code but it gives me this message:

power <- function(1 - (res(i) /1000));
Error: unexpected numeric constant in "kuasa <- function(1

Thank you so much Arun for your help.


From: arun <smartpink111 at yahoo.com>
To: Roslina Zakaria <zroslina at yahoo.com> 
Cc: R help <r-help at r-project.org> 
Sent: Monday, January 21, 2013 2:12 PM
Subject: Re: [R] compare and count data


HI,

Your example data had only 4 columns
mydata<-read.table(text="
  3.014505 62.96425 3.014505 138.0673
  2.817503 56.03400 2.817503 133.3411
  2.976227 47.12192 2.976227 139.2438
  2.825495 75.05284 2.825495 129.2959
  2.793500 52.75190 2.793500 130.9874
  3.006333 54.56210 3.006333 136.2982
  3.245833 71.31358 3.245833 150.0974
  2.559980 64.28300 2.559980 114.5395
  2.548456 69.48727 2.548456 118.7383
 2.715040 58.64134 2.715040 123.0574
 2.799652 58.96281 2.799652 126.9216
 2.900839 75.43427 2.900839 131.6100
 3.183844 63.57657 3.183844 143.5530
 2.792940 63.04534 2.792940 120.1010
 3.229830 71.09790 3.229830 136.6931
 3.528115 68.64850 3.528115 156.0135
 2.973711 58.45141 2.973711 128.6738
 2.518891 53.62788 2.518891 110.6072
 2.848493 79.29142 2.848493 127.6110
 2.971362
57.81008 2.971362 129.5411
",sep="",header=F)


critical<-c(1.61,75.89)
res<-unlist(lapply(1:2,function(i) colSums(mydata[,seq(i,ncol(mydata),by=2)]<critical[i])))
 res
#V1 V3 V2 V4 
# 0  0 19  0 
A.K.


----- Original Message -----
From: Roslina Zakaria <zroslina at yahoo.com>
To: "r-help at r-project.org" <r-help at r-project.org>
Cc: 
Sent: Monday, January 21, 2013 12:39 AM
Subject: [R] compare and count data

Dear r-users,

I have these data below:

I would like to compare each column with a certain value and count how many in each column less than that specified value.  

odd column
(1,3,5) will compare with 1.61 and even column (2,4,6) will compare with 75 and I would like to count how many for each column. 


I tried these below but it give me just one set of data.

critical  <- c(1.61,75.89, 1.61, 75.89, 1.61, 75.89)   ; critical
pilih <- mydata[mydata < critical] ; pilih
jumlah <- length(pilih)   ; jumlah
kuasa <- 1 - (jumlah/1000);  kuasa

> head(mydata,20)
         V1       V1       V1       V1
1  3.014505 62.96425 3.014505 138.0673
2  2.817503 56.03400 2.817503 133.3411
3  2.976227 47.12192 2.976227 139.2438
4  2.825495 75.05284 2.825495 129.2959
5  2.793500 52.75190 2.793500 130.9874
6  3.006333 54.56210 3.006333
136.2982
7  3.245833 71.31358 3.245833 150.0974
8  2.559980 64.28300 2.559980 114.5395
9  2.548456 69.48727 2.548456 118.7383
10 2.715040 58.64134 2.715040 123.0574
11 2.799652 58.96281 2.799652 126.9216
12 2.900839 75.43427 2.900839 131.6100
13 3.183844 63.57657 3.183844 143.5530
14 2.792940 63.04534 2.792940 120.1010
15 3.229830 71.09790 3.229830 136.6931
16 3.528115 68.64850 3.528115 156.0135
17 2.973711 58.45141 2.973711 128.6738
18 2.518891 53.62788 2.518891 110.6072
19 2.848493 79.29142 2.848493 127.6110
20 2.971362 57.81008 2.971362 129.5411
    [[alternative HTML version deleted]]


______________________________________________
R-help at 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
and provide commented, minimal, self-contained, reproducible code.               



More information about the R-help mailing list