Hi!
Suspicion Breeds Confidence!
--Brazil
Take a look at the packages listed on cran.r-project.org
or on the r prompt (if all R-packages installed type)
help.search(differential equations)
Help files with alias or concept or title matching 'differential
equations' using fuzzy
On Mon, 7 Jun 2004, Laura Holt wrote:
Dear R People:
Is it possible to fit an AR model such as:
y_t = fee_1 y_t-1 + fee_2 y_t-9 + a_t,
please?
I know that we can fit an AR(9) model, but I was wondering if we could do a
partial as described.
Yes. See ?arima, and especially the
Rui wrote:
Hi folks,
I have a question about how to load a dll.
First, I use the command
dyn.load(lassofu.dll);
then, I got the message below
NULL
Warning message:
DLL attempted to change FPU control word from 9001f to 90003;
See ?dyn.load.
This message should make you a little bit worried
Sorry, it's simple question:
example :
data.test
label value
1 one21.35746
2 two22.07592
3 three20.74098
I would like the return :
label value
3 three20.74098
1 one21.35746
2 two22.07592
Anyone can help it?
Thanks,
try to use:
data.test[order(data.test[,2]),]
I hope this helps.
Best,
Dimitris
Dimitris Rizopoulos
Doctoral Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/396887
Fax: +32/16/337015
Web:
Spencer Graves [EMAIL PROTECTED] writes:
Data: DF
log-likelihood: -55.8861
Random effects:
Groups NameVariance Std.Dev.
smpl (Intercept) 1.7500e-12 1.3229e-06
Estimated scale (compare to 1) 3.280753
Fixed effects:
Estimate Std. Error z value Pr(|z|)
Hello,
Here is a slightly more sophisticate and fully vectorized, answer.
RecodeChoices - function(mat) {
# Make sure mat is a matrix (in case it is a data.frame)
mat - as.matrix(mat)
# Get dimensions of the matrix
Dim - dim(mat)
Nr - Dim[1]
Nc -
On Mon, 7 Jun 2004, Jingky Lozano wrote:
Good day R help list!!!
I've been trying to do Bootstrap in R on Censored data. I encountered
WARNING/ERROR messages which I could not find explanation.
I've been searching on the literature for two days now and still can't find
answers. I hope
Hi!!
I am a new user of R (just trying to analysis microarrays with some packages from the
bioconductor project).
I would like to import a text-delimeted file containing 20 columns and 22200 rows.
I have tried
read.table;
scan(file=)
matrix(scan(file, n=20*200,20,200, byrow=TRUE));
Doesn't
Please use a more appropriate subject!
- Original Message -
From: [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Tuesday, June 08, 2004 8:25 PM
Subject: [R] (no subject)
Hi!!
I am a new user of R (just trying to analysis microarrays with some
packages from the bioconductor project).
Hi!
The message means that R does not find the file.
Under windows (if you are using windows) you must escape the backslash with a
backslash \.
or give us an example how you provide the path to the file.
Or use:
file.choose()
or.
library(tkWidgets)
fileBrowser()
both functions will
Hi.
First, I think you should address questions about aroma directly to me (the
author) instead of the r-help list since it is not a core R package.
To answer you question, you can not make a GenePixData object out of a
RawData object because the latter contains much less information than a
Marc Schwartz wrote:
On Mon, 2004-06-07 at 15:51, Gavin Simpson wrote:
snip
Thanks Roger and Marc, for suggesting I use ./tools/rsync-recommended
from within the R-patched directory.
This seems to have done the trick as make completed without errors this
time round. The Recommended directory
Hi,
I have been trying to download BradleyTerry and brlr software but in
vain. Can you help me please. I already have R installed in my computer.
Thank you in advance.
Elias Obudho.
-
University of Nairobi Mail Services
You can't afford to stay
Here I got a function :
f-function(x){
res-list()
res$bool=(x=0)
res$tot=x
return(res)
}
I want that a=res$bool and b=res$tot
I know that a possible solution is :
Result=f(x)
a=Result$bool
b=Result$tot
But I don't want to use the variable Result.
I'd like to write directly something like
russell wrote:
Hello,
I'm wondering if anyone has used the msm package to compute the steady
state probabilities for a Markov model?
There's no built-in function in msm to do this, but this would be a
useful feature. For discrete time Markov chains this is a matter of
finding the eigenvector
On Tue, 2004-06-08 at 21:44, zze-PELAY Nicolas FTRD/DMR/BEL wrote:
Here I got a function :
f-function(x){
res-list()
res$bool=(x=0)
res$tot=x
return(res)
}
I want that a=res$bool and b=res$tot
I know that a possible solution is :
Like
foo - function(x,...){
[EMAIL PROTECTED] wrote:
Hi,
I have been trying to download BradleyTerry and brlr software but in
These are packages.
vain. Can you help me please. I already have R installed in my computer.
Thank you in advance.
Version of R?
Operating system and its version?
What was the error message after
Hi,
i would like recocde some numeric variables in one
step, but hanging unexpected in a level asignment problem?
for(i in 2:length(msegmente))
{ msegmente[,i] - as.factor(msegmente[,i])
}
Problem is that not every level is in every variable, so the
asignment is necessary!?
On Tue, 8 Jun 2004, Gavin Simpson wrote:
Marc Schwartz wrote:
On Mon, 2004-06-07 at 15:51, Gavin Simpson wrote:
snip
Thanks Roger and Marc, for suggesting I use ./tools/rsync-recommended
from within the R-patched directory.
This seems to have done the trick as make completed
On Tue, 8 Jun 2004 13:09:25 +0200 Christian Schulz wrote:
Hi,
i would like recocde some numeric variables in one
step, but hanging unexpected in a level asignment problem?
for(i in 2:length(msegmente))
{ msegmente[,i] - as.factor(msegmente[,i])
}
Problem is that not every level is
Hi,
Thanks for your message. I tried to prove that the
R-squared of the averaged model is always greater than
or equals to the average R-squared of individual
models (supposed m=2), Please see the attached r2.pdf.
I hope this can be generalized to general case (m
2).
Any comment would be very
Dear R-users;
I am interested in getting the number of objects in a list, and have thus
far been unsuccessful. Can you please help.
Thank you.
Vumani
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Vumani Dlamini wrote:
Dear R-users;
I am interested in getting the number of objects in a list, and have
thus far been unsuccessful. Can you please help.
Thank you.
Vumani
__
[EMAIL PROTECTED] mailing list
Lists in the S language are like vectors, so length(mylist) would tell you
how many components there are in mylist. Beware, though, that lists can
well be nested...
Andy
From: Vumani Dlamini
Dear R-users;
I am interested in getting the number of objects in a list,
and have thus
far
many thanks!
christian
I'm not sure what exactly you are trying to do, but does replacing
as.factor(msegmente[,i])
by
factor(msegmente[,i],
levels = c(1, 3, 10, 29, 32, 38, 52, 53, 59),
labels = c(AK,GC,OC,AM,IA,ACH,ZBA,A9L,EHK))
yield the desired result?
hth,
On Sun, 6 Jun 2004 14:13:59 -0400, Charles and Kimberly Maner
[EMAIL PROTECTED] wrote :
Hello. I have researched this topic and have found no answers. I am
running R 1.9.0 and am trying to print a lattice graph, (e.g., xyplot(1~1)),
using mouse right click - print. It produces a blank page.
Hello,
I suppose this is a basic question but couldn't find a solution.:
I have a large matrix with let say 3 columns:
V1 V2 V3
a x 2
a x 4
a y 8
b z 16
and I want to compute some statistics based on
the levels resulting form the
?tapply
?aggregate
You probably have a data frame, not a matrix.
Marc Mamin wrote:
I suppose this is a basic question but couldn't find a solution.:
I have a large matrix with let say 3 columns:
V1 V2 V3
a x 2
a x 4
a y 8
b z 16
and I
Marc Mamin wrote:
Hello,
I suppose this is a basic question but couldn't find a solution.:
I have a large matrix with let say 3 columns:
V1 V2 V3
a x 2
a x 4
a y 8
b z 16
and I want to compute some statistics based on
the levels resulting
On Tue, 2004-06-08 at 06:23, Prof Brian Ripley wrote:
On Tue, 8 Jun 2004, Gavin Simpson wrote:
Marc Schwartz wrote:
On Mon, 2004-06-07 at 15:51, Gavin Simpson wrote:
snip
snip
Perhaps I am being dense, but in reviewing the two documents (R Admin
and the CRAN sources
Spencer Graves [EMAIL PROTECTED] writes:
Another GLMM/glmm problem: I simulate rbinom(N, 100, pz),
where logit(pz) = rnorm(N). I'd like to estimate the mean
and standard deviation of logit(pz). I've tried GLMM{lme4},
glmmPQL{MASS}, and glmm{Jim Lindsey's
Dear all,
there has been a thread on missing Rd.sty, when converting Rd-files
to DVI-files, see
http://maths.newcastle.edu.au/~rking/R/help/03b/8179.html.
This problem could be solved by editing Rd2div.sh manually, but then
I got the following error message:
Hi, Peter:
Thanks. The help page on GLMM in lme4 0.6-1 2004/05/31 mentions
GLMM(formula, family, data, random, ...) with additional arguments
subset, method, na.action, control, and model, x logicals. I may try
reading the source code.
On the other hand, my need is sufficiently
Hi, Doug:
Thanks. I'll try the things you suggests. The observed
proportions ranged from roughly 0.2 to 0.8 in 100 binomial random
samples where sigma is at most 0.05. Jim Lindsey's glmm does
Gauss-Hermite quadrature, but I don't know if it bothers with the
adaptive step. With it,
Doug Bates wrote:
As George Box would tell us,
You have a big approximation and a small approximation. The big
approximation is your approximation to the problem you want to
solve. The small approximation is involved in getting the solution
to the approximate problem.
I asked Prof.
Hi
(B
(BI'm trying to use SJava and I have troubles.
(BI try to run examples from "Calling R from Java"
(Bbut,I have an error that
(B"fatal error: enable to open the base package"
(B
(BI heard SJAVA bug,
(Bso,could you send me your compiled SJava package with the modified
(BREmbed.c
I need to make plots similar to those produced by the s-plus rayplot function but
can't seem to find it in R. These 'vector maps' plot a ray or vector at each
specified location. Is there something similar in R ?
--Rich
Richard Kittler
AMD TDG
408-749-4099
Hi,
The problem is that I would like to do an interaction plot with
intervals based on Tukey's honestly significant difference (HSD)
procedure, but I do not know how to do it in R.
I have 3 factors A, B and C and a response variable response.
I would like to study a model where there are main
Hi,
I'm using the function clustalw in packages dna, but every time i have a segmentation
fault!
In your opinion What is the problem?Memory?
Please help me
Daniela
__
[EMAIL PROTECTED] mailing list
Might be a compilation problem, or change of system libraries.
[EMAIL PROTECTED] [EMAIL PROTECTED] writes:
Hi,
I'm using the function clustalw in packages dna, but every time i have a
segmentation fault!
In your opinion What is the problem?Memory?
Please help me
Daniela
Hi Daniela!
From where you got this package?
I even havent found it at cran.
It is not a bioconductor package either.
and
search.help(clustalW)
In any case it is better to contact the package provider directly.
Sincerely Eryk
*** REPLY SEPARATOR ***
On 6/8/2004 at 6:38 PM
At Monday 07:58 PM 6/7/2004, Richard A. O'Keefe wrote:
[snip]
There are three perspectives on programming languages like the S/R family:
(1) The programming language perspective.
I am sorry to tell you that the only excuse for R is S.
R is *weird*. It combines error-prone C-like syntax
Dear All,
I was writing a small wrapper to bootstrap a classification algorithm, but if
we generate the indices in the usual way as:
bootindex - sample(index, N, replace = TRUE)
there is a non-zero probability that all the samples belong to only
one class, thus leading to problems in the
Dear R People:
Please excuse the off topic question, but I
know that I'll get a good answer here.
If a single author is writing a journal article,
should she use We performed a test
or I performed a test,
please?
I had learned to use we without regard to the number
of authors. Is that true,
It's on Jim Lindsey's WWW site. It's a nice package. We've got an
extension that is somewhat incomplete that incorporates BLAST results
(from NCBI, though could config for local use) in a package called
BioSeq1 on Bioconductor's dev site -- I don't think that the CVS
viewer is available, but I
is there a good way to get the following fragment to work when calling
it as wrapper(1) ?
# cut here==
wrapper - function (choose=0)
{
x - seq(0,2*pi,len=100)
y - sin(1.5*x);
y - rnorm(y,y,.1*max(y))
if (choose==0) {
rm(fifu,pos=1)
fifu - function(w,x)
What's the preferred way in R for handling samples with binary data
(like chemical fingerprints encoded as hexadecimal strings with 0's and
1's indicating the absence or presence of chemical features) in methods
such as clustering and MDS. Do you always have to expand the fingerprint
data into
Hi all,
I'm using R to analyze some research and I'm not sure which test would be
appropriate for my data. I was hoping someone here might be able to help.
Short version:
Evaluate null hypothesis that change A1-A2 is similar to change C1-C2, for
continuous, non-normal datasets.
Long
Marc Schwartz [EMAIL PROTECTED] writes:
wget -r -l1 --no-parent -A*.gz -nd -P src/library/Recommended
http://www.cran.mirrors.pair.com/src/contrib/1.9.1/Recommended
The above _should_ be one one line, but of course will wrap here. There
should be a space between the two lines.
Kids these
On Tue, 2004-06-08 at 12:40, Peter Dalgaard wrote:
Marc Schwartz [EMAIL PROTECTED] writes:
wget -r -l1 --no-parent -A*.gz -nd -P src/library/Recommended
http://www.cran.mirrors.pair.com/src/contrib/1.9.1/Recommended
The above _should_ be one one line, but of course will wrap here.
I am trying to add some markers onto a contour map image I have created,
by using the text() function when I have already produced the map using
either image() or filled.contour(). For some reason the points appear to
be shifted considerably to the right of where they should be appearing,
despite
Have you considered qqplot(A1, A2) and qqplot(C1, C2)? If A2,
A2, C1, C2 are more like Poisson, I might try qqplot(sqrt(A1),
sqrt(A2)), etc.: Without the sqrt, the image might be excessively
distorted by largest values, at least in my experience.
hope this helps. spencer graves
No, not lazy evaluation. The explanation has to do with environments
and how callGeneric works.
It's an interesting (if obscure) example.
Here's the essence of it, embedded in some comments on debugging and on
style. (This will be a fairly long discussion, I'm afraid.)
A useful starting
On Tuesday 08 June 2004 10:43 am, you wrote:
If I understand you correctly, you have two set of ***paired***
data, one set from the A population, and one from the C population.
Form the pairwise differences:
A.diff - A1 - A2
C.diff - C1 - C2
Alas, they are not paired. A1 and
Hi everyone,
I'm interested in the analysis of spatial data, and I'm trying out several
R-packages.
Today I was attempting to use the package vardiag (version 0.1):
library(vardiag)
rs4.vo - varobj(rs4[,2:4],trace=2)
[1] 1
Error: couldn't find function nlregb
so far I know nlregb is a S-plus
Hi,
To speed up reading of large (few million lines) CSV files I am writing
custom read functions (in C). By timing various approaches I figured out
that one of the bottlenecks in reading character fields is the mkChar()
function which on each call incurs a lot of garbage-collection-related
It looks like ks.test might help, but it seems to me you need to
make parsimonious models of the change in distributions, so you estimate
a few parameters for the distributions of each of your 4 data sets with
standard errors for all the parameters you estimate. Then you can test
if the
I had assumed that the use of we in articles was either due to formality,
like the distinction between tu and vous in French. The English monarch
never refer to themselves in the singular. Or we as in both the author and
the reader. However a sample of size 4 of articles to hand suggests that
Hello all,
I'm working with different models that where implemented in DOS system each one, but
for doing comparisons between them I require to put some of the results in an unified
data form for all the models. I wonder if some of you have the knowledge of one
function in R that enable me
Hi
Laura Quinn wrote:
I am trying to add some markers onto a contour map image I have created,
by using the text() function when I have already produced the map using
either image() or filled.contour(). For some reason the points appear to
be shifted considerably to the right of where they should
?system
?shell (under Windows, which we must guess you are using if `DOS' means
MS-DOS).
Please use an informative subject line.
On Tue, 8 Jun 2004, Martínez Ovando Juan Carlos wrote:
I'm working with different models that where implemented in DOS system
each one, but for doing comparisons
Vadim Ogranovich [EMAIL PROTECTED] writes:
Hi,
To speed up reading of large (few million lines) CSV files I am writing
custom read functions (in C). By timing various approaches I figured out
that one of the bottlenecks in reading character fields is the mkChar()
function which on each
On Tue, 8 Jun 2004 12:23:58 -0700, Vadim Ogranovich
[EMAIL PROTECTED] wrote :
Hi,
To speed up reading of large (few million lines) CSV files I am writing
custom read functions (in C). By timing various approaches I figured out
that one of the bottlenecks in reading character fields is the
Hello again,
In a previous message I request your help, but I don't have been clear in my problem.
Specifically, I'm trying to create an interface in R for the X-12-ARIMA and TRAMO
SEATS, for the versions that run in MS-DOS. This problem awake in me the interest for
make interfaces to
Hello again,
In a previous message I request your help, but I don't have been clear in my problem.
Specifically, I'm trying to create an interface in R for the X-12-ARIMA and TRAMO
SEATS, for the versions that run in MS-DOS. This problem awake in me the interest for
make interfaces to
implemented in MS-DOS to. The question for you is if you know about
the existence of an R function that allows me to run -in Windows-
executable files in MS-DOS from the R command window.
Does system() do what you're looking for?
__
[EMAIL
Is there a limit to the number of columns that a data.frame can have? For
example, can I read.csv() a file that has 1000 columns and 10,000 rows,
will it break or is it limited by available memory.
...
Philip
All of R is available in source from www.r-project.org. This
includes software for ARIMA, etc.
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html;; it describes several things you can do that will get you closer to what you want and will help you formulate a
I am no expert in memory management in R so it's hard for me to tell
what is and what is not doable. From reading the code of allocVector()
in memory.c I think that the critical part is to vectorize
CLASS_GET_FREE_NODE and use the vectorized version along the lines of
the code fragment below
sorry to impose again.
At the default point size, R seems very good in selecting nice xlim/ylim
parameters, and leaving a little bit of space at the edges of its
xlim/ylim. alas, I now need to create ps graphics that can only occupy
a quarter of a page, so I need to blow up the text for
Vadim Ogranovich [EMAIL PROTECTED] writes:
I am no expert in memory management in R so it's hard for me to tell
what is and what is not doable. From reading the code of allocVector()
in memory.c I think that the critical part is to vectorize
CLASS_GET_FREE_NODE and use the vectorized version
On Wed, 2004-06-09 at 05:11, Erin Hodgess wrote:
If a single author is writing a journal article,
should she use We performed a test
or I performed a test,
please?
Does the particular journal specify a style manual or style sheet? The
Chicago Manual of Style \cite{Chicago2003}, p160, says
On Tue, 08-Jun-2004 at 08:56PM +0100, Phineas Campbell wrote:
| I had assumed that the use of we in articles was either due to
| formality, like the distinction between tu and vous in French. The
| English monarch never refer to themselves in the singular.
Hence the use of the term The Royal
thank you, marc. I will play around with these parameters tomorrow at
my real computer. yes, the idea is to just create an .eps and .pdf
file, which is then \includegraphics[0.25\textwidth]{} in pdflatex. I
need to tweak with the parameter ps.options(pointsize) because
otherwise, I end up
can we put a how to donate money to R on the R webpage? perhaps with
a paypal button?
even better, because I would like to donate some funds from my research
budget, could the R-project possibly sell some trinkets for a high price
for support? it is difficult to explain to a non-profit org
On Tue, 08 Jun 2004 21:18:34 -0400, ivo welch [EMAIL PROTECTED]
wrote:
And once I
do this, I need different R parameter defaults on the axes. With the
advice I have gotten, I think I am all set now. However, I am a little
bit surprised that noone has written a package around this
On Tue, 2004-06-08 at 20:18, ivo welch wrote:
thank you, marc. I will play around with these parameters tomorrow at
my real computer. yes, the idea is to just create an .eps and .pdf
file, which is then \includegraphics[0.25\textwidth]{} in pdflatex. I
need to tweak with the parameter
On Tue, 2004-06-08 at 21:02, Duncan Murdoch wrote:
On Tue, 08 Jun 2004 21:18:34 -0400, ivo welch [EMAIL PROTECTED]
wrote:
And once I
do this, I need different R parameter defaults on the axes. With the
advice I have gotten, I think I am all set now. However, I am a little
bit surprised
On Wed, 2004-06-09 at 13:22, ivo welch wrote:
can we put a how to donate money to R on the R webpage? perhaps with
a paypal button?
Does the R Foundation link meet this need?
http://www.r-project.org/foundation/main.html
Cheers
Jason
__
[EMAIL
I wrote:
The scope rules (certainly the scope rules for S) were obviously
designed by someone who had a fanatical hatred of compilers and
wanted to ensure that the language could never be usefully
compiled.
Drat! I forgot the semi-smiley!
Tony Plate [EMAIL
81 matches
Mail list logo