On Tue, 11 Apr 2006, Brian Quinif wrote:
Perhaps someone will have a solution to my more general problem, but
here is the specific one:
These are issues discussed in all good books on R/S programming.
E.g. section 3.5 in MASS.
I used the round() function to round some estimates to 3 decimal
On Mon, 10 Apr 2006, Tarca, Adi wrote:
Hi all,
I have a dataset in which the output Y is observed on two groups of
patients (treatment factor T with 2 levels).
Every subject in each group is observed three times (not time points but
just technical replication).
I am interested in
Are you looking for trunc.POSIXt? That does the -1 direction, and its
code could easily be adapted (or just add one unit).
I don't see what this has to do with 'align'.
On Mon, 10 Apr 2006, Omar Lakkis wrote:
I use POSIXct for datetimes. Is thee a timeAlign function that I can
use where :
2006/4/10, Tarca, Adi [EMAIL PROTECTED]:
Hi all,
I have a dataset in which the output Y is observed on two groups of
patients (treatment factor T with 2 levels).
Every subject in each group is observed three times (not time points but
just technical replication).
I am interested in
On Mon, 10 Apr 2006, R. A. L. Carter wrote:
Hi Folks,
I sent this message before but I don't think it got through because I
wasn't registered.
It did appear:
https://stat.ethz.ch/pipermail/r-help/2006-April/092079.html
(you could have checked the archives as easily as I did).
We didn't
On Mon, 10 Apr 2006, Dave Armstrong wrote:
Dear R-help Listers,
I am curious if there is some (hopefully easy) way to change the number of
characters that can be converted to a single string via any of deparse,
as.character or toString. It seems that the limit is 500 for all of
these.
Hello,
I have about 2000 data files which I want to analyse.
The file names are all very similar =
pvariable1_tvariable2_cvariable3.txt
There aren´t so much different variables (about 70) , just different
combinations of them.
To allow an easy way of handling those data files I was wondering
I am using using the latex() function in the Hmisc package. I want to
create some table captions that are relatively long and need to take
up more than one line.
When I do this, the LaTeX output has the majority of the text on one
line with the remainder spilling over into the second line as if
*on Windows*
This is how the Windows file system works. You cannot have a.png and
A.png in a directory on a Windows file system. If you unlink() the file
with the clashing name it will create it with the case you want (but since
asking Windows for B.png will give you b.png, does it matter?).
Timo Becker wrote:
Dear R users,
when I create files with the png graphics device, it overwrites existing
files which differ concerning upper/lower case file names.
For example, I want to create four files named a.png, A.png, b.png
and B.png which contain the corresponding letters:
Dear R users,
when I create files with the png graphics device, it overwrites existing
files which differ concerning upper/lower case file names.
For example, I want to create four files named a.png, A.png, b.png
and B.png which contain the corresponding letters:
png(a.png)
plot(1, type=n)
Prof Brian Ripley schrieb:
*on Windows*
This is how the Windows file system works. You cannot have a.png and
A.png in a directory on a Windows file system. If you unlink() the file
with the clashing name it will create it with the case you want (but
since asking Windows for B.png will
Dear R users,
I´m student of Master in Statistic and Data analysis, in New University of
Lisbon. And now i´m writting my dissertation in variance estimation.So i´m
using Survey Package to compute the principal estimators and theirs variances.
My data is from Incoming and Expendire Survey.
'letters' is an argument to your function, with default 'letters', that it
itself. (Default arguments are evaluated in the evaluation frame of the
function, not in the caller.)
Try letters=base::letters in your argument list.
On Tue, 11 Apr 2006, Matthieu Dubois wrote:
Dear useRs,
I have
Carlos Creva Singano (M2004078) wrote:
1. How to compute Bootstrap and Jackknife Bias of estimates, like mean?
Have you had a look at packages boot and bootstrap? E.g. you can
compute the bias and s.e. of an estimate theta using bootstrap
library(boot)
a - boot( data, theta, R=1000 )
where
Hi there.
Have a look at this. Want an ipod nano?
http://mygift321.com/?r=MTGERHcEFhJJAWcKAy4Mi=gmailz=1[EMAIL PROTECTED]tc=21
Cheers
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Or use letters. = letters and similarly for the other arguments of the
form x = x.
On 4/11/06, Prof Brian Ripley [EMAIL PROTECTED] wrote:
'letters' is an argument to your function, with default 'letters', that it
itself. (Default arguments are evaluated in the evaluation frame of the
I don't know wether R does variable substitution. In order to help
solve the problem, I suggest:
- you might use cat in order to concatenate the files (if this fits in
the problem)
- you might write a shell or perl script which would write an .R file
by doing the required looping, given the three
Hi All,
I'm running garch models for different combinations, like (1,1),(1,2) etc.
in a for loop. But, when i'm running it, R is showing a text *
ESTIMATION WITH ANALYTICAL GRADIENT * . How can delete it? Any options
for this?
--
SUMANTA BASAK.
[[alternative HTML version
Here are some different approaches:
x - 1
y - 2
paste(p, x, y, .dat, sep = ) # p12.dat
sprintf(p%s%s.dat, x, y)
file.path(a, b.dat) # a/b.dat
library(gsubfn)
gsubfn(,,p$x$y.dat) # p12.dat
On 4/11/06, Damaris Zurell [EMAIL PROTECTED] wrote:
Hello,
I have about 2000 data files which I
Dear list members
I have been using RWinEdt for an extended time and on different machines
(all under Windows XP). There is a strange problem that comes up
sometimes on one of them, always on another and never on a third one.
(However, I do think I have followed the same procedures on all
Dear R-users,
I am wondering if anyone can tell me how Garch coefficients are calculated in
R. What is algorithm for that? If anyone give me a detail desceiption I will be
very grateful.
Thanks and Regards,
stat74
thanks in advance
Try:
junk - capture.output(print(abc))
On 4/11/06, Sumanta Basak [EMAIL PROTECTED] wrote:
Hi All,
I'm running garch models for different combinations, like (1,1),(1,2) etc.
in a for loop. But, when i'm running it, R is showing a text *
ESTIMATION WITH ANALYTICAL GRADIENT * . How can
Christian Gold wrote:
Dear list members
I have been using RWinEdt for an extended time and on different machines
(all under Windows XP). There is a strange problem that comes up
sometimes on one of them, always on another and never on a third one.
(However, I do think I have followed the
Brian Quinif wrote:
I am using using the latex() function in the Hmisc package. I want to
create some table captions that are relatively long and need to take
up more than one line.
When I do this, the LaTeX output has the majority of the text on one
line with the remainder spilling over
On Tue, 2006-04-11 at 07:54 -0500, Frank E Harrell Jr wrote:
Brian Quinif wrote:
I am using using the latex() function in the Hmisc package. I want to
create some table captions that are relatively long and need to take
up more than one line.
When I do this, the LaTeX output has the
Dear all,
I have a result my experiment like this below (here my toy example):
foo1 - list()
foo1[[1]] - c(10, 20, 30)
foo1[[2]] - c(11, 21, 31)
foo2 - list()
foo2[[1]] - c(100, 200, 300)
foo2[[2]] - c(110, 210, 310)
foo3 - list()
foo3[[1]] - c(1000, 2000, 3000)
foo3[[2]] - c(1100, 2100, 3100)
I got strange error messages when I had other packages opened the same
time as lme4 (e.g. as nlme).
Try to not have other libraries loaded.
/A
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Are you trying to turn each of the components of foo into a matrix?
If that's it then:
lapply(foo, function(x) matrix(unlist(x), nrow = length(x)))
On 4/11/06, Muhammad Subianto [EMAIL PROTECTED] wrote:
Dear all,
I have a result my experiment like this below (here my toy example):
foo1 -
Brian Quinif bquinif at gmail.com writes:
I used the round() function to round some estimates to 3 decimal
places. I then sent put the rounded estimates in a matrix and used
latex() to make a LaTeX table from them. However, in my table, there
are estimtes which only have 2 decimal places.
I forgot the byrow = TRUE. Try this instead:
lapply(foo, function(x) matrix(unlist(x), nrow = length(x), byrow = TRUE))
On 4/11/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:
Are you trying to turn each of the components of foo into a matrix?
If that's it then:
lapply(foo, function(x)
stat stat wrote:
Dear R-users,
I am wondering if anyone can tell me how Garch coefficients are calculated
in R. What is algorithm for that? If anyone give me a detail desceiption I
will be very grateful.
Thanks and Regards,
stat74
Goto www.rmetrics.org and download the
R-help,
I'm using a lm model for some data.
The code is below:
fitData - lm(formula = log(f1) ~ year + mon + pair + rek, weights = w,
data = data)
I want to get the fitted values of log(f1) for every year so I do:
tapply(fitted(fitData), list(data96$year), mean)
but there is no parameter
Hi,
On a new pc, I am trying to 'CHECK' an R package
containing only R code but the process is hanging and I cant see why.
The 00check.log file shows
* using log directory 'c:/temp/opRisk.Rcheck'
* using R version 2.2.1, 2005-12-20
* checking for file 'opRisk/DESCRIPTION' ... OK
* this is
Hi
you should use the package boot to compute bootsrap and jacknife
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
On Tue, 11 Apr 2006, Renaud Lancelot wrote:
2006/4/10, Tarca, Adi [EMAIL PROTECTED]:
Hi all,
I have a dataset in which the output Y is observed on two groups of
patients (treatment factor T with 2 levels).
Every subject in each group is observed three times (not time points but
just
If I understand what you want correctly:
foo1 - list(c(10,20,30), c(11,21,31))
foo2 - list(c(100,200,300), c(110,210,310))
foo3 - list(c(1000,2000,3000), c(1100,2100,3100))
fool - list(foo1, foo2, foo3)
lapply(fool, function(x) do.call('rbind', x))
[[1]]
[,1] [,2] [,3]
[1,] 10
On Tue, 11 Apr 2006, Carlos Creva Singano (M2004078) wrote:
Dear R users,
I´m student of Master in Statistic and Data analysis, in New University
of Lisbon. And now i´m writting my dissertation in variance
estimation.So i´m using Survey Package to compute the principal
estimators and theirs
On Tue, 11 Apr 2006, justin bem wrote:
Hi
you should use the package boot to compute bootsrap and jacknife
Not if he wants to get the right answers. This was for a stratified
multistage cluster-sampled survey.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
On Tue, 11 Apr 2006, Thomas Lumley wrote:
On Tue, 11 Apr 2006, Renaud Lancelot wrote:
2006/4/10, Tarca, Adi [EMAIL PROTECTED]:
Hi all,
I have a dataset in which the output Y is observed on two groups of
patients (treatment factor T with 2 levels).
Every subject in each group is observed
On Tue, 11 Apr 2006, Prof Brian Ripley wrote:
On Tue, 11 Apr 2006, Thomas Lumley wrote:
He has a linear model with the same number of observations for each person
Not so: some have 3 and some have 2, and the two levels of T are not quite
balanced (29/28).
and no covariates that vary
I am just starting with R and I am having a stupid
problem understanding an if else statement because I
am thinking in terms of something like SAS or Fortran.
I am completely misunderstanding what I am reading in
the Intro to R, the R-Language Definitions and the
Help Mailing List. I have
Thank you very much for your useful suggestions.
These are exactly what I was looking for.
foo - list(foo1, foo2, foo3)
lapply(foo, function(x) matrix(unlist(x), nrow = length(x), byrow = TRUE))
or
lapply(foo, function(x) do.call('rbind', x))
Best, Muhammad Subianto
On 4/11/06, Muhammad
Dear All,
When working with composite plots many times I have to manually adjust the
graphics display window size to get the best results (in Windows XP). I can
then find the size of the adjusted window in inches using par(din).
Now I would like to export the adjusted graphs as png's and have
John Kane wrote:
I am just starting with R and I am having a stupid
problem understanding an if else statement because I
am thinking in terms of something like SAS or Fortran.
I am completely misunderstanding what I am reading in
the Intro to R, the R-Language Definitions and the
Dear All,
When working with composite plots many times I have to manually adjust the
graphics display window size to get the best results (in Windows XP). I can
then find the size of the adjusted window in inches using par(din).
Now I would like to export the adjusted graphs as png's and have
Thank you very much for your useful suggestions.
These are exactly what I was looking for.
foo - list(foo1, foo2, foo3)
lapply(foo, function(x) matrix(unlist(x), nrow = length(x), byrow = TRUE))
or
lapply(foo, function(x) do.call('rbind', x))
Best, Muhammad Subianto
On 4/11/06, Muhammad
On Tue, 11 Apr 2006, Christos Hatzis wrote:
Dear All,
When working with composite plots many times I have to manually adjust the
graphics display window size to get the best results (in Windows XP). I can
then find the size of the adjusted window in inches using par(din).
Now I would like
Dear all
I need to run an anova from a factorial model
y_{ijk}=\alpha_i+\beta_j+(\alpha\beta)_{ij}+e_{ijk}
and calculate type II and III sums of square, but I have an empty
cells, so anova function from package car fail. (I believe)
y-c(7,13,6,10,8,11,8,3,7,5,65)
On Tue, 11 Apr 2006, Prof Brian Ripley wrote:
On Tue, 11 Apr 2006, Christos Hatzis wrote:
Dear All,
When working with composite plots many times I have to manually adjust the
graphics display window size to get the best results (in Windows XP). I can
then find the size of the adjusted
Dear Mario,
When there are empty cells, some of the usual contrasts defining
interactions (and possibly even main effects) can't be estimated, and one
should be especially careful about what hypotheses are being tested. There's
a good introduction to these problems (in my opinion) by Searle,
One other thought. If they are all of the same dimension you could
alternately consider putting them into a 3d array:
library(abind)
abind(lapply(foo, function(x) do.call(rbind, x)), along = 3)
which may or may not have some advantage to you.
On 4/11/06, Muhammad Subianto [EMAIL PROTECTED]
Hi there, I encountered a weird problem using cph()
with Design package:
I have 2 datasets, say dat1 and dat2, both data
frames with 3 columns time,status and scores,
all numeric
If I run the following:
dd-datadist(dat1)
options(datadist='dd')
dd
time status scores
sorry this problem only occurs in S-Plus, not in R.
--- array chip [EMAIL PROTECTED] wrote:
Hi there, I encountered a weird problem using cph()
with Design package:
I have 2 datasets, say dat1 and dat2, both data
frames with 3 columns time,status and scores,
all numeric
If I run the
Hello.
I have n files in a directory: file1, ..., filen.
I read them with the following commands:
list=scan(file=list,what=list(nom=))
# in the file list, I have all the filenames.
n=length(list[[1]])
for (i in 1:n){
aux - paste(p,i,sep=)
assign(aux,
On 4/11/2006 2:30 PM, Arnau Mir wrote:
Hello.
I have n files in a directory: file1, ..., filen.
I read them with the following commands:
list=scan(file=list,what=list(nom=))
# in the file list, I have all the filenames.
n=length(list[[1]])
for (i in 1:n){
aux - paste(p,i,sep=)
From version 2.01-17, splancs, a package for spatial point pattern
analysis, depends on package sp.
Because both packages have defined a bounding box function bbox(), the
bbox() method in sp has been given priority, and the default object for
dispatch is an array with two or more columns (as in
Actually it does.
Thank you, Prof. Ripley.
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Tuesday, April 11, 2006 1:26 PM
To: Christos Hatzis
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Graphics device size
On Tue, 11 Apr 2006, Prof Brian Ripley wrote:
Dear All,
Probably it is not highly relevant question: Why do stepwise regression
functions in R (step() or stepAIC()) add/delete categorical variables as a
set? For example, I have a four-level factor variable d, so dummies are
d1,d2,d3, as stepwise regression operates d, adding or removing,
Mike Wolfgang wrote:
Dear All,
Probably it is not highly relevant question: Why do stepwise regression
functions in R (step() or stepAIC()) add/delete categorical variables as a
set? For example, I have a four-level factor variable d, so dummies are
d1,d2,d3, as stepwise regression operates
I am using R for prediction. I want to use VARIMA. I could not not
find any package which would provide multivariate ARIMA that is
VARIMA. I looked at the list of CRAN packages and found DSE package
and also fracdiff package but they do not contain VARIMA modeling
function.
DSE had only VAR
Hi all --
So I'm working through my statistics homework again, and trying to
reproduce the examples in the book (Kirk's _Experimental Design_,
third edition) in R. This is a completely randomized hierarchical
design (CRH-28(A)). The B factor is completely nested within the A
factor. Pages
Thanks for the numerous responses. I was asked to post the solution so that
others could see it.
The solution is to use the estimable function from the gregmisc or gmodels
package. First you create a matrix (the cm matrix) of the particular linear
combinations of parameters you are concerned
Farrel Buchinsky [EMAIL PROTECTED] wrote:
2)Storing all the SNP data as a string seems quite clever and a space-saving
way of doing it. However, if you were to analyze a whole chromosome at a
time you would still be creating one almighty big table albeit only
temporarily. Do you use R to run
On 4/11/06, Fredrik Karlsson [EMAIL PROTECTED] wrote:
Dear list,
I used to have this code in a Sweave document:
sel - placeDF$Place == Velar placeDF$manner == Plosive
table(placeDF$agem[sel], placeDF$place[sel]) - pd
print(
dotplot(
prop.table(pd,margin=1),
Do some of your observations have weights=0? If yes, that might be
the problem. I don't know the details, but some functions like glm have
an argument 'subset' and might give inappropriate answers with weights =
0.
Beyond that, I just got 69 hits for RSiteSearch(glmmPQL
RSiteSearch(multicomp) produced 42 hits for me just now, including
mentions of a multcomp package (with no i) and functions TukeyHSD,
pairwise.t.test and pairwise.wilcox.test in ctest, including a comment
from Ripley that, If you look in the archives you will see this has
cropped up
Trying to define group methods for the expression class, running
into some problems.
setMethod(Arith, signature(e1 = expression, e2 = expression),
+ function(e1, e2) 5)
[1] Arith
expression(2)+expression(5)
Error in expression(2) + expression(5) : non-numeric argument to
binary operator
I have found out that the way to have a break apart a long caption in
the way that one desires is to have a short caption for use in the
list of tables that is not broken with \\
I can manually adjust the LaTeX output from the Hmisc function latex()
to get the tables how I want them. However, I
Try this way
i - 1
if (i==1) {years - Freshmen}
if (i==1) {years - Sophomores}
it should work!
Brian Quinif wrote:
I am ashamed to be asking this question, but I couldn't find the
solution anywhere. Searching for if and R is not very
productive...
I cannot get a simple if statement to
70 matches
Mail list logo