[R] superimpose graphs
jim holtman
jholtman at gmail.com
Sun Apr 17 03:47:37 CEST 2011
This will add the regression line to it:
> x <- read.table(textConnection(" id invest payout
+ 1 10 -1
+ 1 33 33
+ 1 20 -5
+ 2 200 33
+ 2 33 -20
+ 3 5 -5
+ 3 5 -5"), header = TRUE)
> closeAllConnections()
> # normalize 'invest' by the mean
> x$investNorm <- ave(x$invest, x$id, FUN = function(a) a / mean(a))
> x
id invest payout investNorm
1 1 10 -1 0.4761905
2 1 33 33 1.5714286
3 1 20 -5 0.9523810
4 2 200 33 1.7167382
5 2 33 -20 0.2832618
6 3 5 -5 1.0000000
7 3 5 -5 1.0000000
> # plot the points
> plot(x$payout, x$investNorm)
> abline(lm(x$investNorm ~ x$payout))
>
On Sat, Apr 16, 2011 at 1:17 PM, 苏江东Su Jiangdong <sujiangdong at gmail.com> wrote:
> Hi there,
>
> I have a data frame DF of over 600 people's short term trade data in time
> order. Below is the simplified structure of the data.
>
> id invest payout
> [1] 1 10 -1
> [2] 1 33 33
> [3] 1 20 -5
> [4] 2 200 33
> [5] 2 33 -20
> [6] 3 5 -5
> [7] 3 5 -5
>
> id is each person's id. Each person have invested many times in the sampling
> period, in temporal order.
>
> What I want to check is the correlation between invest and payout.
>
> 1. How do I run the regression for each person, with the "invest" being
> devided by the mean or medium of the person's "invest"?
> 2. How do I plot a graph with y axis being invest/mean(invest) and x axis
> being payout, all 600 people's dots superimposed on one graph?
>
> I tried to use
>
> for (i in 1:(dim (DF)[1]-1))
> {
> if (DF[i,1]=DF[i+1,1]) id.lm <- lm(invest ~ payput, data=DF)
> }
>
> But I don't know how to superimpose graphs onto each other.
>
> Thanks a lot!
>
> Su
>
> [[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.
>
--
Jim Holtman
Data Munger Guru
What is the problem that you are trying to solve?
More information about the R-help
mailing list