Displaying 20 results from an estimated 30000 matches similar to: "reading file names for R batch process"
2001 Nov 22
4
changing the magnification of axis annotation
Hi all,
how do I change the magnification of xaxis annotation? I need a smaller
axis text size in some of my plots. I tried cex.axix=0.5 in plot(), but
this doesn't work.
Thanks, Sven
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or
2002 May 29
4
Why is.integer() doesn't work with single values?
Hi all,
I don't understand the behavior of is.integer():
> x <- integer()
> is.integer(x)
[1] TRUE
> x <- 10
> is.integer(x)
[1] FALSE
> x <- 1:10
> is.integer(x)
[1] TRUE
Why is.interger() returns FALSE if x has only one element? And how can
someone check if x is an integer but contains only one value? (R 1.5.0
on Linux i386)
Thanks, Sven
2001 Jun 01
1
Read to the end of a stream (package Rstreams)
Hi there,
I want to read all data from a binary file (with package Rstreams). I
tried the following function (small piece of code):
readdata <- function (file) {
# open file
s <- openstream(file, "read")
# read int values
while (s$position <= s$size) {
readint(s, 1, 2) -> type
[...]
}
return()
}
and got the following error:
> Error in
2003 Aug 27
2
How to test a model with two unkown constants
Hi all,
suppose I've got a vector y with some data (from a repeated measure
design) observed given the conditions in f1 and f2. I've got a model
with two unknown fix constants a and b which tries to predict y with
respect to the values in f1 and f2. Here is an exsample
# "data"
y <- c(runif(10, -1,0), runif(10,0,1))
# f1
f1 <- rep(c(-1.4, 1.4), rep(10,2))
# f2
f2 <-
2002 Sep 18
2
comparing elements in a vector
How about --
d is a vector as you described it, then define
d1 <- d[1:(length(d)-1)]
d2 <- d[2:length(d)]
d3 <- d1!=0 & d2!=0 & (d1-d2)!=0
I think d3 is the final vector you wanted.
I didn't test it, so I'm not sure my logic is right, but I think this is
definitely a path you can take...
-----Original Message-----
From: owner-r-help at stat.math.ethz.ch
2001 Jun 13
2
Maybe OT: large fonts in eps-figures
Hi there,
if I copy an x11() graphics device to an eps-file (with dev.copy2eps())
the font in the legend is very large and doesn't fit to the legend box
in the eps-file (same with a postscript file). I'm not sure if this is a
R problem rather than a ghostscript one.
But is there a way to solve this problem in R or depends this on my
ghostscript installation?
System: R Version 1.2.3 on
2001 Jul 09
1
Error plotting time series
Hi there,
when plotting a time series I got this error message:
Error in xy.coords(x, y, xlabel, ylabel, log) :
x and y lengths differ
In addition: Warning messages:
1: longer object length
is not a multiple of shorter object length in: vx.1[int.lr] + (-1) *
vx.2[int.rl]
vx.1[int.rl] and vx.2[int.rl] have the same length. Does anybody know
what the warning message means?
Thanks, Sven
2002 Mar 21
1
Wilks Lamba
Hi all,
I can't figure out how to compute Wilks Lambda in a one way repeated
measure design. My matrix looks like:
> t2.m
Blank ECR ENC UEA UED
1 -0.15 0.14 0.16 0.09 0.14
2 0.30 0.08 0.14 0.14 0.14
[...]
where each row is a case and the columns are levels of one factor (named
trial):
> t2.fit <- manova(t2.m ~ 1)
> summary(t2.fit, intercept=T,
2000 Apr 13
1
how convert an array to a matrix?
Dear R-list,
is it possible to convert a threedimensional array to a matrix? The
array has the form
, , -4
-4 -3 -2 -1 0 1 2 3 4
-80 1588 NA 1171 NA 1121 NA 852 NA 1580
0 1497 NA 1311 NA 762 NA 1185 NA 1466
80 1960 NA 1257 NA 941 NA 1435 NA 1636
[...]
I need a matrix like
-80 -4 -4 1588
-80 -3 -4 NA
-80 -2 -4 1171
etc.
I?ve tried write.table() with no succes:
2002 Oct 09
1
Large F-value and small P-value
Hi all,
I computed a Wilks Lambda Test with manova:
> df.man <- manova(df.mul ~ 1, na.rm=TRUE)
> summary(df.man, intercep=T, test="Wilks")
Df Wilks approx F num Df den Df Pr(>F)
(Intercept) 1 0.0002824 393.3 9 1 0.03911 *
Residuals 9
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*'
2001 Jul 26
6
replacing values in a vector
Hi all,
there is a vector v with several NAs. I want to create a new vector n of
the same length as v and the same NAs as in v and tried this:
n <- vector(length=length(v), mode="numeric")
replace(n, which(is.na(v)), NA)
but this does't work, all values in n are 0. What went wrong?
Thanks, Sven
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
2002 Jan 09
2
Creating subsets with factors
Hi all,
I don't understand the following output. I've created a data subset from
a data frame by
> p1.sub <- subset(p1.dat, vp!="p1")
this is ok. But
> attach(p1.sub)
> vp
[1] p1ab p1ab p1ab p1ab p1ab p1br p1br p1br p1br p1br p1kf p1kf p1kf
p1kf p1kf
[16] p1mg p1mg p1mg p1mg p1mg p1mw p1mw p1mw p1mw p1mw
Levels: p1 p1ab p1br p1kf p1mg p1mw
shows me that the
2000 Apr 10
2
Newbie: how to calculate group averagege?
Dear R-List,
I?m new in R, so I hope my question is not to primitive, but I haven?t
found a solution in the R-help. I?ve got a datatframe with 3 factors,
called xf (9 levels), yf (9 levels), zf (3 levels) and one variable
(rt):
> xf yf zf rt
> 1 1 1 67
> 1 1 1 56
> 1 1 1 60
[...]
> 1 1 2 58
> 1 1 2 61
[...]
> 9 1
2004 Oct 25
2
box() and hist()
Hi,
does anybody know why the following is not working:
> hist(rnorm(200))
> box(bty="o")
gives me a box without rounded corners.
System:
> R.Version()
$platform
[1] "i386-pc-mingw32"
$arch
[1] "i386"
$os
[1] "mingw32"
$system
[1] "i386, mingw32"
$status
[1] ""
$major
[1] "2"
$minor
[1] "0.0"
$year
2000 Jan 31
1
Is there a "by()" function in R like in S?
Hallo R-users,
I'm a R (release 0.65) novice, but have read some about S. Is there a
similar function in R like the S function "by()" ?
This function groups data typically from a matrix by a variable and
applies a function on the grouped data. Example in S:
> by(data, year, summary)
where "data" is a matrix containing the variable "year". After grouping
2004 Jul 16
1
Fixed and random factors in aov()
Hi,
I'm confused about how to specify random and fixed factors in an aov()
term. I tried to reproduce a textbook example: One fixed factor (Game, 4
levels) and one random factor (Store, 12 levels), response is Points.
The random factor Store is nested in Game. I tried
> str(kh.df)
`data.frame': 48 obs. of 4 variables:
$ Subj : Factor w/ 48 levels
2002 Jun 24
1
Covarianz matrix from an aov model
Try vcov() from the MASS library.
-----Original Message-----
From: Sven Garbade [mailto:garbade at psy.uni-muenchen.de]
Sent: Monday, June 24, 2002 11:19 PM
To: R-Help
Subject: [R] Covarianz matrix from an aov model
Hi all,
how can I extract the covariance matrix from an object of class aov?
Thanks, Sven
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
-.-
r-help
2010 Apr 09
2
computation of dispersion parameter in quasi-poisson glm
Hi list,
can anybody point me to the trick how glm is computing the dispersion
parameter in quasi-poisson regression, eg.
glm(...,family="quasipoisson")?
Thanks ®ards, Sven
2006 May 03
1
Inverse X'WX matrix from weighted linear regression
Dear list,
how can I compute the inverse of the X'WX matrix ("inverse of the weighted sum of squares and crossproducts matrix") from an object of class "lm" from a weigthed linear regression?
Thanks, Sven
2006 Dec 04
2
background color in strip.custom()
Hi all,
how can I change the background color in lattice strips according to a
factor level, eg:
library(lattice)
x <- rnorm(100)
y <- sqrt(x)
f <- gl(2, 50, c("A", "B"))
xyplot(y ~ x | f)
I like to change the background color of the strips according to the
levels in f and tried several things like this with no success:
xyplot(y ~ x | f,