Sorry here's working code, I wrote it from my iPhone before!
x = runif(100, 0, 10)
plot(x)
abline(lsfit(1:100, x))
abline(lsfit(1:100, x + sd(x)))
fit = lsfit(1:100, x + sd(x))
y = which(residuals(fit) > 0)
points(y, x[y], col="red")
Notice the points that are red are all the ones above the first lsfit. I'd
like only those above the second lsfit that was drawn on the plot.
David Winsemius wrote:
On Oct 25, 2009, at 1:20 PM, Tanya Cashorali wrote:
I'm trying to extract the points above and below a particular lsfit. I can
only get the residuals from the original fit though.
x = runif(100, 0, 10)
plot(x)
abline(lsfit(1:100, test))
Error in as.matrix(y) : object 'test' not found
I ran the first three lines and was not in the least bit surprised that I
got an error. I was expecting that lsfit needed to see some object (perhaps
"x") rather than an integer series, and I didn't see "test" anywhere
defined.
abline(lsfit(1:100, test + sd(test))) #I want the points above THIS line.
Did you really have running code that produced a line?
Is there a way to use the coefficients from the fit to do this?
Thanks for any help.
If you give lsfit a proper input, and assign the result to a fit object, say
"xfit", then coef(xfit) will give you the coefficients, or if you wanted to
apply some sort of test on the residuals, you could get them with either
residuals(xfit) or xfit$residuals.
On Mon, Oct 26, 2009 at 7:48 AM, David Winsemius <dwinsem...@comcast.net>wrote:
>
> On Oct 25, 2009, at 1:20 PM, Tanya Cashorali wrote:
>
> I'm trying to extract the points above and below a particular lsfit. I can
>> only get the residuals from the original fit though.
>>
>> x = runif(100, 0, 10)
>> plot(x)
>> abline(lsfit(1:100, test))
>>
>
> Error in as.matrix(y) : object 'test' not found
>
> I ran the first three lines and was not in the least bit surprised that I
> got an error. I was expecting that lsfit needed to see some object (perhaps
> "x") rather than an integer series, and I didn't see "test" anywhere
> defined.
>
>
> abline(lsfit(1:100, test + sd(test))) #I want the points above THIS line.
>>
>
> Did you really have running code that produced a line?
>
>
>
>> Is there a way to use the coefficients from the fit to do this?
>> Thanks for any help.
>>
>
> If you give lsfit a proper input, and assign the result to a fit object,
> say "xfit", then coef(xfit) will give you the coefficients, or if you wanted
> to apply some sort of test on the residuals, you could get them with either
> residuals(xfit) or xfit$residuals.
>
> --
>
> David Winsemius, MD
> Heritage Laboratories
> West Hartford, CT
>
>
[[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.
