

Begin forwarded message:

> From: stephen sefick <ssefick@gmail.com>
> Date: January 28, 2011 9:11:07 AM CST
> To: "Claudia Paladini" <claudiapaladini@web.de>
> Subject: Re: [R] mvoutlier
>
> Try this.
>
> #your code
> library(mvoutlier)
> library(robustbase)
> data(bsstop)
> x=bsstop[1:100,5:14]
> x.out=sign1(x,makeplot=TRUE)
>
> #my code
> y <- data.frame(distance=x.out$x.dist, constant=x.out$const)
> y <- data.frame(y, row=rownames(y))
>
> qplot(row, distance, data=y, label=row, geom="text") 
> +geom_hline(yintercept=y[1,"constant"])
>
>
> On Jan 28, 2011, at 3:29 AM, Claudia Paladini wrote:
>
>> library(mvoutlier)
>> library(robustbase)
>> data(bsstop)
>> x=bsstop[1:100,5:14]
>> x.out=sign1(x,makeplot=TRUE)
>


	[[alternative HTML version deleted]]

