Displaying 20 results from an estimated 6000 matches similar to: "List seems to drop empty levels of factors when containin g them"
2005 Apr 19
2
Odd diagnostic plots in mixed-effects models
Dear R community,
In the excellent nlme package the default diagnostic plot graphs the innermost residuals against innermost fitted values. I recently fit a mixed-effects model in which there was a very clear positive linear trend in this plot.
I inferred that this trend occurred because my fixed effect was a two-level factor, and my random effect was a 12-level factor. The negative residuals
2003 Mar 08
2
Looking for non-central F quantile
Greetings all,
I'm trying to figure out how to calculate the inverse CDF (i.e. a
quantile) for a non-central F distribution. I could put together a quick
numerical solver routine using the CDF, but I wonder if there's a function
that I've missed that would be more efficient?
Thank-you,
Andrew
Andrew Robinson Ph: 208 885 7115
Department of Forest Resources Fa: 208 885
2002 Jul 18
3
Oddity with names
Hi all,
I'm using R 1.5.1 on Windows 2000. The following snippet of code doesn't
seem to do anything - no error is reported, and there is no name change.
names(myFrame[,c(1:3)]) <- c("name1", "name2", "name3")
This code however works nicely:
names(myFrame)[c(1:3)] <- c("name1", "name2", "name3")
Can anyone suggest why
2004 Nov 10
1
List seems to drop empty levels of factors when containing them
Greetings R community,
I am curious about the following behaviour: if I define a factor, and
then store a subset of it in a list, the stored version seems to drop
levels that were not included in the subset. E.g. ..
> mask <- c(T, F)
> grp.1 <- factor(c(1,2))
> list(grp.1)
[[1]]
[1] 1 2
Levels: 1 2
> list(grp.1[mask])
[[1]]
[1] 1
Levels: 1
It is as though the list were
2005 Jun 15
2
Plotting second axes outside xyplot
Hi all,
I'm trying to find a way to get xyplot to produce a second set of axes outside the right hand side of the graph. This is my progress so far:
EE <- equal.count(ethanol$E, number=9, overlap=1/4)
xyplot(NOx ~ C | EE, data = ethanol,
prepanel = function(x, y) prepanel.loess(x, y, span = 1),
xlab = "Compression Ratio", ylab = "NOx (micrograms/J)",
2004 Mar 23
2
Coefficients and standard errors in lme
Hello,
I have been searching for ways to obtain these for combinations of fixed
factors and levels other than the 'baseline' group (contrasts coded all
0's) from a mixed-effects model in lme. I've modelled the continuous
variable y as a function of a continuous covariate x, and fixed factors
A, B, and C. The fixed factors have two levels each and I'd like to know
whether
2004 Sep 21
2
Bootstrap ICC estimate with nested data
I would appreciate some thoughts on using the bootstrap functions in the
library "bootstrap" to estimate confidence intervals of ICC values
calculated in lme.
In lme, the ICC is calculated as tau/(tau+sigma-squared). So, for instance
the ICC in the following example is 0.116:
> tmod<-lme(CINISMO~1,random=~1|IDGRUP,data=TDAT)
> VarCorr(tmod)
IDGRUP = pdLogChol(1)
2004 Mar 10
3
converting lists got by tapply to dataframes
I have two lists:
xa <- list( X=c(1,2,3), Y=c(4,5,6), Z=c(7,8,9) )
xb <- with( barley, tapply( X=seq(1:nrow(barley)), INDEX=site
, FUN=function(z)yield[z]))
I can convert xa to a dataframe easily with:
as.data.frame(xa)
But if i try the same with xb I get:
as.data.frame(xb)
Error in as.data.frame.default(xb) :
can't coerce array into a data.frame
What
2004 Mar 26
8
stop() vs. error() ?
Why does stop("we are done") print
"Error in eval.with.vis(expr, envir, enclos) :" ?
It would seem to me that a plain stop() is not an error, and that it
would make more sense to have an error() function that is different from
a stop(). Is there a rationale here that I am missing?
sincerely, /iaw
2002 Jul 07
2
R on Sharp Zaurus 5500?
Hello,
has anyone tried to get R running on a Sharp Zaurus 5500?
Andrew
Andrew Robinson Phone: 208-885-7115
Department of Forest Resources Fax: 208-885-6226
University of Idaho E: andrewr at uidaho.edu
Po Box 441133 WWW: http://www.uidaho.edu/~andrewr
Moscow, ID 83843 and: http://www.biometrics.uidaho.edu/
No statement above
2003 Sep 01
1
par(new=T) works differently in pdf vs postscript if applied before a plot statement. (PR#4037)
If I place par(new=T) before I create a plot in a script that is sent to a pdf
device, the pdf is unopenable and reports itself as having no pages. The
postscript device seems to ignore the par instruction.
I guess one of these is a bug, but I don't know which one!
--please do not edit the information below--
Version:
platform = i386-pc-linux-gnu
arch = i386
os = linux-gnu
system =
2003 Apr 13
2
Peculiarity in non-central qchisq for ncp > 294.92 ...
Hello all,
Here's my query:
Running R 1.6.2 on FreeBSD 5.0, and on WinXP, and I find that the
following hangs the process:
dchisq(alpha=0.01, df=1, ncp=295)
but it does work for ncp < 294.92.
Is this general?
Best wishes to all,
Andrew
Andrew Robinson Ph: 208 885 7115
Department of Forest Resources Fa: 208 885 6226
University of Idaho E : andrewr at uidaho.edu
PO
2005 Jan 02
2
How to quieten axis() for Sweave: avoid echoing NULL?
Dear R-community,
I'm using Sweave to produce reports. The reports require the "axis"
command. When I run axis the program returns NULL as well as creating
the axis.
> plot(1:4, rnorm(4), axes=FALSE)
> axis(1, 1:4, LETTERS[1:4])
NULL
>
So, my Sweave tex files have
\begin{Schunk}
\begin{Soutput}
NULL
\end{Soutput}
\end{Schunk}
in front of each graphic that requires
2004 Jan 03
1
R doesn't run on g4 iBook with Panther
I have a 12" G4 iBook with the latest os (Panther) with the recent
updates. When I start R by double clicking, it expands for a fraction of
a second and then, nothing happens.
There is mention in the readme about using something called "i-install"
to uninstall a mistakenly installed library, but I do not know how to do
this.
Any help would be appreciated.
Thanks,
Martin
2004 Oct 04
1
Could anyone tell me how to extract pvalue from "lm" fitting?
Dear R people,
I have a naive question: after fitting "lm" to a data, I can't extract
the pvalue corresponding to a specific covariate in a direct way.
Could anyone give me a hint?
Thank you very much.
Frank
2005 Jan 04
1
Histogram
I have data on a single variable LOGT. It has about
300,000 observations. I am trying to make a Histogram
out of this data set. Following is my effort. Could
anyone help me to solve this error.
> hist(x)
Error in hist.default(x) : `x' must be numeric
> class(x)
[1] "data.frame"
> is.object(x)
[1] TRUE
> is.vector(x)
[1] FALSE
> is.numeric(x)
[1] FALSE
>
2004 Oct 03
1
help working with persp plot function
Hi,
I don't understand why this is not working. Help is appreciated.
I need to plot the following as a surface, but persp returns an error.
tpcp_xy[1:10,]
X_COORD Y_COORD TPCP
1 465459.7 175924.1 0.85
2 466145.8 324017.3 2.30
3 467720.2 372143.1 1.56
4 470293.2 348064.8 2.87
5 476566.9 205501.8 0.94
6 477774.9 142561.0 1.31
7 479207.0 162919.6 3.04
8 480890.8 290641.3 2.20
9
2004 Mar 12
2
Problems startin RAqua for Mac
I am a user of a i book G4 with OS X 10.3 and I installed RAqua 1.8.1
without any problems. When I want to start RStart the "R" Symbol appers
for one second and then disappers again and nothing else happens. I've
tried to reinstall the program and I also reinstalled my OS X wich did
not help. I really need this program for my studies. If someone knows
whats wrong please eMail
2004 Dec 13
2
Farmating printed Numbers and Text
Hello,
is there any function to print (using "cat" or "print" etc.) a table of
values row per row in a exact order, even if the value has 3 digit and in
the next row eg. 5 digits? This means I want to print a number with 3 digits
with 2 spaces in front and a number with 4 digits only with one space in
front.
regards Andreas
2005 Jan 01
1
lme: Variances
Hi R users!
I will try to state my question again. I have longitudinal data and fitted the
following model with lme:
Y = X*beta + U + W(t) + Z
where
U ~ N(0, nu*I) I is the identity matrix, so this is the random intercept
W(t)~ N(0, sigma*H) and H is a matrix which incorporates a Gaussian serial
correlation (covariance) in the offdiagonal