Hi
Im attempting to build a sample variogram for 300 obersvations
of longitudinal data. So what I need to do is compute the half
squared differences between pairs of residuals (for instance
if a subject has 4 obersvations, this is 4 choose 2 paird differences)
for each subject.
Also, then I
1.The simple problem ist that i have different Samples , from
which i make repeated measurements (each sample is measured 6
times) and i repeat this experiment over several Days, so i
get the lme grouping term random=~1|Days/Sample.
I would rather specify this as: random= ~ 1 |
A variable is considered to be nested when it is unique within that
combination. i.e. if your measuring Apparatus was different at each sample.
Hence, as Lorenz suggested it is reasonable to treat Apparatus as fixed,
specially if you are considering making comparisons between the 2 apparatus.
JohnF == John Fox [EMAIL PROTECTED]
on Sun, 24 Oct 2004 12:37:09 -0400 writes:
JohnF Dear Sarah, If the data are allegedly bivariate
JohnF normal, then they are probably two vectors, not
JohnF one. Assuming that this is the case, I know of
JohnF nothing quite as neat as a
Actually, you are unable to connect to the Internet from R on your
computer. Perhaps you are not connected, or need to set a proxy?
Please read the rw-FAQ, as the posting guide asks.
Note that right now you there are no packages on CRAN and Bioconductor has
not yet released packages for R
Hi Everybody:
I want to draw some chart using command histogram and add another curve in it.
Example:
require(stats)
data(singer)
library(lattice)
histogram( ~ height | voice.part, data = singer, nint = 17, endpoints = c(59.5,
76.5), layout = c(2,4), aspect = .5, xlab = Height (inches))
Now
Hi YiYao,
you need the `?panel.abline()' function, somehing like:
panel=function(x, breaks, equal.widths, type, ...){
panel.histogram(x, breaks, equal.widths, type, ...)
panel.abline(v=mean(singer$height))
}
inside the `histogram()' function could do the trick.
I hope
Hi Dimitris:
Thanks for your help, I will try.
BR
Yiyao
--
: Dimitris Rizopoulos [mailto:[EMAIL PROTECTED]
be]
: 20041025 15:39
: YiYao_Jiang
: [EMAIL PROTECTED]
: Re: [R] Multiple formula in one block
Hi YiYao,
you need the `?panel.abline()' function, somehing like:
panel=function
Johannes == Johannes Graumann [EMAIL PROTECTED]
on Sat, 23 Oct 2004 11:04:25 -0700 writes:
Johannes == Johannes Graumann [EMAIL PROTECTED]
on Sat, 23 Oct 2004 11:04:25 -0700 writes:
Johannes Hello,
Johannes I seem unable to construct a legend which contains
Johannes a
I write to ask you an help about the package gregmisc.
I saw the instructions, and I need some functionalities of this package,
but I am unable ti download it.
On friday I was able to download thte .zip, but R does not install this
package, today there is no possibility to download it.
What I
Gregmisc is a *bundle*, not a library, not a package. It unpacks into
four packages: see
http://cran.r-project.org/src/contrib/Descriptions/gregmisc.html
I guess you are using Windows (see the posting guide, and tell us).
There is a temporary problem on the CRAN master. Set
Prof Brian Ripley wrote:
Actually, you are unable to connect to the Internet from R on your
computer. Perhaps you are not connected, or need to set a proxy?
Please read the rw-FAQ, as the posting guide asks.
Note that right now you there are no packages on CRAN and Bioconductor has
not yet
On Mon, 25 Oct 2004 10:10:06 +0200,
Uwe Ligges (UL) wrote:
Prof Brian Ripley wrote:
Actually, you are unable to connect to the Internet from R on your
computer. Perhaps you are not connected, or need to set a proxy?
Please read the rw-FAQ, as the posting guide asks.
Thanks for your reply.
Considering the reproducibility of R this is an important issue.
Regarding my particular problem, i.e. the package build
failure with R CMD check package in the examples section,
I do not know how to proceed. I have tried the following
strategy: direct installation of the
Hi,
On Mon, 25 Oct 2004, Anna Maria Paganoni wrote:
I write to ask you an help about the package gregmisc.
I saw the instructions, and I need some functionalities of this package,
but I am unable ti download it.
On friday I was able to download thte .zip, but R does not install this
Hello,
am I using 'setIs' in the correct way in the subsequent (artifical) example?
Do I have to specify explicit 'setAs' for 'list' and 'vector' or
should this work automatically, since getClass(List1) states
an explicit coerce also for these classes.
I'm working with R 2.0.0 Patched
Hi ALL,
I am trying to compute adjusted SS estimated component variance
in GLM with un-balanced data using R. Can anyone advise me? Thanks in
advance.
Best Regards,
WeiQiang Li
[[alternative HTML version deleted]]
__
[EMAIL
Hi Matthias,
A similar problem to yours (with one level of inheritance less) was
disccussed this month on the r-devel list.
You find an answer from JChambers here:
https://stat.ethz.ch/pipermail/r-devel/2004-October/030980.html
And yes specifying _setAs_ to each _setIs_ with the coerce and
Thank you,
Matthias
Hi Matthias,
A similar problem to yours (with one level of inheritance less) was
disccussed this month on the r-devel list.
You find an answer from JChambers here:
https://stat.ethz.ch/pipermail/r-devel/2004-October/030980.html
And yes specifying _setAs_ to each
Martin Maechler wrote:
Johannes == Johannes Graumann [EMAIL PROTECTED]
on Sat, 23 Oct 2004 11:04:25 -0700 writes:
Johannes == Johannes Graumann [EMAIL PROTECTED]
on Sat, 23 Oct 2004 11:04:25 -0700 writes:
Johannes Hello,
Johannes I seem unable to construct a legend which contains
Hi Matt,
there are several R packages that will compute the sample variogram for you.
Check out GeoR, sgeostat, nlme, spatial. There's no point in recoding the
whole lot yourself, unless as a learning excercise.
D
p.s. For time series autocorrelations, you could use acf in package stats.
Atte Tenkanen wrote:
Is it possible to input data from sound card of the computer to R? What
do I need to know about my computer (it's linux pc)? Can I get some real
time graphical information this way, spektrum for example?
I guess it is possible to get data directly from the sound card using
hi,
i need a function that can replace na's in univariate ts objects with a value
interpolated linearly
using the adjacent non-na values. for example, the function i need (myfun) would do
the following
x - ts(10,11.4,NA,12,9.7)
y - myfun(x)
y
10 11.4 11.7 12 9.7
i can code an na.action
Dear all
I was looking at the aov documentation page and came across the
following which seems like a contradiction to me:
This provides a wrapper to |lm| for fitting linear models to balanced
or unbalanced experimental designs. (I presume 'This' refers to aov)
and
|aov| is designed for
On Mon, 25 Oct 2004, AJ Smit wrote:
I was looking at the aov documentation page and came across the
following which seems like a contradiction to me:
This provides a wrapper to |lm| for fitting linear models to balanced
or unbalanced experimental designs. (I presume 'This' refers to aov)
AJ Smit [EMAIL PROTECTED] writes:
Dear all
I was looking at the aov documentation page and came across the
following which seems like a contradiction to me:
This provides a wrapper to |lm| for fitting linear models to
balanced or unbalanced experimental designs. (I presume 'This' refers
Since installing R 2.0 I've had the following 2 issues
I can't find solutions for. I'm using Windows XP.
1) R 1.9 continues to read my .Rprofile file in my
home directory but R 2.0 does not. How do I get R 2.0
to do this?
2) The following function no longer works in R 2.0 but
continues to work
how about
hist(rep(x,w))
i think it works.
I have the numeric vector `x' and the weights `w' for it of the same
length. I'd like to have hist of `x' weighted with `w'.
How to do that? I found nothing about weights in the hist help page..
__
[EMAIL
I am about to write general functions to read the output of simulations
models.
These model generate output files with different sections which I want
to analyze plot etc.
Since this will be used many people at the department I wanted to make
sure that will do this in the best way.
For
Hi All:
A couple of weeks back, I asked a question on the list that I was invited to provide
an introductory lecture on R to a group of academicians in Kolkata. I thank all of you
who had generously guided me in providing me web links and words to the wise.
Time to give back. I did the
Hello,
I'm new to the installation R on Linux and I've followed the
instructions on the R Installation and Administration Manual and on the
FAQs. I also have changed the path at the .bash_profile file but when I
try to type R at the shell prompt, it says:
Fatal error: R home directory is not
Henrik Andersson h.andersson at nioo.knaw.nl writes:
:
: I am about to write general functions to read the output of simulations
: models.
:
: These model generate output files with different sections which I want
: to analyze plot etc.
:
: Since this will be used many people at the
I presume R when executed from its bin directory works fine.
Are you actually using bash shell. Try 'echo $SHELL'.
If yes, then try 'which R' and see if it pointing to the right path.
On Mon, 2004-10-25 at 15:55, Alicia Amadoz wrote:
Hello,
I'm new to the installation R on Linux and I've
On Mon, 2004-10-25 at 16:04, Peter Dalgaard wrote:
Alicia Amadoz [EMAIL PROTECTED] writes:
What happens if you type 'which R'? What is the version of R and what
is in your path?
I suspect you are somehow trying to run the R binary, bypassing the R
shell script, which sets things like
Please excuse my previous email. I was going to follow up Peter's post
showing that if I run the R binary directly I get exactly the error
message reported by the original poster. Unfortunately, after I decided
that was too trivial to post, as my mouse headed towards the X button to
kill the email
Suzette Blanchard [EMAIL PROTECTED] writes:
I am not able to download packages (xtable) from CRAN
not even just a .zip file, I get that the page can not be
found. Please could you help?
Suzette
Which CRAN mirror? There was a messup over the weekend where a router
in Dortmund went
Suzette Blanchard wrote:
I am not able to download packages (xtable) from CRAN
not even just a .zip file, I get that the page can not be
found. Please could you help?
Suzette
=
Suzette Blanchard, Ph.D.
Research Scientist
Frontier Science
Hi Suzette,
Hi Suzette,
did you try several CRAN mirrors, for instance I'm able to download it
from a UK mirror:
http://www.stats.bris.ac.uk/R/src/contrib/Descriptions/xtable.html
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic
Same behavior with Firefox 1.0Pr.
/E
Thomas Schönhoff wrote:
Hello Arin,
Arin Basu schrieb:
A couple of weeks back, I asked a question on the list that I was
invited to provide an introductory lecture on R to a group of
academicians in Kolkata. I thank all of you who had generously guided
me in
Hi,
I ran into an odd problem with the print command for R-2.0 on a windows
machine.
The icelandic character thorn (þ,Þ) which is included in in the Latin-1
character set [iso/8859-1 char# 222 (upper case) and #254 (lower case)]
prints out incorrectly. Instead of getting the correct
character I
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
[1]
Hi,
I have a package of my own which seems to work fine under linux.
I want to make a compiled version for windows.
(I work with windows 2000 and R 2.0)
I have followed the steps described in the file
readme.packages (in the top-level directory of the binary installation)
and I had a trouble at
On Mon, 25 Oct 2004, Halldor Björnsson wrote:
Hi,
I ran into an odd problem with the print command for R-2.0 on a windows
machine.
The icelandic character thorn (þ,Þ) which is included in in the Latin-1
character set [iso/8859-1 char# 222 (upper case) and #254 (lower case)]
prints out
Thank you so much ... works now ... sooo much to learn ...
Joh
On Mon, 25 Oct 2004 09:56:35 +0200
Martin Maechler [EMAIL PROTECTED] wrote:
Johannes == Johannes Graumann [EMAIL PROTECTED]
on Sat, 23 Oct 2004 11:04:25 -0700 writes:
Johannes == Johannes Graumann [EMAIL PROTECTED]
Halldor Björnsson [EMAIL PROTECTED] writes:
Hi,
I ran into an odd problem with the print command for R-2.0 on a
windows machine.
The icelandic character thorn (þ,Þ) which is included in in the Latin-1
character set [iso/8859-1 char# 222 (upper case) and #254 (lower case)]
prints out
Dear R users,
I need to fit an ARMA model. As far as I've seen, EACF (extended ACF)
is not available in R.
1. Let's say I fit a series of ARMA models in a loop. Given the
code/output included below, how do I pull 'Model' and 'Fit' (AIC)
from each summary() so that I can combine them into an
Okay thanks!
No I have the following difficulty implementing a variogram in the
nlme package:
cd1 - lme(count ~ time, data=cd4,random= ~ time | id)
plot(Variogram(cd1, form= ~ time | id, robust=TRUE))
Error in as.array(X) : attempt to set an attribute on NULL
Im not sure how to fix this, and
[EMAIL PROTECTED] writes:
Hi,
does anybody know why the following is not working:
hist(rnorm(200))
box(bty=o)
gives me a box without rounded corners.
Because that isn't what it is supposed to do. Did you expect 7 to
give you a slanting right edge? And which letter should represent a
On Mon, 25 Oct 2004 15:56:37 +0200, Henrik Andersson
[EMAIL PROTECTED] wrote :
I am about to write general functions to read the output of simulations
models.
These model generate output files with different sections which I want
to analyze plot etc.
Since this will be used many people at the
On Mon, 25 Oct 2004, Gilles GUILLOT wrote:
Hi,
I have a package of my own which seems to work fine under linux.
I want to make a compiled version for windows.
(I work with windows 2000 and R 2.0)
There is no such version of R.
I have followed the steps described in the file
On Mon, 2004-10-25 at 12:02, [EMAIL PROTECTED] wrote:
Hi,
does anybody know why the following is not working:
hist(rnorm(200))
box(bty=o)
gives me a box without rounded corners.
You will not get rounded corners with box().
The 'bty' argument determines how many sides of the plot
Gilles GUILLOT wrote:
Hi,
I have a package of my own which seems to work fine under linux.
I want to make a compiled version for windows.
(I work with windows 2000 and R 2.0)
I have followed the steps described in the file
readme.packages (in the top-level directory of the binary installation)
Date: Mon, 25 Oct 2004 14:10:15 -0400
MIME-Version: 1.0
X-Mailer: Internet Mail Service (5.5.2654.89)
Content-Type: text/plain; charset=iso-8859-1
It is possible to save all of these items using the save.session and
restore.session functions provided by the 'session' package available from
CRAN.
Hello,
I'm blotting a series of growth curves into a multiplot environment
created with layout().
since I want the four plots to be easily visually comparable, I do the
following:
#first plot
plot(x,y,stuff)
standarduser-par()$usr
...
some fitting
...
lines(spline(x, fitted_equation))
Dear R- helpers
Following a draft structure of the R script for which I am facing problem
Step 1
x - of type array with original values
y - of type array with original values
Step 2
for (ctr in 1:10) {
# my problem here the both x and y still show the original values from step 1
# in
On Mon, 25 Oct 2004, bogdan romocea wrote:
Dear R users,
I need to fit an ARMA model. As far as I've seen, EACF (extended ACF)
is not available in R.
But PACF is as well as ACF.
1. Let's say I fit a series of ARMA models in a loop. Given the
code/output included below, how do I pull
Just a little addendum to Martin's comments below. It is well known that
using LS centers and covariances for the M-distances is generally not a good
way to do this, as these statistics, themselves, are distorted by the long
tails (do 1D distributions have tails?) so that the problems are
Hello,
I want to do feature selection for classification purpose (using lda). Can someone
point me to any R package or S-plus package for this? Something like SFS or SFFS
method would be useful for me.
Thanks.
Raj
[[alternative HTML version deleted]]
Johannes Graumann wrote:
Hello,
I'm blotting a series of growth curves into a multiplot environment
created with layout().
since I want the four plots to be easily visually comparable, I do the
following:
#first plot
plot(x,y,stuff)
standarduser-par()$usr
...
some fitting
...
lines(spline(x,
Kunal Shetty wrote:
Dear R- helpers
Following a draft structure of the R script for which I am facing problem
Step 1
x - of type array with original values
y - of type array with original values
Step 2
for (ctr in 1:10) {
# my problem here the both x and y still show the original values from
On Mon, 2004-10-25 at 13:19, Johannes Graumann wrote:
Hello,
I'm blotting a series of growth curves into a multiplot environment
created with layout().
since I want the four plots to be easily visually comparable, I do the
following:
#first plot
plot(x,y,stuff)
standarduser-par()$usr
You were right! 'ylim' does what I want!
Thanks for deciphering my cryptic scribble and helping out!
Joh
On Mon, 25 Oct 2004 14:25:55 -0500
Marc Schwartz [EMAIL PROTECTED] wrote:
On Mon, 2004-10-25 at 13:19, Johannes Graumann wrote:
Hello,
I'm blotting a series of growth curves into a
Uwe Ligges
thank you for u prompt reply
my problem was in step3
where my function returns two different arrays.
Yes i did try returning an object from the array.
but the problem became...i tired returning the two arrays in a dataframe such as
newXY -
Robert Sams wrote:
hi,
i need a function that can replace na's in univariate ts objects with a value interpolated linearly
using the adjacent non-na values. for example, the function i need (myfun) would do the following
x - ts(10,11.4,NA,12,9.7)
y - myfun(x)
y
10 11.4 11.7 12 9.7
i can code
Sundar Dorai-Raj wrote:
Robert Sams wrote:
hi,
i need a function that can replace na's in univariate ts objects with
a value interpolated linearly using the adjacent non-na values. for
example, the function i need (myfun) would do the following
x - ts(10,11.4,NA,12,9.7)
y - myfun(x)
y
10
Hi everyone,
My purpose is to read a .CEL file into R.
The .CEL file was created from a .CAB by using DTT software found on
Affymetrix website
I read the .CEL file in R using ReadAffy as follows:
d2=ReadAffy(widget=T)
and I complete the fields as required.
It does not complain. For example I
Much improved (patched Firefox 1.0pre on RH9)
On 25 Oct 2004, Arin Basu wrote:
Hi All:
This follows my earlier post on webized slides on lecture presentation on
introducing R. I learned that in Mozilla (Firefox) browsers, the slides did not show
up. Sorry for the no show. As a reluctant
A colleague of mine is trying to use nls() to effect an optimization,
and is encountering a scoping problem. I should know how to solve it
for him but well, I just don't.
I also had a quick scrounge of the archives --- I know I've seen this
topic addressed before --- but I couldn't track
This is a question that should be asked on the bioconductor mailing list
rather than R-help.
Quick answer:
library(reposTools)
install.packages2(whateverchiptypethisiscdf)
Where you have to replace whateverchiptypethisiscdf with e.g.,
hgu133acdf.
HTH,
Jim
James W. MacDonald
Affymetrix and
Dear Bert,
The data.ellipse() function in the car package optionally uses the
covariance matrix and location vector returned by cov.trob() (from the MASS
package). I believe that any 1D function of the two variables is potentially
problematic. As to why do it -- comparing the bivariate
Raj,
look at academic.uprm.edu/eacuna/softw.htm, I built a library of R
functions for data preprocessing tasks including feature selection
for supervised classification.
Please send me your comments.
Edgar
On Mon, 25 Oct 2004, Rajdeep Das wrote:
Hello,
I want to do feature selection for
From: Berton Gunter
Just a little addendum to Martin's comments below. It is well
known that
using LS centers and covariances for the M-distances is
generally not a good
way to do this, as these statistics, themselves, are
distorted by the long
tails (do 1D distributions have
Hi,
I am looking for a means to compute empirical distribution function for
a very large data set and evolution of that edf with time.
Here are some specifics. Each day I have an estimate of a distribution
function and a new sample of about 1e4 points from the distribution in
question. I want
Rolf Turner rolf at math.unb.ca writes:
:
: A colleague of mine is trying to use nls() to effect an optimization,
: and is encountering a scoping problem. I should know how to solve it
: for him but well, I just don't.
:
: I also had a quick scrounge of the archives --- I know I've seen
Kunal Shetty kshe4 at student.monash.edu writes:
:
: Uwe Ligges
:thank you for u prompt reply
: my problem was in step3
: where my function returns two different arrays.
: Yes i did try returning an object from the array.
: but the problem became...i tired returning
75 matches
Mail list logo