You have forgotten drop=FALSE in mat[1,]: take a look at what it produces.
See ?Extract for more details.
On Mon, 3 Jan 2005, Patrick Giraudoux H wrote:
Would a specialist of the point or plot functions try the following:
mat-matrix(c(1.836767,4.025989,0.6396777,0.376),ncol=2)
plot(mat)
Patrick Giraudoux H wrote:
Hi,
Would a specialist of the point or plot functions try the following:
mat-matrix(c(1.836767,4.025989,0.6396777,0.376),ncol=2)
plot(mat)
points(mat[1,],col=red)
..A lag appears on x for mat [1,1] between the two displays.
I wonder if this example may be due to a
Hello,
R exits every time I want to fit a GEE with AR-1 correlation
structure. It just closes itself without any error.
My Data has about 3300 observation (using the big data with
about 20 000 observation it doesn't work at all). I'm a
Windows 2000 User of R 2.0.1 but it was the same problem on
Dear Jim,
many thanks for your reply and support.
It seems to be that with your help I could solve my problem with the
plotting of the data. The only thing that does not work is to see the
coloured lines, maybe because of the crowd of curves.
With the density of curves I ment to distinguish in
UweL == Uwe Ligges [EMAIL PROTECTED]
on Fri, 31 Dec 2004 08:51:06 +0100 writes:
UweL Jeroen Van Houtte wrote:
To display many small graphics one above the another
(like Edward Tufte's sparklines), I'd like to use layout
in the graphics package with more than 15 rows.
In S-Plus, I can look at the structure of a function (for example, hist)
simply by entering
hist RETURN
however, if I do this in R, I get the response
function (x, ...)
UseMethod(hist)
environment: namespace:graphics
How can I inspect the structure of a function in
There's a new page on the Burns Statistics website
http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html
that looks at spreadsheets from a quality assurance perspective. It
presents R as a suitable alternative to spreadsheets. Also there are
several specific problems with Excel that
Package gee does that -- it is rather old and abandoned by its author a
long time ago. Perhaps you could try yags or geepack instead?
On Mon, 3 Jan 2005, Johanna BRANDT wrote:
R exits every time I want to fit a GEE with AR-1 correlation structure. It
just closes itself without any error.
My
Richard Dybowski wrote:
In S-Plus, I can look at the structure of a function (for example, hist)
simply by entering
hist RETURN
however, if I do this in R, I get the response
function (x, ...)
UseMethod(hist)
environment: namespace:graphics
How can I inspect the
Richard,
There are usually several methods available - depending on the class of
the object you apply 'hist' on.
In this case two are 'non-visible'
##
methods(hist)
[1] hist.Date* hist.default hist.POSIXt*
Non-visible functions are asterisked
##
Probably you are interested in
On Mon, 3 Jan 2005, Richard Dybowski wrote:
In S-Plus, I can look at the structure of a function (for example, hist)
simply by entering
hist RETURN
however, if I do this in R, I get the response
function (x, ...)
UseMethod(hist)
environment: namespace:graphics
How can
On Mon, 03 Jan 2005 10:50:22 +, Richard Dybowski
[EMAIL PROTECTED] wrote:
In S-Plus, I can look at the structure of a function (for example, hist)
simply by entering
hist RETURN
however, if I do this in R, I get the response
function (x, ...)
UseMethod(hist)
Christoph Buser wrote:
Hi all
I tried to reproduce an example with lme and used the Orthodont
dataset.
library(nlme)
fm2a.1 - lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 | Subject)
anova(fm2a.1)
numDF denDF F-value p-value
(Intercept) 180 4123.156 .0001
age
Hello list (happy new yeaR),
Here's a copy of a message i just send to Christian Hennig (who wrote
the fpc package).
That may interrest some of you, and maybe someone could have a better
solution than mine.
Romain.
The dataset lutenhorm in package bootstrap has four
time series, while the book has only one. The one from the book is
V4. V1 is just the time indeces. Anybody knows wht is
V2, V3, V5
?
Kjetil
--
Kjetil Halvorsen.
Peace is the most effective weapon of mass construction.
-- Mahdi
Hi all -
R2.0.1 on OSX;MASS library;nlme library
I am trying to emulate the solution to a problem set that has normally
been run in Genstat, using R. The problem that I am having at the
moment is with the following glmm question (using glmmPQL from the MASS
library):
We have two different
Dear List:
I am having to build a block-diagonal matrix (vl) and am currently using
the following code.
I-diag(sample.size)
vl-kronecker(I,vl.mat)
This code works fine, but for large N, it is a huge memory hog. Is there
a more efficient method for constructing vl?
Thanks,
Harold
Hello,
My name is Sivan and I am a master degree student in statistics,my problem is
as follows:
I have a dataset containing gene sequences and I would like to create a
phylogenetic tree from it.
The problem that I can't seem to find a function to do this kind of operation.
I read the ape
Douglas Bates wrote:
Christoph Buser wrote:
Hi all
I tried to reproduce an example with lme and used the Orthodont
dataset.
library(nlme)
fm2a.1 - lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 |
Subject)
anova(fm2a.1)
...
Regards,
Christoph Buser
No. The calculation of denominator
Hello:
I am attempting to use R to analyze amino acid frequencies in aligned
protein sequences and need some help. So far, I have imported my sequence
alignment into a data frame (lets call it alignment) with each site in one
column, so that I have a data frame consisting of columns of letters
Hi Silvan - I'd love to proven wrong, but I think unless
you're interested in a cluster dendrogram
(neighbour-joining or UPGMA etc.) that the most standard
phylogeny-generation algorithms aren't implemented as R
functions (yet). There are very efficient third party
programs for tree
On Mon, 3 Jan 2005, Doran, Harold wrote:
Dear List:
I am having to build a block-diagonal matrix (vl) and am currently using
the following code.
I-diag(sample.size)
vl-kronecker(I,vl.mat)
This code works fine, but for large N, it is a huge memory hog. Is there
a more efficient method for
Estarei de férias no período de 20 de dezembro à 03 de fevereiro de 2005.
Durante minha ausência, favor procurar Josi ([EMAIL PROTECTED]) 3235-5109.
Obrigada.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
My best advice:
Yours is a complicated question, but it is probably the wrong question. I
suspect that you are trying to reinvent wheels: I believe there has been a
lot of work in the statistical community on (OK, maybe DNA not amino acid)
sequence alignment. Undoubtedly much more published in
Hi,
I had the same opportunity and ran into the same problem. The error
occurred when trying to read 25 CEL-files with the ReadAffy-package and
happened in R1.9 as well as R2.0. The problem seemed to be due to an error
in the Tcl/Tk package and might relate to the installation of Tcl/Tk on
our
On Sat, 1 Jan 2005, Gabor Grothendieck wrote:
Look at the source of expand.model.frame.
The second line of expand.model.frame uses the data component of mod$call
and if its not there expand.model.frame chokes. As a workaround check
if its missing and supply it yourself:
if
Hi there,
I was wondering if you could help me to write the following procedure. I
need to plot x against y.
The x is read directly from a file. The y is delta from power.t.test:
y - as.vector(power.t.test(n=4, sd =z, power = 0.8)$delta), where z
values are also read from the file.
Thank you for
Hi,
I run into the following problem:
data(iris)
as.array(iris[,-5])
Error in dimnames-.data.frame(`*tmp*`, value = list(c(Sepal.Length, :
invalid dimnames given for data frame
as.array(as.matrix(iris[,-5]))
When trying to convert numeric data frame to an array. Conversion to matrix
works
Happy new year to all;
A few days ago, I posted similar problem. At that time, I found out that our
R program had been 32-bit compiled, not 64-bit compiled. So the R program
has been re-installed in 64-bit and run the same job, reading in 150
Affymetrix U133A v2 CEL files and perform dChip
Have you checked whether there are limits set? What does `ulimit -a' say?
Do you know how much memory the R process is using when the error occurred?
We've had R jobs using upwards of 13GB on a box with 16GB of RAM (SLES8 on
dual Opterons) and never had problems.
Andy
From: Tae-Hoon Chung
Thanks Peter and Andy;
I just found it was not due to memory problem. It was false alarm ...
64-bit compiled program works fine!
On 1/3/05 3:39 PM, Peter Dalgaard [EMAIL PROTECTED] wrote:
Tae-Hoon Chung [EMAIL PROTECTED] writes:
Happy new year to all;
A few days ago, I posted similar
To anyone who can help:
I have two brief questions concerning sapply. Following below is the
code for my example. The two problems are described at the end of the
code:
site - rep(2:6, each = 12)
tillage - rep(c(1,-1), each = 6, times = 5)
carbon - c(18.23, 16.06, 17.81, 16.07, 17.26,
Philip,
1) right now, tillage is being defined as a numeric vector and then
being coerced into a factor. So, try defining tillage in this way:
tillage - factor(rep(c(1,-1), each = 6, times = 5), levels=c(1,-1))
2) Why do you want to do this?
I hope this helps,
Andrew
On Mon, Jan 03,
Hi, Christoph:
As documented in Pinheiro and Bates (2000) Mixed-Effects Models
for S and S-Plus (Springer), the nlme package includes a function
simulate.lme. They used that function to study the performance of the
likelihood ratio statistic under the null hypothesis of no effect. The
Try something like:
nr.carbon.2 - tapply(carbon, list(tillage, site), mean)
nr.carbon.2
2345 6
-1 14.61667 18.52667 27.83833 21.42167 20.085
1 17.08500 19.76000 32.87000 24.70833 24.120
as.vector(nr.carbon.2[2:1,])
[1] 17.08500 14.61667 19.76000
IT,IT
ITHRCTOIT
IThttp://www.hbitjob.org
http://www.51itjob.org/bbs
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
I have a problem building an extension using ISNAN() on R version 2.0.x.
In R 1.9.1 Arith.h and Rmath.h contained code like
#ifdef IEEE_754
# define ISNAN(x) (isnan(x)!=0)
#else
# define ISNAN(x) R_IsNaNorNA(x)
#endif
#define R_FINITE(x)R_finite(x)
int R_IsNaNorNA(double);
int
On Tue, 4 Jan 2005, Bill Northcott wrote:
I have a problem building an extension using ISNAN() on R version 2.0.x.
In R 1.9.1 Arith.h and Rmath.h contained code like
#ifdef IEEE_754
# define ISNAN(x) (isnan(x)!=0)
#else
# define ISNAN(x) R_IsNaNorNA(x)
#endif
#define R_FINITE(x)
Hi
I'd like to know if the R can do the DCCA? Because I can't find the package
about
the DCCA in R. If it can not, please introduce a free software to me, which you
think can easily do the work , thanks!
jeff
I see you are quite correct that IEEE_754 is defined in Rconfig.h on
MacOS X. However, I was building against a standalone libRmath v1.9.1.
So I was including Rmath.h only and IEEE_754 was not defined. The
result was that I got R_IsNaNorNA in the preprocessed source.
I finally found it.
40 matches
Mail list logo