Hi everyone,
-
Coefficients:
ar1 ar2 ma1 ma2sar1 intercept drift
1.5283 -0.7189 -1.9971 0. 0.3982 0.0288 -9e-04
s.e. 0.0869 0.0835 0.0627 0.0627 0.1305NaN NaN
sigma^2 estimated as 0.04383: log
Dear list,
I would like to make a Sweave document that, for a given package,
automatically inserts the graphical outputs obtained when runing
the example code of high level plot functions. That is, for the
graphics base package, something like this:
http://pbil.univ-lyon1.fr/R/fichestd/tdr79.pdf
I've rolled up R-2.3.1.tar.gz a short while ago. This version contains
several minor fixups and removes a couple of bad bugs. See the full list
of changes below.
You can get it (in a short while) from
http://cran.r-project.org/src/base/R-2/R-2.3.1.tar.gz
or wait for it to be mirrored at a CRAN
On Thu, 1 Jun 2006, Jean lobry wrote:
Dear list,
I would like to make a Sweave document that, for a given package,
automatically inserts the graphical outputs obtained when runing
the example code of high level plot functions. That is, for the
graphics base package, something like this:
Good day R-Users,
I have a problem accessing some values in the output from the summary of an
lme fit.
The structure of my data is as shown below (I have attached a copy of the full
data).
id trials endp Z.sas ST
1 1 -1 -142.42884
1 1
Hi,
I'm a group of functions and I would like to create a package for load in R.
I have created a directory named INE and a directory below that named R, for
the files of R functions. A have created the files DESCRIPTION and INDEX in
the INE directory. The installation from local zip files, in
The minimum is to create a DESCRIPTION file, plus R and man
directories containing R code and .Rd files respectively.
It might help to run Rcmd CHECK mypkg before installation
and fix any problems it finds.
Googling for creating R package will locate some tutorials.
On 6/1/06, Rita Sousa
?package.skeleton
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gabor
Grothendieck
Sent: 01 June 2006 12:20
To: Rita Sousa
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] FW: How to create a new package?
The minimum is to create a DESCRIPTION file,
Pryseley:
You asked this question last week and received answers to that question
along with other suggestions from Doug Bates and myself. Is there some
reason the post by Andrew Robinson is not working?
Harold
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On
michael watson (IAH-C) wrote:
?package.skeleton
Folks, please!
Rita Sousa told you she already has a DESCRIPTION file.
Obviously, Rita forgot to build and *install* the package using
R CMD build
and
R CMD INSTALL
(the Meta directory is creating during the installation procedure)
Please
In that case, I have found the following useful:
http://www.biostat.jhsph.edu/~kbroman/Rintro/Rwinpack.html
-Original Message-
From: Uwe Ligges [mailto:[EMAIL PROTECTED]
Sent: 01 June 2006 12:43
To: michael watson (IAH-C)
Cc: Gabor Grothendieck; Rita Sousa; r-help@stat.math.ethz.ch
Dear Rita,
Do you want a package just for yourself, or something useful for others, with
docs, etc? I think the rest of the answers in this thread will help you
create a full fledged package. See also the detailed explanation
in Writing R extensions.
If you just want something quick and dirty
michael watson (IAH-C) wrote:
In that case, I have found the following useful:
http://www.biostat.jhsph.edu/~kbroman/Rintro/Rwinpack.html
This page is really outdated and not comprehensive, since some relevant
information is missing.
The manual contains even more *relevant* information,
hi list,
This question must be very basic but I am just 3 days
old to R. I am trying to find a function to merge two
tables of data in two different files as one.
Does merge function only fills in colums between two
table where data is missing or is there a way that
merge can be used to merge
hi list,
This question must be very basic but I am just 3 days
old to R, so I think i can ask. I am trying to find a
function to merge two
tables of data in two different files as one.
Does merge function only fills in colums between two
table where data is missing or is there a way that
merge
Ahamarshan jn wrote:
hi list,
This question must be very basic but I am just 3 days
old to R, so I think i can ask. I am trying to find a
function to merge two
tables of data in two different files as one.
Does merge function only fills in colums between two
table where data is missing
Dear R-users,
Could any give me a guidence to solve the following problem?
I have three balls A, B and C. I can pick A, B or C at the first
selection. At the second and third selections, I can either (i) pick
one of the remaining balls, or (ii) stop to pick any balls. Therefore,
I have a
Good day,
Yeah i tried most of the suggestions and still had the same problem. Some
people also suggested that i supplied a data set so they could verify what i am
saying. That is why i am still asking the same question with more details as
required.
My main focus now is to
Hi
I am struggling to split up the plots of the grouped objects in nlme in a
usable way. The standard plot command generates plots for each group on a
single page. When there are many groups however this does not look so good.
I have discovered the layout option which allows one to split up
OK, in the future it might be helpful to note what you tried and what
error messages you received when you worked with the syntax provided.
Why did VarCorr() not work?
In general, there seems to be a larger probelm with your model as Doug
Bates noted. I would also suggest that you go about this
Hi
I have found some postings referring to the fact that one can try and
understand why a particular model is failing to solve/converge from the
verbose output one can generate when fitting a nonlinear mixed model. I am
trying to understand this output and have not been able to find out much:
Could you tell us what version (and variant) of Linux you're using, and the
GCC version (see gcc -v)? Are you logged in as root when you built R, or as
regular user? What does the summary output of configure look like (i.e.,
the last 25 lines or so of output from running configure)?
Andy
On Thu, 2006-06-01 at 09:16 -0400, Chuck Cleland wrote:
Ahamarshan jn wrote:
hi list,
This question must be very basic but I am just 3 days
old to R, so I think i can ask. I am trying to find a
function to merge two
tables of data in two different files as one.
Does merge
I want to get index number of 0s, like
x - c(1, 2, 0, 3, 0, 4)
I want an output of
3, 5
Thank you.
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
I am having difficulty understanding the syntax of the nlme() function for
nonlinear mixed models. The data frame is called Marouane.chlorophyll. The
model involves a dependent variable (Absorb) and an independent variable
(Ch.surf), which are both numeric variables in the data frame. The data
On Thu, 2006-06-01 at 11:06 -0400, Mu Tian wrote:
I want to get index number of 0s, like
x - c(1, 2, 0, 3, 0, 4)
I want an output of
3, 5
Thank you.
See ?which
which(x == 0)
[1] 3 5
HTH,
Marc Schwartz
__
R-help@stat.math.ethz.ch
try:
which(x == 0)
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
Hello;
I'm preparing a simple function that plots sequentially a number of time
series for visual inspection, after some previous transformations.
To allow the user to have time to see the series I'm using
Sys.sleep(1.5) between calls to line(). But is it possible to add some
mechanism that
I am using the following command from the Fcalendar Package :
x = timeSequence(1992-12-31,1994-12-31)
and then y = as.character(x) is a vector of character strings
[1] 1992-12-31 1993-01-31 1993-03-03 1993-03-31 1993-05-01 etc
This is very close to what I need and thank you very much
to
Under the below section it talks about X86_64 architecture for Linux.
1. tar -zxvf R-2.3.0.tar.gz
2. changed the directory to the newly created directory R-2.3.0
3. Typed ./configure 'r_arch=name'
4. Typed make
5. Make check
6. make check-all
7. Typed make install
8. Typed R
I still get the same
Michael Dondrup [EMAIL PROTECTED] writes:
Hi Karin,
I would like to help with this, but it's not completely clear to me
what your vioplots look or should look like. Could you post a little
reproducible code example to the list, and then specify what should be
different?
OK.
Adapted from
Dear R-users,
I have longitudinal data and would like to fit a model where both the
variance-covariance matrix of the random effects and the residual variance are
conditional on a (binary) grouping variable.
I guess the model would have the following form (in hierarchical notation)
Yi|bi,k ~
I'm an R novice, so I hope my question is a valid one. I'm trying to
run the following script in the current version of R.
for (i in 1:1640){for (j in (i+1):1641){
if (i == 1 j == 2){x - cor(sage[i,],sage[j,],method=spearman); y
- cor(frie[i,],frie[j,],method=spearman)}
if (i != 1 || j !=
You should try to figure out how large x and y need to be, and allocate them
before the computation. Growing objects like that in R is very, very,
very bad, as you've found out.
Andy
From: Tim Alcon
I'm an R novice, so I hope my question is a valid one. I'm trying to
run the following
Compare
system.time({
v - vector()
for (i in 1:10^5) v - c(v,1)
})
with
system.time({
v - vector(length=10^5)
for (i in 1:10^5) v[i] - 1
})
If you don't know exactly how long v will be, use a value that's large
enough, then throw away what's extra.
-Original Message-
Hi Karin,
I think it's just setting horizontal to FALSE.
vioplot(normal,uniform,horizontal=FALSE)
For the axis:
How about this:
vioplot(normal, uniform)
axis(2, pos=1.5)
sometimes the second axis might overlap, so try this one too:
ylim - range(normal, uniform)
par(mfrow=c(1,2), bty='n')
I have the following problem which seems really strange but it's more likely
something that I don't understand ?
The code below works in that it adds 1 day to all of the dates.
( NB : the Fcalendar package needs to be loaded for this code to run )
startDate=1992-12-31
endDate=1994-12-31
Hi All.
(This is probably an incredibly stupid question :))
I have just noticed that every time I start an R2.3 session and request
a random number (say with rnorm) that I receive the same number. AFAIK I
have not set a seed at any point. Returning to an older version of
R(v1.8) I get the
Hi,
I'm creating some lattice plots that have a key and I'd like to put a
label on the key. The problem is that the text label for the key
prevents the values of the group variables from being shown (see example
below). I don't think that this is a feature, but I might be abusing the
title arg
On 6/1/2006 3:13 PM, Carl wrote:
Hi All.
(This is probably an incredibly stupid question :))
I have just noticed that every time I start an R2.3 session and request
a random number (say with rnorm) that I receive the same number. AFAIK I
have not set a seed at any point. Returning to an
according to ?.Random.seed, `Note' section, different R sessions give
differen simulation results; maybe each time you start R you load a
previously saved workspace, with an existing .Random.seed, i.e., check
if the .Random.seed value is the same each time you start R.
I hope it helps.
Best,
Kuhn, Max wrote:
Hi,
I'm creating some lattice plots that have a key and I'd like to put a
label on the key. The problem is that the text label for the key
prevents the values of the group variables from being shown (see example
below). I don't think that this is a feature, but I might be
That certainly sounds like the behaviour you would get if you
had a .Random.seed in your work space. If you do not save your
workspace at the end of the session then the random seed will
be in exactly the same state every time you start a new R session,
and you will get identical simulations from
Hi Max,
You have 'title' as a component of the 'text' component of key.
According to the help page for xyplot, 'title' is a component
of key. AFAICS, 'text' does not have a 'title' component.
Peter Ehlers
Kuhn, Max wrote:
Hi,
I'm creating some lattice plots that have a key and I'd like to
Quicker and dirtier is to simply save a workspace with your desired
functions and attach() it automatically (e.g. via .First or otherwise) at
startup (?Startup for details and options). Of course none of the tools and
benefits of package management are available, but then I did say quicker and
Hi All.
(Just to stem the growing tide of helpful replies pointing out my
stupidity).
It as as everyone suspects - I have at some point saved my workspace so
there is a .Random.seed lurking. This is something I don't normally do,
hence haven't encountered this before... it also explains why
Peter and Sundar,
Thanks, I did miss that. I used:
key = list(
title = More Text,
cex.title = 1,
columns = 2,
text=list(
to get what I wanted.
Max
-Original Message-
From: P Ehlers [mailto:[EMAIL PROTECTED]
Sent: Thursday, June 01, 2006 4:14 PM
To: Kuhn, Max
Cc:
Greetings,
I am having a very hard time with a nonlinear regression. The last chance is
that maybe somebody can spot something wrong
The data and the model are
described below:
number of observations = 3030
y = [0,
,~16]
D1969 = [.16,
,~70,000]
mod=nls(log(D1969)~d-log(1+d1*exp(-gt+g1*y)),
Hello all:
I have searched through the help files and I have been unsuccessful in
solving this problem.
I am trying to create a small wrapper function that will go around a
call to a plot function and create a windows metafile in the directory
and also write the name of the file to a text
Dear List:
I have the follow code:
y - replicate(10,replicate(8,sum(rnorm(rpois(1,5)
Now I need to apply the following condition to _every_ randomly generated
Normal number in the code above:
x - max(0,x-15) + max(0,x-90), where x represents the individual Normal
numbers.
In other words,
A very simple one but I cannot figure it out. I have a matrix and I want to
display it in a seprate window. I would really like it if it worked on both
linux and windows.
cheers
[[alternative HTML version deleted]]
__
What can I do?
Thanks a lot!
[[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! http://www.R-project.org/posting-guide.html
Hi
This is a second announcement and call for abstracts.
Please forward and circulate to other interested parties.
DSC 2007, a conference on systems and environments for statistical
computing, will take place in Auckland, New Zealand on February 15 16,
2007.
We invite abstracts on the
Hi.
My data contains information for 10 hospitals for 12 different measures.
Let's call it x1-x12. I need to create a boxplot for each one of this
measures and put them into one page. Each plot also needs to be independent,
i.e. cannot use the group feature because of different scales for each
Michael == Michael [EMAIL PROTECTED] writes:
Michael What can I do?
Implement this feature and send it to the person responsible for the
code for inclusion?
Michael Thanks a lot!
Indeed, you contribution would be very much appreciated.
Cheers,
Berwin
55 matches
Mail list logo