cumulative sum by group and under some criteria - R-help

arun

Wed, Feb 27, 2013 5:43 AM #

dat1<- read.table(text="
m1 n1 m n??? A????? B??? C?? D
2??? 2??? 4? 5? 0.1??? 0.2? 0.2 0.3
2?? 2???? 4? 4? 0.2?? 0.1?? 0.3? 0.4
2?? 3???? 4? 5? 0.5?? 0.6?? 0.2? 0.2
2? 4????? 3? 4? 0.2?? 0.8?? 0.5? 0.3
",sep="",header=TRUE) 


dat1[which.max(dat1$D),]
#? m1 n1 m n?? A?? B?? C?? D
#2? 2? 2 4 4 0.2 0.1 0.3 0.4
A.K.

From: Joanna Zhang <zjoanna2013 at gmail.com>
To: arun <smartpink111 at yahoo.com>
Sent: Wednesday, February 27, 2013 1:25 AM
Subject: Re: [R] cumulative sum by group and under some criteria

Thanks very much! it works!

suppose that I have a data:
m1 n1 m n??? A???? ?B?? ?C? ?D
2??? 2??? 4? 5? 0.1??? 0.2? 0.2 0.3
2?? 2???? 4??4? 0.2?? 0.1?? 0.3? 0.4
I want to identify the max value of D and extract the row , in this example, it is the 2nd row.


On Tue, Feb 26, 2013 at 11:13 PM, arun <smartpink111 at yahoo.com> wrote:

res2<- join(res1,d3,by=c("m1","n1"),type="inner")
>
>p0L<-0.05
>p0H<-0.05
>p1L<-0.20
>p1H<-0.20
>?
>res2<- within(res2,{p1<- x/m; p2<- y/n;term2_p0<-dbinom(x1,m1, p0L, log=FALSE)* dbinom(y1,n1,p0H, log=FALSE)*dbinom(x-x1,m-m1, p0L, log=FALSE)* dbinom(y-y1,n-n1,p0H, log=FALSE);term2_p1<- dbinom(x1,m1, p1L, log=FALSE)* dbinom(y1,n1,p1H, log=FALSE)*dbinom(x-x1,m-m1, p1L, log=FALSE)* dbinom(y-y1,n-n1,p1H, log=FALSE)})?
>
>res4<-do.call(rbind,lapply(seq_len(nrow(res2)),function(i) {Pm2<-rbeta(1000,0.2+res2[i,"x"],0.8+res2[i,"m"]-res2[i,"x"]);Pn2<- rbeta(1000,0.2+res2[i,"y"],0.8+res2[i,"n"]-res2[i,"y"]); Fm2<- ecdf(Pm2); Fn2<- ecdf(Pn2); Fmm2<- Fm2(res2[i,"p1"]); Fnn2<- Fn2(res2[i,"p2"]);R2<- (Fmm2+Fnn2)/2; Fmm_f2<- min(R2, Fmm2); Fnn_f2<- max(R2, Fnn2); Qm2<- 1-Fmm_f2; Qn2<- 1-Fnn_f2;data.frame(Fmm2,Fnn2,R2,Fmm_f2,Fnn_f2,Qm2,Qn2)}))
>
>res5<-cbind(res2,res4)
>?head(res5,5)
>#? m1 n1 x1 y1 m n x y cterm1_P0L cterm1_P1L cterm1_P0H cterm1_P1H?? term2_p1
>#1? 2? 2? 0? 0 4 4 0 0???? 0.9025?????? 0.64???? 0.9025?????? 0.64 0.16777216
>#2? 2? 2? 0? 0 4 4 0 1???? 0.9025?????? 0.64???? 0.9025?????? 0.64 0.08388608
>#3? 2? 2? 0? 0 4 4 0 2???? 0.9025?????? 0.64???? 0.9025?????? 0.64 0.01048576
>#4? 2? 2? 0? 0 4 4 1 0???? 0.9025?????? 0.64???? 0.9025?????? 0.64 0.08388608
>#5? 2? 2? 0? 0 4 4 1 1???? 0.9025?????? 0.64???? 0.9025?????? 0.64 0.04194304
>#???? term2_p0?? p2?? p1 Fmm2? Fnn2???? R2 Fmm_f2 Fnn_f2?? Qm2?? Qn2
>#1 0.663420431 0.00 0.00 0.00 0.000 0.0000? 0.000? 0.000 1.000 1.000
>#2 0.069833730 0.25 0.00 0.00 0.601 0.3005? 0.000? 0.601 1.000 0.399
>#3 0.001837730 0.50 0.00 0.00 0.612 0.3060? 0.000? 0.612 1.000 0.388
>#4 0.069833730 0.00 0.25 0.59 0.000 0.2950? 0.295? 0.295 0.705 0.705
>#5 0.007350919 0.25 0.25 0.60 0.566 0.5830? 0.583? 0.583 0.417 0.417
>
>
>A.K.
>
>
>
>
>>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.
>

Zjoanna

Wed, Feb 27, 2013 9:38 AM #

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130227/55827b8c/attachment.pl>

Jeff Newmiller

Wed, Feb 27, 2013 4:06 PM #

Read The Fine Manual... the document "An Introduction to R" is supplied with R, and section 10 discusses functions. (Macros are not supported.) Depending on what you want to parameterize, a review of indexing and lists (section 2.7, 5.3, and 6) is likely to be needed as well. This is a valuable document to be familiar with, though I recommend avoiding the use of attach/detach discussed in 6.3.2.
---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
--------------------------------------------------------------------------- 
Sent from my phone. Please excuse my brevity.

Zjoanna <Zjoanna2013 at gmail.com> wrote:

________________________________
From: Joanna Zhang <[hidden

email]<http://user/SendEmail.jtp?type=node&node=4659801&i=0>>

To: arun <[hidden

email]<http://user/SendEmail.jtp?type=node&node=4659801&i=1>>

Sent: Wednesday, February 27, 2013 1:25 AM
Subject: Re: [R] cumulative sum by group and under some criteria


Thanks very much! it works!

suppose that I have a data:
m1 n1 m n    A      B    C   D
2    2    4  5  0.1    0.2  0.2 0.3
2   2     4  4  0.2   0.1   0.3  0.4
I want to identify the max value of D and extract the row , in this
example, it is the 2nd row.


On Tue, Feb 26, 2013 at 11:13 PM, arun <[hidden

email]<http://user/SendEmail.jtp?type=node&node=4659801&i=2>>

wrote:

res2<- join(res1,d3,by=c("m1","n1"),type="inner")

p0L<-0.05
p0H<-0.05
p1L<-0.20
p1H<-0.20

res2<- within(res2,{p1<- x/m; p2<- y/n;term2_p0<-dbinom(x1,m1, p0L,

log=FALSE)* dbinom(y1,n1,p0H, log=FALSE)*dbinom(x-x1,m-m1, p0L,

log=FALSE)*

dbinom(y-y1,n-n1,p0H, log=FALSE);term2_p1<- dbinom(x1,m1, p1L,

log=FALSE)*

dbinom(y1,n1,p1H, log=FALSE)*dbinom(x-x1,m-m1, p1L, log=FALSE)*
dbinom(y-y1,n-n1,p1H, log=FALSE)})

res4<-do.call(rbind,lapply(seq_len(nrow(res2)),function(i)

{Pm2<-rbeta(1000,0.2+res2[i,"x"],0.8+res2[i,"m"]-res2[i,"x"]);Pn2<-
rbeta(1000,0.2+res2[i,"y"],0.8+res2[i,"n"]-res2[i,"y"]); Fm2<-

ecdf(Pm2);

Fn2<- ecdf(Pn2); Fmm2<- Fm2(res2[i,"p1"]); Fnn2<-

Fn2(res2[i,"p2"]);R2<-

(Fmm2+Fnn2)/2; Fmm_f2<- min(R2, Fmm2); Fnn_f2<- max(R2, Fnn2); Qm2<-
1-Fmm_f2; Qn2<-

1-Fnn_f2;data.frame(Fmm2,Fnn2,R2,Fmm_f2,Fnn_f2,Qm2,Qn2)}))

res5<-cbind(res2,res4)
head(res5,5)
#  m1 n1 x1 y1 m n x y cterm1_P0L cterm1_P1L cterm1_P0H cterm1_P1H

term2_p1

#1  2  2  0  0 4 4 0 0     0.9025       0.64     0.9025       0.64

0.16777216

#2  2  2  0  0 4 4 0 1     0.9025       0.64     0.9025       0.64

0.08388608

#3  2  2  0  0 4 4 0 2     0.9025       0.64     0.9025       0.64

0.01048576

#4  2  2  0  0 4 4 1 0     0.9025       0.64     0.9025       0.64

0.08388608

#5  2  2  0  0 4 4 1 1     0.9025       0.64     0.9025       0.64

0.04194304

#     term2_p0   p2   p1 Fmm2  Fnn2     R2 Fmm_f2 Fnn_f2   Qm2   Qn2
#1 0.663420431 0.00 0.00 0.00 0.000 0.0000  0.000  0.000 1.000 1.000
#2 0.069833730 0.25 0.00 0.00 0.601 0.3005  0.000  0.601 1.000 0.399
#3 0.001837730 0.50 0.00 0.00 0.612 0.3060  0.000  0.612 1.000 0.388
#4 0.069833730 0.00 0.25 0.59 0.000 0.2950  0.295  0.295 0.705 0.705
#5 0.007350919 0.25 0.25 0.60 0.566 0.5830  0.583  0.583 0.417 0.417


A.K.

[hidden email]

<http://user/SendEmail.jtp?type=node&node=4659801&i=3>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.


______________________________________________
[hidden email]
<http://user/SendEmail.jtp?type=node&node=4659801&i=4>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.


------------------------------
 If you reply to this email, your message will be added to the
discussion
below:


http://r.789695.n4.nabble.com/cumulative-sum-by-group-and-under-some-criteria-tp4657074p4659801.html
 To unsubscribe from cumulative sum by group and under some criteria,
click

here<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4657074&code=WmpvYW5uYTIwMTNAZ21haWwuY29tfDQ2NTcwNzR8LTE3NTE1MDA0MzY=>
.

NAML<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>





--
View this message in context:
http://r.789695.n4.nabble.com/cumulative-sum-by-group-and-under-some-criteria-tp4657074p4659813.html
Sent from the R help mailing list archive at Nabble.com.
[[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.

David Winsemius

Wed, Feb 27, 2013 6:11 PM #

This is no longer a question for which the Subject is accurate. Please do not continue using inaccurate subject lines.

On Feb 27, 2013, at 9:38 AM, Zjoanna wrote:

Try:

[[1]]
  m1 n1 m n   A   B   C   D
3  2  3 4 5 0.5 0.6 0.2 0.2

[[2]]
  m1 n1 m n   A   B   C   D
4  2  4 3 4 0.2 0.8 0.5 0.3

[[3]]
  m1 n1 m n   A   B   C   D
4  2  4 3 4 0.2 0.8 0.5 0.3

[[4]]
  m1 n1 m n   A   B   C   D
2  2  2 4 4 0.2 0.1 0.3 0.4

David.


>> 
>> 
>> 
>> 
>> ________________________________
>> From: Joanna Zhang <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659801&i=0>>
>> 
>> To: arun <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659801&i=1>>
>> 
>> Sent: Wednesday, February 27, 2013 1:25 AM
>> Subject: Re: [R] cumulative sum by group and under some criteria
>> 
>> 
>> Thanks very much! it works!
>> 
>> suppose that I have a data:
>> m1 n1 m n    A      B    C   D
>> 2    2    4  5  0.1    0.2  0.2 0.3
>> 2   2     4  4  0.2   0.1   0.3  0.4
>> I want to identify the max value of D and extract the row , in this
>> example, it is the 2nd row.
>> 
>> 
>> On Tue, Feb 26, 2013 at 11:13 PM, arun <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659801&i=2>>
>> wrote:
>> 
>> res2<- join(res1,d3,by=c("m1","n1"),type="inner")
>> 
>>> 
>>> p0L<-0.05
>>> p0H<-0.05
>>> p1L<-0.20
>>> p1H<-0.20
>>> 
>>> res2<- within(res2,{p1<- x/m; p2<- y/n;term2_p0<-dbinom(x1,m1, p0L,
>> log=FALSE)* dbinom(y1,n1,p0H, log=FALSE)*dbinom(x-x1,m-m1, p0L, log=FALSE)*
>> dbinom(y-y1,n-n1,p0H, log=FALSE);term2_p1<- dbinom(x1,m1, p1L, log=FALSE)*
>> dbinom(y1,n1,p1H, log=FALSE)*dbinom(x-x1,m-m1, p1L, log=FALSE)*
>> dbinom(y-y1,n-n1,p1H, log=FALSE)})
>>> 
>>> res4<-do.call(rbind,lapply(seq_len(nrow(res2)),function(i)
>> {Pm2<-rbeta(1000,0.2+res2[i,"x"],0.8+res2[i,"m"]-res2[i,"x"]);Pn2<-
>> rbeta(1000,0.2+res2[i,"y"],0.8+res2[i,"n"]-res2[i,"y"]); Fm2<- ecdf(Pm2);
>> Fn2<- ecdf(Pn2); Fmm2<- Fm2(res2[i,"p1"]); Fnn2<- Fn2(res2[i,"p2"]);R2<-
>> (Fmm2+Fnn2)/2; Fmm_f2<- min(R2, Fmm2); Fnn_f2<- max(R2, Fnn2); Qm2<-
>> 1-Fmm_f2; Qn2<- 1-Fnn_f2;data.frame(Fmm2,Fnn2,R2,Fmm_f2,Fnn_f2,Qm2,Qn2)}))
>>> 
>>> res5<-cbind(res2,res4)
>>> head(res5,5)
>>> #  m1 n1 x1 y1 m n x y cterm1_P0L cterm1_P1L cterm1_P0H cterm1_P1H
>> term2_p1
>>> #1  2  2  0  0 4 4 0 0     0.9025       0.64     0.9025       0.64
>> 0.16777216
>>> #2  2  2  0  0 4 4 0 1     0.9025       0.64     0.9025       0.64
>> 0.08388608
>>> #3  2  2  0  0 4 4 0 2     0.9025       0.64     0.9025       0.64
>> 0.01048576
>>> #4  2  2  0  0 4 4 1 0     0.9025       0.64     0.9025       0.64
>> 0.08388608
>>> #5  2  2  0  0 4 4 1 1     0.9025       0.64     0.9025       0.64
>> 0.04194304
>>> #     term2_p0   p2   p1 Fmm2  Fnn2     R2 Fmm_f2 Fnn_f2   Qm2   Qn2
>>> #1 0.663420431 0.00 0.00 0.00 0.000 0.0000  0.000  0.000 1.000 1.000
>>> #2 0.069833730 0.25 0.00 0.00 0.601 0.3005  0.000  0.601 1.000 0.399
>>> #3 0.001837730 0.50 0.00 0.00 0.612 0.3060  0.000  0.612 1.000 0.388
>>> #4 0.069833730 0.00 0.25 0.59 0.000 0.2950  0.295  0.295 0.705 0.705
>>> #5 0.007350919 0.25 0.25 0.60 0.566 0.5830  0.583  0.583 0.417 0.417
>>> 
>>> 
>>> A.K.
>>> 
>>> 
>>> 
>>> 
>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659801&i=3>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.
>>> 
>> 
>> ______________________________________________
>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659801&i=4>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.
>> 
>> 
>> ------------------------------
>> If you reply to this email, your message will be added to the discussion
>> below:
>> 
>> http://r.789695.n4.nabble.com/cumulative-sum-by-group-and-under-some-criteria-tp4657074p4659801.html
>> To unsubscribe from cumulative sum by group and under some criteria, click
>> here<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4657074&code=WmpvYW5uYTIwMTNAZ21haWwuY29tfDQ2NTcwNzR8LTE3NTE1MDA0MzY=>
>> .
>> NAML<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>> 
> 
> 
> 
> 
> --
> View this message in context: http://r.789695.n4.nabble.com/cumulative-sum-by-group-and-under-some-criteria-tp4657074p4659813.html
> Sent from the R help mailing list archive at Nabble.com.
> 	[[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.

David Winsemius
Alameda, CA, USA