Could Petr Pikal's peaks function
(http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html) be of any
use?
Ulrik
On 7/24/06, Tauber, Dr E. [EMAIL PROTECTED] wrote:
Dear R-users,
We are monitoring the activity of animals during a few days period. The
data from each animal (crossing of
I am trying to use nlme but instead of using one of the identity variance or
covariance matrixes such as compsymm or ar1. Instead I want the covariance
matrix to be represented in the following manor. Is it possible to define my
own covariance matrix?
I have search and found papers saying I
On Monday 24 July 2006 20:16, Giuseppe Pagnoni wrote:
Dear all,
First of all I apologize if you received this twice: I was checking the
archive and I noticed that the text was scrubbed from the message,
probably due to some setting in my e-mail program.
I am unsure about how to specify a
Hi , hi all,
Am I correct to understand from the previous discussions on this topic (a
few years back) that if I have a matrix with missing values my PCA options
seem dismal if:
(1) I dont want to impute the missing values.
(2) I dont want to completely remove cases with missing
Nair, Murlidharan T wrote:
Can I run jobs in the background and the check the status of it from
time to time in Windows version of R?
Depends on your version of Windows. I am running the automatical
building of R packages in the background on a Windows Server 2003.
Unfortunately, if you are
This is a very late rply to this topic, but I had a similar problem and
I got different results with using ps.options() then using
postscript(...) and the one with ps.options produced the results
expected. The problem were the size of the graph and the paper and the
rotation.
If you want to know
Dear R-helpers:
My question is how do I efficient and valid correct for multiple tests in a
repeated measurement design:
Suppose we measure at two distinct visits with repeated subjects a treatment
difference on the same variable.
The treatment differences are assessed with a mixed model and
On 25-Jul-06 Ulrik Stervbo wrote:
Could Petr Pikal's peaks function
(http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html)
be of any use?
Ulrik
On 7/24/06, Tauber, Dr E. [EMAIL PROTECTED] wrote:
Dear R-users,
We are monitoring the activity of animals during a few days
period.
Good morning,
in a 3D plot based on persp() the axis title (of dimension z) overlaps with
the axis labels.
How can the distance (between axis labels and axis title) be increased?
I would be grateful if anybody could help me.
Kilian
Hi list,
I am interested in cluster analysis of microarray data. The data was generated
using cDNA method and a loop design.
I was wondering if any one has a suggestion about which package I can use to
analyse such data.
Many thanks in advance
Mahdi
--
---
Dear R users,
seqinR 1.0-5 has been released yesterday on CRAN, so that the source code
of the package should be available on all CRAN mirrors within the next 24h.
The updated package vignette is here:
http://pbil.univ-lyon1.fr/software/SeqinR/seqinr_1_0-5.pdf
User level visible changes are:
I tried odfWeave to create an OpenOffice file and found that it
exhausted the memory of my large linux machine and took a long time to
run. LaTeX with Sweave is blazing fast and extremely flexible. Most of
the time I can give clients a pdf file. Sometimes I'd like to make
Sweave LaTeX files
On Tue, 25 Jul 2006, Uwe Ligges wrote:
Nair, Murlidharan T wrote:
Can I run jobs in the background and the check the status of it from
time to time in Windows version of R?
Depends on your version of Windows. I am running the automatical
building of R packages in the background on a
On Mon, 24 Jul 2006, johan Faux wrote:
I have a question about deparse function in R
What is the reason that deparse use an argument like width.cutoff ?
Why the maximum cutoff is 500?
I was manipulating an R formula and used deparse. Since the length of
user's formula was greater then
Spencer Graves spencer.graves at pdf.com writes:
comments in line
Grathwohl, Dominik, LAUSANNE, NRC-BAS wrote:
Dear r-helpers,
I have a question about multiple testing.
Here an example that puzzles me:
All matrixes and contrast vectors are presented in treatment contrasts.
Hi Mixed Modelers,
I was wondering if there are citations relevant to the concerns that
Bates and others have regarding the inappropriateness of significance
tests in mixed models that use comparison F-ratios to theoretical F
distributions. It would just make my life a little easier
Using regular expression matching for this case may be overkill (the RE
engine will be doing a lot of backtracking looking at a lot of
non-matches). Here is an alternative that splits the text into a vector
of words, extracts the last 2 letters of each word (remember if the last
3 letters match,
Dear all,
Based on my question a few months ago
https://stat.ethz.ch/pipermail/r-help/2006-January/086952.html
and solved with
https://stat.ethz.ch/pipermail/r-help/2006-January/086955.html
https://stat.ethz.ch/pipermail/r-help/2006-January/086956.html
and from
Frank E Harrell Jr f.harrell at vanderbilt.edu writes:
Has anyone
been able to get tth to work with Sweave, defining Sweave.sty and other
needed .sty files to be accepted by tth?
Thanks.
I use tth and ttm (Tex to MathML) -- not for Word users,
just for easier web posting.
I just use
Regarding having to do a lot of backtracking one can just
look at the relative comparison of speeds and we see
that they are comparable in speed.
In fact the bottleneck is not the backtacking but strapply.
I had coded the regexp version for compactness of code but if we replace
the strapply with
I just sent this to you in a personal response, but for purposes of
archives, the following is one reference:
@book{mccu:sear:2002,
author ={Charles E. McCulloch and Shayle Searle},
year={2002},
title ={Generalized, Linear,
Dear Kurt,
I've uploaded a new version (0.9-4) of the sem package to CRAN.
I hope that everything is well with you.
John
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox
Dear all,
Suppose h1.lm1 is my multiple regression object. When I tried to use
plot(h1.lm1), I got the following warning message.
-
plot(hla4.lm1)
Warning messages:
1: NaNs produced in: sqrt(crit * p * (1 - hh)/hh)
2: NaNs produced in: sqrt(crit * p * (1 - hh)/hh)
Before comparing times we should make sure that they functions return
the same thing. My original function (f1 below) labels the potential
rymes with match numbers as well as finding possible rymes, if you just
want the r flag then the for loop can be eliminated giving f4 as
follows:
f4 -
Ben Bolker wrote:
Frank E Harrell Jr f.harrell at vanderbilt.edu writes:
Has anyone
been able to get tth to work with Sweave, defining Sweave.sty and other
needed .sty files to be accepted by tth?
Thanks.
I use tth and ttm (Tex to MathML) -- not for Word users,
just for easier web
Hi,
I've got a problem with RODBC and saving (sqlSave) of a dataframe in Access.
R 2.0.1 is running on windows XP.
When executing the examples in R help for the USArrests data set sqlSAve
works fine, but running sqlSave() for a dataframe Adat
str(Adat)
`data.frame': 1202 obs. of 18
Hi All,
I sent the following message to R-help on July 14th, 2006:
Let's say I have the following formula:
a.formula - x ~ y + z
I want to extract the left and right-hand sides of the function so
that I have two character vectors like the ones you would create using
the following assignments:
Hi
Kilian Plank wrote:
Good morning,
in a 3D plot based on persp() the axis title (of dimension z) overlaps with
the axis labels.
How can the distance (between axis labels and axis title) be increased?
Unfortunately, there is not the same level of control that is available
for 2D plots,
Hello,
I want to have a title that will look something like:
Results for \theta=2.1, given that I have a variable theta=2.1, and
\theta should show on the screen like the greek letter.
I've tried a lot of things:
theta - 2.1
plot(1:10, main=expression(paste(Results for, theta, =, eval(theta
Le Mardi 25 Juillet 2006 15:48, Daniel Gerlanc a écrit :
Hi All,
I sent the following message to R-help on July 14th, 2006:
Let's say I have the following formula:
a.formula - x ~ y + z
I want to extract the left and right-hand sides of the function so
that I have two character vectors
This'll work.
theta - 2.1
plot(NA, xlim=c(0,1), ylim=c(0,1), xlab=bquote(theta == .(theta)),
ylab=bquote(theta == .(theta)), main=bquote(paste(Results for ,theta ==
.(theta
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
On Tue, 2006-07-25 at 13:12 -0700, Adrian Dragulescu wrote:
Hello,
I want to have a title that will look something like:
Results for \theta=2.1, given that I have a variable theta=2.1, and
\theta should show on the screen like the greek letter.
I've tried a lot of things:
theta - 2.1
Try:
plot(1:10, main = bquote(Results for ~ theta == .(theta)))
On 7/25/06, Adrian Dragulescu [EMAIL PROTECTED] wrote:
Hello,
I want to have a title that will look something like:
Results for \theta=2.1, given that I have a variable theta=2.1, and
\theta should show on the screen like
Paul Murrell p.murrell at auckland.ac.nz writes:
Hi
Kilian Plank wrote:
Good morning,
in a 3D plot based on persp() the axis title (of dimension z) overlaps with
the axis labels.
How can the distance (between axis labels and axis title) be increased?
Paul
Another way to
Can anyone share experience with opteron versus the xeon (woodcrest) for R
under linux? I am looking at using 16-32Gb of ram in a workstation (as
opposed to a server).
Thanks in advance
Sean
__
R-help@stat.math.ethz.ch mailing list
Hi everone,
I have a question on using lme on a mixed effects model.
The linear mixed model is in the form of:
y = bX +Zu + e
where X and Z are the matrices, b is the coefficient vector of fixed
effects, u is the coefficient vector of random effects, and e is an error
vector.
I would like
Thanks to Ben Bolker, I have put up the following page with tips and
linux/unix tools for creating html reports from Sweave:
http://biostat.mc.vanderbilt.edu/SweaveConvert
Reports can be opened in Word, or you can copy and paste from a browser
into OpenOffice or Word, keeping the formatting.
At first I would like to express the gratitude toward Oosting, J. and James
MacDonal for their help solving the previous problem I posted on this forum.
Now I have a new problem. I was trying to normalize set of microarray data
on HGU133 paltform. I put all the HGU133A platform chips in one
Patrick Drechsler wrote on 11 Jul 2006 02:10:21 MET:
[...]
I am now confronted with the problem that I have data which
requires a modelII regression (also called reduced major axes
regression (RMA) or geometric mean regression). For this I use
the function modelII (see below).
What
Hello!
Is there any possibility to use command completion in R-WinEdt?
Thanks
Franco
__
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
and
Nantachai
It seems as though you have created the model matrices from the matrix
representation of the model. This is unecessary as lme will construct those for
you from your data frame.
The first thing I recommend you do is use lmer instead of lme. Second, you
should look in the vignette in
Hi,
I was very much hoping someone could help me with the following.
I am trying to convert some SAS NLMIXED code to NLME in R (v.2.1),
but I get an error message. Does anyone have any suggestions?
I think my error is with the random effect u which seems to be
parametrized differently in the SAS
Gabor,
Your suggestion:
library(lattice)
xyplot(Consumption ~ Quarter, group = Year, data, type = o)
is very elegant indeed.
Thanks,
Jack.
Gabor Grothendieck [EMAIL PROTECTED] wrote: And if lattice is ok then try
this:
library(lattice)
xyplot(Consumption ~ Quarter, group = Year, data, type
Dear all,
I am currently analyzing a set of arrays hybe on the lattest affy Drosophila
2.0 GeneChip. I am trying to run simple annaffy analysis but can¹t find what
is the name of the annotation file I need to use.
Here is the output of the AffyBatch object I am using:
expData
AffyBatch object
Hi Marco,
I'm pretty sure you want to resend your question to the Bioconductor
mailing list. AFAICT your question is rather BioC specific and you
will likely get a more helpful response there.
+ seth
__
R-help@stat.math.ethz.ch mailing list
Mahdi Osman [EMAIL PROTECTED] writes:
Hi list,
I am interested in cluster analysis of microarray data. The data was
generated using cDNA method and a loop design.
I was wondering if any one has a suggestion about which package I can
use to analyse such data.
There are many packages
Here is yet another solution. This one consists only of
two gsubs and a function to reverse a string. It runs
at about the same speed as f3 but its main advantage
is how compact it is.
pat could be the same as before however we have made
use of Greg's discussion to use \\w\\w to
avail ourself
I need sequences that have gaps in them, such as the
following:
4 5 6 | 12 13 14 | 20 21 22
a simple question, I've been scratching my head for a
R function that will do this
The : and seq do not allow this, and the c() can be
used although not in an automatic way. I'm sure there
is a way to
Here are a few possibilies:
x - c(4, 12, 20)
rep(x, each = 3) + 0:2
rep(x, each = 3) + sequence(rep(3, length(x))) - 1
c(sapply(x, seq, length = 3))
On 7/25/06, etienne [EMAIL PROTECTED] wrote:
I need sequences that have gaps in them, such as the
following:
4 5 6 | 12 13 14 | 20 21 22
Hello sir:
How can I get S curve function via R?
For SPSS,the function is:y=exp(b0+b1/x)
Thanks a lot!
--
***
Xin Meng
Capitalbio Corporation
National Engineering Research Center
for Beijing Biochip Technology
50 matches
Mail list logo
