Message-ID: <125442c7-2f5c-4591-a6f1-03a768c5ea75@email.android.com>
Date: 2013-02-28T00:06:49Z
From: Jeff Newmiller
Subject: cumulative sum by group and under some criteria
In-Reply-To: <CACexUksKSOPB5V9e=n1jeqhRN4fQA6tEKaNvA15kxNCckCseoQ@mail.gmail.com>
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:
>Hi,
>I need to change one or several parameters and run the same code
>several
>times, how to create a macro with parameters in R?
>
>On Wed, Feb 27, 2013 at 7:45 AM, arun kirshna [via R] <
>ml-node+s789695n4659801h42 at n4.nabble.com> wrote:
>
>> 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 <[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.