Displaying 20 results from an estimated 900 matches similar to: "outlier identify in qqplot"
2012 May 11
1
identify() doesn't return "true" numbers
Dear R community.
I am using the identify() function to identify outliers in my dataset.
This is the code I am using:
####################################################################
# Function to allow identifying points in the QQ plot (by mouseclicking)
qqInteractive <- function(..., IDENTIFY = TRUE)
{
qqplot(...) -> X
abline(a=0,b=1)
if(IDENTIFY) return(identify(X))
2010 Mar 27
5
producing a QQ plot.
Hello everyone I'm a beginner in Stats and R, I'm using R 2.10.1. I need to
create a multivariate qq plot, there is 8 variable group with each has 55
number of input. An example of what I did so far, just to get my point out:
> data=read.csv(file.choose(),header=T)
> data
country village group av_expen P2ary_ed no_fisher
1 Cook Islands Aitutaki D
2012 Feb 10
1
making multiple lines using qqplot
Hi Everyone,
I want to make 3 lines on the same graph (not as subplots, all in the same
window, one on top of each other) and I want them to be quantile-quantile
plots (qqplot). Essentially, I am looking for the equivalent of Matlab's
"hold on" command to use with qqplot. I know I can use 'points' or 'lines',
but these do not give me a qqplot (only appear to work
2006 Feb 06
2
qqplot
Hello,
I would like to use qqplot() to compare two
experimental distributions. But I do not understand
how qqplot() compute quantiles. In fact, quantile() do
not return the same results.
Thank you for your help.
Vincent.
2011 Mar 25
1
multiple plots with QQplot of PerformanceAnalytics
Hi All,
I am trying to plot 4 graphs on to 1 page using layout(...), or par(mfcol =
c(...)); with the function QQplot from the package PerformanceAnalytics.
The problem is that, no matter what order I use, it only plots 3 graphs on to 1
page and the last QQplot is shunted to the next page.
Also, this only happens to the QQplot, i.e. there is no problem with 4
Histograms.
set.seed(1033)
data
2001 Jan 05
1
pairs(NxK_Matrix,panel=qqplot) (PR#803)
Full_Name: Matthias von Davier
Version: 1.2.0
OS: linux
Submission from: (NULL) (144.81.31.148)
pairs(NxK_Matrix,panel=qqplot)
produces a message
Error in pairs.default(NxK_Matrix,panel=qqplot) :
The panel function made a new plot
best regards and a happy new year
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read
2009 Dec 23
2
how to create normal qqplot with the 95% confidence interval
hi everyone!
season's greetings!
is there any way that i can create a normal qqplot showing, aside from the
qqline, the 95% confidence limits? thank you very much..
happy holidays!
--
View this message in context: http://n4.nabble.com/how-to-create-normal-qqplot-with-the-95-confidence-interval-tp977727p977727.html
Sent from the R help mailing list archive at Nabble.com.
2010 Jun 24
4
Simple qqplot question
I am a beginner in R, so please don't step on me if this is too
simple. I have two data sets datax and datay for which I created a
qqplot
qqplot(datax,datay)
but now I want a line that indicates the perfect match so that I can
see how much the plot diverts from the ideal. This ideal however is
not normal, so I think qqnorm and qqline cannot be applied.
Perhaps you can help?
Ralf
1999 Aug 31
1
Suggestion for qqplot() improvement
>>>>> On Tue, 31 Aug 1999 14:57, Werner Stahel <stahel@stat.math.ethz.ch> said:
WSt> Here is a suggestion. It seems that qqplots, comparing a sample
WSt> to a distribution other than the normal, are not explicitly
WSt> available in S or R. I found
(in S-plus / Trellis it is, see below)
WSt> qqplot(y, rt(300, df = 5))
WSt> as an
2009 Nov 02
7
qqplot
Hi,
We could use qqplot to see how two distributions are different from each other. To show better how they are different (departs from the straight line), how is it possible to plot the straight line that goes through them? I am looking for some thing like qqline for qqnorm. I thought of abline but how to determine the slope and intercept?
Best wishes,
Carol
2004 Jan 21
1
outlier identification: is there a redundancy-invariant substitution for mahalanobis distances?
Dear R-experts,
Searching the help archives I found a recommendation to do multivariate
outlier identification by mahalanobis distances based on a robustly estimated
covariance matrix and compare the resulting distances to a chi^2-distribution
with p (number of your variables) degrees of freedom. I understand that
compared to euclidean distances this has the advantage of being scale-invariant.
2012 Nov 29
1
QQplot
Hi!
<http://r.789695.n4.nabble.com/file/n4651293/qq.png>
We are stuck with a problem considering the qqplot of a dataset.
We are trying to discover what kind of distribution this is. We already
tried to normal, exponential or the logaritmical distribution but none of
those are able to solve our problem. Is there someone able to tell us what
kind deformation we should try?
(I'm
2006 Mar 15
1
(newbie) Weighted qqplot?
Folks,
Normally, in a data frame, one observation counts as one observation
of the distribution. Thus one can easily produce a CDF and (in Splus
atleast) use cdf.compare to compare the CDF (BTW: what is the R
equivalent of the SPlus cdf.compare() function, if any?)
However, if each point should not count equally, how can I weight the
points before comparing the distributions? I was thinking of
2005 Aug 04
1
some thoughts on outlier detection, need help!
Dear listers:
I have an idea to do the outlier detection and I need to use R to
implement it first. Here I hope I can get some input from all the
guru's here.
I select distance-based approach---
step 1:
calculate the distance of any two rows for a dataframe. considering
the scaling among different variables, I choose mahalanobis, using
variance as scaler.
step 2:
Let k be the number of
2012 Sep 14
2
when to use "I", "as is" caret
Dear community,
I've check it while working, but just to reassure myself. Let's say we have
2 models:
model1 <- lm(vdep ~ log(v1) + v2 + v3 + I(v4^2) , data = mydata)
model2 <- lm(vdep ~ log(v1) + v2 + v3 + v4^2, data = mydata)
So in model1 you really square v4; and in model2, v4*^2 *doesn't do
anything, does it? Model2 could be rewritten:
model2b <- lm(vdep ~
2011 Sep 26
2
Mahalanobis Distance
Hello R helpers,
I'm trying to use Mahalanobis distance to calculate distance of two time
series, to make some comparations with euclidean distance, DTW, etc, but I'm
having some dificults.
I have, for example, two objects:
s.1 <- c( 5.6324702, 1.3994353, -3.2572327, -3.8311846, -1.2248719,
0.9894694, -2.2835332, -5.1969285, -5.2823988, -3.1499400, -1.7307950,
2.8221209,
2010 Jan 30
2
Questions on Mahalanobis Distance
Hello,
I am a new R user and trying to learn how to implement the mahalanobis
function to measure the distance between to 2 population centroids. I
have used STATISTICA to calculate these differences, but was hoping to learn
to do the analysis in R. I have implemented the code as below, but my
results are very different from that of STATISTICA, and I believe I may not
have interpreted the help
2004 Mar 26
1
Mahalanobis
Dear all
Why isn'it possible to calculate Mahalanobis distances with R for a matrix
with 1 row (observations) more than the number of columns (variables)?
> mydata <- matrix(runif(12,-5,5), 4, 3)
> mahalanobis(x=mydata, center=apply(mydata,2,mean), cov=var(mydata))
[1] 2.25 2.25 2.25 2.25
> mydata <- matrix(runif(420,-5,5), 21, 20)
> mahalanobis(x=mydata,
2011 Mar 22
1
Using the mahalanobis( ) function
Hello all,
I am a 2 month newbie to R and am stumped. I have a data set that I've run multivariate stats on using the manova function (I included the data set). Now it comes time for a table of effect sizes with significance. The univariate tests are easy. Where I run into trouble filling in the table of effect sizes is the Mahalanobis D as an effect size. I've included the table so
2009 Jul 20
2
mahalanobis distance
http://www.nabble.com/file/p24569511/mahalanobis.txt mahalanobis.txt
http://www.nabble.com/file/p24569511/concentrations.txt concentrations.txt
Dear Forum members,
I have a problem calculating mahalanobis distances. My data file
mahalanobis.txt and categories file concentrations.txt are attached. I do
the following steps:
x <- as.matrix(read.table("mahalanobis.txt", header=TRUE))