Displaying 20 results from an estimated 5000 matches similar to: "qqnorm.lme & pairs.lme"
2000 Jul 04
1
nlme errors ?
Dear friends. Below is ouput directly from the help on qqnorm.lme.html
It do not seem to work as expected - on win98, R 1.1, although the
regression is undertaken, so qqnorm
misunderstands, or what ?
library(nlme)
data(Orthodont)
fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject)
# normal plot of standardized residuals by gender
qqnorm(fm1, ~ resid(., type = "p") | Sex,
2010 Apr 15
1
classes and functions for qqnorm and stem
Referring to "Using R for Data Analysis and Graphics" by J H Maindonald,
and available from the R site, I found the example on p.30 non-working:
> stem(qqnorm(possum$hdlngth))
Error in stem(qqnorm(possum$hdlngth)) : 'x' must be numeric
Since qqnorm(possum$hdlngth) plots, and
> class(possum$hdlngth)
[1] "numeric"
, the problem must be here:
>
2005 Apr 25
2
residuals in lmer
Does anyone know how to extract residuals in lmer?
Here's the error I get:
>
crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep:
pesticide)+(1|rep:pesticide:irrigation), crop.data)
> qqnorm(crop.lme)
Error in qqnorm.default(crop.lme) : y is empty or has only NAs
> resid(crop.lme)
NULL
Thanks!
--Jake
2007 Jun 29
2
why this doesn't work for qqnorm
I want to qqnorm every column in a table. When I try the first column
using
qqnorm(table$column1), it worked.
But when I use
qqnorm(table[1]), it tells me "Error in stripchart(x1, ...) : invalid
plotting method".
What happen? How can I make a function that qqnorms every column?
thanks a lot.
-jiong
The email message (and any attachments) is for the sole use of the intended
2001 Nov 27
2
overlaying qqnorm plots...
I know this topic has had plenty of discussion in the last couple of days,
but....
I've been trying to compare the effects of different fitted methods for
systems of equations (OLS, SUR, 2SLS, 3SLS ) and would like to compare the
residual plots (easy) and the qqnorm/qqplot of the fits for the different
fitted methdos. For example,
qqnorm( residuals( lm( q ~ p + f + a ) ) )
par( new = TRUE )
2000 Sep 21
2
qqnorm(), is it "backwards"?
Hello R friends,
I'm wondering why I get funny qqnorm() results. It seems that they should
all be reflected in the normal qqline().
For instance: if I qqnorm() bimodal or uniform data I get a sigmoidal in
which the qqnorm() points lie above the qqline() at -ve theoretical
quantiles, and the qqnorm() points lie below the qqline() at +ve
theoretical quantiles. Yet I expect such platykurtic
1999 Mar 16
1
qqnorm in R-0.63.3
Dear List,
invoking qqnorm-plots in Version 63.3 produces funny things:
using the option `type="s"ยด on qqnorm should give a nice *line* of
observed quantiles. Now, the line is walking along in order to the
points index instead from lowest to highest, wich makes funny slopes.
try x <- table(rnorm(1000) # or similar and
qqnorm(x,type="s") # in 0.63.2 and 63.3
Well, the
2011 Jul 25
1
scripting qqplot and qqnorm
Hi all,
I am an R newbie, and I have a question about scripting. I have the
following lines which I want to put int\
o a script which I can call from the shell of a Mac/Linux machine :
myrns <- read.csv(file="/Users/vihan/test.csv",sep="",header=FALSE)
qqnorm(myrns)
qqline(myrns)
This works fine on an interactive R session.
However, as I understand the basic framework of
2009 Sep 08
1
Plotting two qqnorm plots:
Hi all,
Does anyone know how to plot overlapping qqnorm plots on the same window?
Suppose I have data in the vector x and y:
qqnorm(x)
lines(qqnorm(y))
I though these two lines will do the job... However, lines doesn't seem to
work. Anyways, thanks in advance!
--
View this message in context: http://www.nabble.com/Plotting-two-qqnorm-plots%3A-tp25352893p25352893.html
Sent from the R help
2009 Apr 26
3
Flipping axes of qqnorm
Hi all,
I have just started using R to produce qqnorm plots. I am trying to
switch the x and y axes so that the theoretical values are plotted on the y
axis and my data on the x axis. Can anyone help me with this?
Thanks
--
View this message in context: http://www.nabble.com/Flipping-axes-of-qqnorm-tp23248007p23248007.html
Sent from the R help mailing list archive at Nabble.com.
2005 Apr 03
2
how to draw a 45 degree line on qqnorm() plot?
# I can not draw a 45 degree line on a qqnorm() plot,
jj <- sample(c(1:100), 10)
qqnorm(jj)
abline() don't work.
Thank you.
2008 Oct 22
1
R 2.8.0 qqnorm produces error with object of class zoo?
Dear list-reader,
by running the following script:
library(zoo)
sessionInfo()
search()
packageDescription("zoo")
data(EuStockMarkets)
dax <- as.zoo(EuStockMarkets[1:10, "DAX"])
daxr <- diff(log(dax))
identical(as.vector(qnorm(daxr)), qnorm(coredata(daxr)))
qqnorm(coredata(daxr))
qqnorm(daxr)
qqnorm() produces an error:
> qqnorm(daxr)
Fehler in if (xi == xj) 0L
2007 Feb 08
1
the plotting position of theoretical quantile for qqnorm
Hello,
I have a doubt about the plotting position of the theoretical quantile for
the qqnorm
command in R.
Let F be the theoretical distribution of Y, we observed a sample of size n,
y1,y2, ...,
yn. We then sort it and comspare these empirical quantiles to the expected
ones
from F. For the plotting poition, there are several options:
1. i/(n+1)
2. (i-.375)/(n+.25)
3. (i- .3175)/ (n + .365)
etc.
2005 Apr 22
1
an interesting qqnorm question
Hi, r-gurus:
I happened to have a question in my work:
I have a dataset, which has only one dimention, like
0.99037297527605
0.991179836732708
0.995635340631367
0.997186769599305
0.991632565640424
0.984047197106486
0.99225943762649
1.00555642128421
0.993725402926564
....
the data is saved in a file called f392.txt.
I used the following codes to play around :)
2003 May 12
1
plot.ranef.lme (PR#2986)
library(nlme)
data(Phenobarb)
na.include <- function(x)x
phe1 <- nlme(conc~phenoModel(Subject, time, dose, lCl, lV),
data = Phenobarb,
fixed = lCl+lV~1,
random= pdDiag(lCl+lV~1),
start = c(-5,0),
na.action = na.include,
naPattern = ~!is.na(conc))
phe.ranef <- ranef(phe1,augFrame=TRUE)
plot(phe.ranef, form=lCl~Wt+ApgarInd)
[Error in max(length(x0),
2006 Apr 13
2
Plotting positions in qqnorm?
Do you know of a reference that discusses alternative choices for
plotting positions for a normal probability plot? The documentation for
qqnorm says it calls ppoints, which returns qnorm((1:m-a)/(m+1-2*a))
with "a" = ifelse(n<=10, 3/8, 1/2)? The help pages for qqnorm and
ppoints just refer to Becker, Chambers and Wilks (1988) The New S
Language (Wadsworth & Brooks/Cole),
2006 Jul 24
3
standardized random effects with ranef.lme()
Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can
obtain random effects for intercept and slope of a certain level (say:
1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm
mistaken here, but the results are identical.
However, if I try to get the standardized random effects adding the
paramter "standard=T" to the
2006 Dec 08
1
How to plot two variables in the same qqnorm-plot?
Dear all,
I have two variables called c2 and c3 and want to plot these variables in the same qqnorm-plot with two different symbols or colors to distinguish them so I can easily compare the variables aginst each other. How can I do in R? I only manage to do two separated qqnorn-plots.
Thanks for your help,
All the best,
Jenny
c2=
-0.1545775144 -0.0601161235
2012 May 25
1
difference between qnorm and qqnorm
dear all,
it will just take you a minute to tell me the difference
between qnorm and qqnorm. are they same or is there any
difference between them??
regards
2011 Dec 21
4
qqnorm & huge datasets
Hi,
When qqnorm on a vector of length 10M+ I get a huge pdf file which
cannot be loaded by acroread or evince.
Any suggestions? (apart from sampling the data).
Thanks.
--
Sam Steingold (http://sds.podval.org/) on Ubuntu 11.10 (oneiric) X 11.0.11004000
http://mideasttruth.com http://honestreporting.com http://camera.org
http://openvotingconsortium.org http://pmw.org.il