Furthermore, if you want to add lines to a plot, I'd advise against the curve() function, rather use lines() to add lines (or points() for, well, points) on a plot.
E.g., x <- seq(0,10,by=0.025) y <- sin(x) y2 <- cos(x) plot(x,y,type="l",col="red4",lwd=3) lines(x,y2) curve() wouldn't work for this Michael On Sun, Sep 18, 2011 at 1:23 PM, Patrick Burns <pbu...@pburns.seanet.com>wrote: > Sarah has told you how to get from where > you are to where you want to be. > > However, it is easier to start somewhere > else. The more "R" way to do what you are > doing is to use a list to put the results > of your regressions into. It is then very > easy to use that list, and it is easier to > keep track of. > > The two documents mentioned in my signature > may help you get up and running in R. Since > you are coming from Stata, you should also > look at http://r4stats.com > > > On 18/09/2011 17:05, Francesco Sarracino wrote: > >> Dear all, >> >> I am a stata user and I am moving my first steps in R. >> I am dealing with a silly issue concerning looping over variables. I read >> previous posts on similar topics in the R help archive, but I did not find >> a >> solution. >> Here is my case: >> >> I run a simple bivariate linear regression saving the output in objects >> called: pd.memb.0, pd.memb.1, pd.memb.2, pd.memb.3 >> >> pd.memb.0<- lm(E.fs.memb ~ M.fs.memb, data, subset = (tag2 == 1)) >> pd.memb.1<- lm(E.fs.memb ~ M.fs.memb, data, subset = (M.fs.memb<= >> quantile(M.fs.memb, .25)& tag2 == 1)) >> pd.memb.2<- lm(E.fs.memb ~ M.fs.memb, data, subset = (M.fs.memb> >> quantile(M.fs.memb, .25)& M.fs.memb<= quantile(M.fs.memb, .75)& tag2 == >> 1)) >> pd.memb.3<- lm(E.fs.memb ~ M.fs.memb, data, subset = (M.fs.memb> >> quantile(M.fs.memb, .75)& tag2 == 1)) >> >> Subsequently, I wish to plot my data superimposing a different line for >> each >> regression model. Here I am trying to apply my loop. Looking at the info I >> found around, the following code sound to me as correct, but I am wrong. >> >> plot(M.fs.memb , E.fs.memb) >> for(i = 1:3){ >> curve(coef(pd.memb.i)[1] + coef(pd.memb.i)[2]*x, add = T, cex = 100, col >> = >> "red") >> } >> >> The plot is plotted, but the curves are not printed and I get the >> following >> error message: >> >> >> Error in coef(pd.memb.i) : >> error in evaluating the argument 'object' in selecting a method for >> function 'coef'> }Error: unexpected '}' in "}" >> >> >> What am I doing wrong? >> Thanks a lot for your help, >> f. >> >> [[alternative HTML version deleted]] >> >> ______________________________**________________ >> R-help@r-project.org mailing list >> https://stat.ethz.ch/mailman/**listinfo/r-help<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. >> >> > -- > Patrick Burns > pbu...@pburns.seanet.com > twitter: @portfolioprobe > http://www.portfolioprobe.com/**blog <http://www.portfolioprobe.com/blog> > http://www.burns-stat.com > (home of 'Some hints for the R beginner' > and 'The R Inferno') > > > ______________________________**________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/**listinfo/r-help<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. > [[alternative HTML version deleted]] ______________________________________________ R-help@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.
