Alas, speaking before I check the optional arguments....c'est la vie.

MW

On Sun, Sep 18, 2011 at 3:00 PM, Jorge I Velez <jorgeivanvelez@gmail.com>wrote:

> Why wouldn't?  It does for me ;-)
>
> curve(sin(x), 0, 10, col = 4, lwd = 2, las = 1)
> curve(cos(x), 0, 10, add = TRUE)
>
> Best,
> Jorge
>
>
> On Sun, Sep 18, 2011 at 2:07 PM, R. Michael Weylandt <> wrote:
>
>>  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 <>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
>> > pburns@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.
>>
>
>

	[[alternative HTML version deleted]]

