Hi Steve,
Just test whether y is greater than the predicted y (i.e., your line).
## function using the model coefficients*
f <- function(x) {82.9996 + (.5589 * x)}
## Find group membership
group <- ifelse(y > foo(x), "A", "B")
*Note that depending how accurate this needs to be, you will probably
want to use the model itself rather than just reading from the
printout like I did. If you need to do that, take a look at ?predict
For future reference, it would be easier for readers if you provided
your data via something like: dput(x) that can be copied directly into
the R console. Also, if you are generating random data (rnorm()), you
can use set.seed() so that we can replicate exactly what you get.
HTH,
Josh
On Fri, Mar 4, 2011 at 1:39 PM, Steve Hong <empti...@gmail.com> wrote:
> Hi R-list,
>
> I have a data set with plot locations and observations and want to label
> them based on locations. For example, I have GPS information (x and y) as
> follows:
[snip]
>> (fm1 <- lm(ysim~xsim))
> Call:
> lm(formula = ysim ~ xsim)
> Coefficients:
> (Intercept) xsim
> 82.9996 0.5589
>
> I overlapped fitted line on the plot.
>
>> abline(fm1)
> My question is:
> As you can see in the plot, how can I label (or re-group) those in upper
> diagonal as (say) 'A' and the others in lower diagonal as 'B'?
>
> Thanks a lot in advance!!!
>
> Steve
>
> [[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.
>
--
Joshua Wiley
Ph.D. Student, Health Psychology
University of California, Los Angeles
http://www.joshuawiley.com/
______________________________________________
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.
