On Wed, 12 Oct 2005, David Zhao wrote:
Does anybody have experience in running jpeg device using xvfb-run on
linux? I've been having sporadic problem with: /usr/X11/bin/xvfb-run
/usr/bin/R --no-save Rinput.txt, with error saying: error in X11
connection. Especially when I run it from a perl
I suspect you mean `Quasi-Monte Carlo', but that is used to do integration
not simulation.
Using help.search led ne to
LowDiscrepancy(fOptions)
Low Discrepancy Sequences
QUnif(sfsmisc) Quasi Randum Numbers via Halton Sequences
both of which generate the
This is actually quadratic programming, so why do you want to use optim()?
There are packages specifically for QP, e.g. quadprog.
A more general approach is to eliminate one variable, which gives you an
inequality constrained problem in n-1 variables to which you could apply
contrOptim().
I am puzzled by the warning message in the output below. It appears
whether or not I fit the seasonal term (but the precise point of doing
this was to fit what is effectively a second seasonal term). Is there
some deep reason why AR parameters
(Warning message: some AR parameters were fixed:
Prof Brian Ripley wrote:
On Wed, 12 Oct 2005, David Zhao wrote:
Does anybody have experience in running jpeg device using xvfb-run on
linux? I've been having sporadic problem with: /usr/X11/bin/xvfb-run
/usr/bin/R --no-save Rinput.txt, with error saying: error in X11
connection. Especially
hy all
I wish to draw a graph against a date,
I have a set of date like this DD/MM/, corresponding to it a set of integer
, i wish to draw on x side the dates (the space between the dates have to be
constant, not based on the time between 2 dates but on the number of dates) and
on y side
Hello all,
A colleague at work set me the challenge to convert some MATLAB
code into R, to see which is faster. We'd seen that benchmark comparing
MATLAB 6.5 to R1.90 (and others), and so I thought that I should be able
to get roughly comparable speeds. The code has lots of multiplications
Hi all,
I want to make all possible combination from dataset below:
V1 - c(0,1,2)
V2 - c(0,1)
V3 - c(0,1)
V4 - c(0,1)
V5 - c(0,1)
V6 - c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
V7 - c(0,1,2,3,4,5,6)
V8 - c(0,1)
V9 - c(0,1)
V10 - c(0,1)
V11 -
On Thu, 13 Oct 2005, [EMAIL PROTECTED] wrote:
hy all
I wish to draw a graph against a date, I have a set of date like this
DD/MM/, corresponding to it a set of integer , i wish to draw on x
side the dates (the space between the dates have to be constant, not
based on the time between
The first thing is to ensure that you are using an optimized BLAS. On
Windows, use Goto's BLAS if you have it (is not currently available and
redistribution is not allowed) or one of the pre-built ATLAS-based
Rblas.dll on CRAN or (best of all) optimize your own build of ATLAS.
The Matrix
Thanks Prof. Ripley for your prompt reply. With regards to Rblas.dll my
current situation is that I have taken the Rblas.dll from CRAN:
contrib/ATLAS/P4/, and replaced the default Rblas.dll in my R /bin with
this one.
Cheers,
Rob.
-Original Message-
From: Prof Brian Ripley
Hi all R users,
I have problems with my second loop for drawing the three curves in the
same graphic. I need help please
Thank you in advance
#
simulation - function(k, n){
conc - seq(0,100,by=0.5)
#choixg - seq(1,
On Thu, 2005-10-13 at 12:19 +0200, Muhammad Subianto wrote:
Hi all,
I want to make all possible combination from dataset below:
V1 - c(0,1,2)
V2 - c(0,1)
V3 - c(0,1)
V4 - c(0,1)
V5 - c(0,1)
V6 - c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
V7 -
Dear all, Martyn Plummer and Jim Holtman (offlist) thanks you for quick
respons. Now I understand. I need more machine and memory.
Thanks a lot.
Muhammad Subianto
--- 20 columns and 54 billion rows? O:-)
On this day 13/10/2005 01:45 PM, Martyn Plummer wrote:
It's all about memory. In your
The Matrix package is under active development and the documentation
has not caught up with the code. Examples of usage can be found in
the tests subdirectory of the source package. At present we are
concentrating on the class hierarchy and writing methods and test
cases for those methods.
R-help,
I use the code below to plot some data by applying apply function.
But I don't know how I can get the argument type or col on the
plot function to distinguish the different lines
in the graph:
apply ( my.data, 2, function ( x ) lines ( dimnames ( my.data ) [[1]] ,
x ) )
Thank you in
I'm doing a big slow computation, and profiling shows that it is
spending a lot of time in match(), apparently because I have code like
x %in% listofxvals
Both x and listofxvals are factors with the same levels, so I could
probably speed this up by stripping off the levels and just treating
Hi dears,
I have difficulty to build a package !
I use Windows XP HOME, with and Intel PVI 2.66Ghz processor and 512MB of memory
I use RCMD check pyra1 and I got this :
* using log directory 'C:/DOCS/R/pyra1.Rcheck'
* using
* checking for file 'pyra1/DESCRIPTION' ... OK
* checking extension type
Sorry, a typo in my previous message (parens in the wrong place in the
conversion).
Here it is corrected:
I'm doing a big slow computation, and profiling shows that it is
spending a lot of time in match(), apparently because I have code like
x %in% listofxvals
Both x and listofxvals are
Hello
I'm fitting a gls model with a variance-covariance structure and an
getting an error message I don't understand
I'm using gls() from the nlme library with the structure defined by
correlation = corSymm(form = ~1|Subject), weights = varIdent(form=~1|strata)
I get the error
Error in
Hi,
I am puzzeled with a differing result of princomp in R and FACTOR in
SPSS. Regarding the amount of explained Variance, the two results are
the same. However, the loadings differ substantially, in the unrotated
as well as in the rotated form.
In both cases correlation matrices are analyzed.
I'd put the extra columns in their own data frame, and save that to disk
(use dates/times/process ids or some other unique identifier in the
filenames to distinguish them). When you need access to a mixture of
columns, load (or source, depending how you did the save) the columns you
need,
Hi friends, I'm beginning in R and I have simple question.
I have this piece of my program and how you see, that's ok (whit num-
0.002)
num-0.002 # ok, but not when I change whit num-0... ?
factor1-1;
while(1)
{
if (num*factor11)
factor1-factor1*10
else
{
print(out ok!!);
break;
}
}
[1] out
On Thu, 13 Oct 2005, justin bem wrote:
Hi dears,
I have difficulty to build a package !
I use Windows XP HOME, with and Intel PVI 2.66Ghz processor and 512MB of
memory
You need to review Appendix F of the R Installation and Administration
manual, and be sure you have installed the Windows
On Thu, 2005-10-13 at 16:13 +0200, Andreas Cordes wrote:
Hi,
I am puzzeled with a differing result of princomp in R and FACTOR in
SPSS. Regarding the amount of explained Variance, the two results are
the same. However, the loadings differ substantially, in the unrotated
as well as in the
Arturo Coral Alamo wrote:
Hi friends, I'm beginning in R and I have simple question.
I have this piece of my program and how you see, that's ok (whit num-
0.002)
num-0.002 # ok, but not when I change whit num-0... ?
factor1-1;
while(1)
{
if (num*factor11)
factor1-factor1*10
else
{
Andreas Cordes [EMAIL PROTECTED] writes:
Hi,
I am puzzeled with a differing result of princomp in R and FACTOR in
SPSS. Regarding the amount of explained Variance, the two results are
the same. However, the loadings differ substantially, in the unrotated
as well as in the rotated form.
ugh!
scan(what= does this...
thx anyway,
__
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
Sorry for the confusing term, just as Prof Ripley mentioned, I am referring to
do the integration via quasi-random number. Thanks.
Best regards,
Tony
Prof Brian Ripley [EMAIL PROTECTED] wrote:
I suspect you mean `Quasi-Monte Carlo', but that is used to do integration
not simulation.
Fair enough. To clarify what I'm trying to achieve I've pasted below a
small piece of the larger data frame with the hierarchical structure of
factors POPULATION and LOCID and the ascending order of YEARS and the
variable DBC that I would like to transform to another variable that is a
lag of
On Thu, 2005-10-13 at 14:50 +0100, Luis Ridao Cruz wrote:
R-help,
I use the code below to plot some data by applying apply function.
But I don't know how I can get the argument type or col on the
plot function to distinguish the different lines
in the graph:
apply ( my.data, 2, function
Dimitris: Thank you for the suggestion but I get an error just as when I
did similar commands using by(), The error given is
Error in $-.data.frame(`*tmp*`, LAGDBC, value =
tapply(csss3lagm81$DBC, :
replacement has 1089 rows, data has 8314
So I'm not sure what the problem is - why
On 10/13/2005 10:04 AM, justin bem wrote:
Hi dears,
I have difficulty to build a package !
I use Windows XP HOME, with and Intel PVI 2.66Ghz processor and 512MB of
memory
I use RCMD check pyra1 and I got this :
* using log directory 'C:/DOCS/R/pyra1.Rcheck'
* using
* checking for file
Le 12 Octobre 2005 18:03, Dirk Eddelbuettel a écrit :
i) In general, and especially between 'testing' and 'unstable', use
apt-pinning, explained in the apt-howto packages, esp apt-howto-en for
English; and on various places across the Net; try Google'ing for
apt-pinning.
Dear Dirk,
I'll jump
On Thu, 2005-10-13 at 10:02 -0400, Duncan Murdoch wrote:
Sorry, a typo in my previous message (parens in the wrong place in the
conversion).
Here it is corrected:
I'm doing a big slow computation, and profiling shows that it is
spending a lot of time in match(), apparently because I have
Hi,
I'm trying to get ...as a list of unevaluated arguments, ie.
substitute(list(...)) gives me an unevaluated list of the arguments,
but I want a list of the unevaluated arguments.
My attempts so far:
(function(...) substitute(...))(a=1, b=a) # Only returns first
(function(...)
Try this:
cl - as.list(match.call())
On 10/13/05, hadley wickham [EMAIL PROTECTED] wrote:
Hi,
I'm trying to get ...as a list of unevaluated arguments, ie.
substitute(list(...)) gives me an unevaluated list of the arguments,
but I want a list of the unevaluated arguments.
My attempts so
Is there any routine to generate for one vector matrix.
If I have X I want to generate start from zero to maximum value each vector.
For example, I have a vector x = (4,2,3,1,4)
I want to generate n=6 times, for 4, start 0 to 4, then 2 start 0 to 2, ect.
The result something like this:
Gabor Grothendieck [EMAIL PROTECTED] writes:
Try this:
cl - as.list(match.call())
or match.call(expand.dots=FALSE)$...
On 10/13/05, hadley wickham [EMAIL PROTECTED] wrote:
Hi,
I'm trying to get ...as a list of unevaluated arguments, ie.
substitute(list(...)) gives me an
On 10/13/2005 1:07 PM, Marc Schwartz (via MN) wrote:
On Thu, 2005-10-13 at 10:02 -0400, Duncan Murdoch wrote:
Sorry, a typo in my previous message (parens in the wrong place in the
conversion).
Here it is corrected:
I'm doing a big slow computation, and profiling shows that it is
On Thu, 2005-10-13 at 14:31 -0400, Duncan Murdoch wrote:
On 10/13/2005 1:07 PM, Marc Schwartz (via MN) wrote:
On Thu, 2005-10-13 at 10:02 -0400, Duncan Murdoch wrote:
Sorry, a typo in my previous message (parens in the wrong place in the
conversion).
Here it is corrected:
I'm
Have you looked at the g.data package? It might be useful
(but may still require some redesign of your dataset).
Greg Snow, Ph.D.
Statistical Data Center, LDS Hospital
Intermountain Health Care
[EMAIL PROTECTED]
(801) 408-8111
Ken Termiso [EMAIL PROTECTED] 10/13/05 08:14AM
I'd put the
Perfect! Thanks Peter and Gabor.
Hadley
On 13 Oct 2005 19:59:49 +0200, Peter Dalgaard [EMAIL PROTECTED] wrote:
Gabor Grothendieck [EMAIL PROTECTED] writes:
Try this:
cl - as.list(match.call())
or match.call(expand.dots=FALSE)$...
On 10/13/05, hadley wickham [EMAIL PROTECTED]
Hi,
I use the code below to aggregate / cnt my test data. It works fine,
but the problem is with my real data (33'000 rows) where the function
is really slow (nothing happened in half an hour).
Does anybody know of other functions that I could use?
Thanks,
Hans-Peter
--
dat -
On Thu, 2005-10-13 at 19:47 +0200, Jan Sabee wrote:
Is there any routine to generate for one vector matrix.
If I have X I want to generate start from zero to maximum value each vector.
For example, I have a vector x = (4,2,3,1,4)
I want to generate n=6 times, for 4, start 0 to 4, then 2
A not-so-clever solution:
n - 6
x - c(4,2,3,1,4)
y - matrix(0, nrow=n, ncol=length(x))
for (i in 1:n){
for (j in 1:length(x)){
y[i,j] - round(runif(1,min=0,max=x[j]))}}
Hope it helps
Xiaohua
On 10/13/05, Jan Sabee [EMAIL PROTECTED] wrote:
Is there any routine to generate for one vector matrix.
OK, I've tried to be a good citizen and use the searchable archives,
but with three search strings I haven't found the answer to what must
really be a simple question.
I want to create a list of objects from a set of matrices (in this
case, 300x300). Suppose the first matrix is A. I tried:
Try
x - list(A)
instead. Coercion (as.list) sees A as a vector of length 90,000 (forgetting
the dim attribute) and assumes you want to make it a generic vector,
which conceptually just changes its mode to list, whereas list
constructs a list with entries you pass as arguments.
Reid Huntsinger
Do Users of Nonlinear Mixed Effects Models Know
Whether Their Software Really Works?
Lesaffre et. al. (Appl. Statist. (2001) 50, Part3, pp 325-335)
analyzed
some simple clinical trials data using a logistic random effects
model. Several packages and
Hi all,
Here is a screenshot of the live cd someone mentioned here on the list.
I lanched it under the virtual pc qemu. the first thing that appears is th
R-php 0.99 in a firefox window (not shown on the screenshot).
http://img449.imageshack.us/my.php?image=encietylivecd1fc.png
--
Benedict
Yes, that is what I am looking for.
Many thanks to Marc Schwartz and Xiaohua for your help.
Sincerely,
Jan Sabee
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Ok so I see the problem that I'm having creating a new variable (LAG1DBC)
in the example data transformation below is that tapply() is creating a
list that is not dimensionally consistent with the data frame (data). So
how do I go from the list output of tapply() to create a dimensionally
Hi Spencer: Thanks for your interest! Also, the posting guide was helpful.
I think my problem might be solved if I could find a way to terminate nlm or
optim runs from within the user-given minimization function they call.
Optimization is unconstrained.
I'm essentially using normal like curves
On Thu, 13 Oct 2005, dave fournier wrote:
Internal Virus Database is out-of-date.
Talk about not being careful!
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861
On Thu, 2005-10-13 at 14:28 -0600, Brian S Cade wrote:
Ok so I see the problem that I'm having creating a new variable (LAG1DBC)
in the example data transformation below is that tapply() is creating a
list that is not dimensionally consistent with the data frame (data). So
how do I go from
Hi,
How can I execute some scripts from within R.
I have a large data file which I process (for instance with gawk, but not
only) before performing some statistics.
I would like to do this in R, so that I do not have to save many data files
and then making analysis on them (which proved to be
Marco,
use the system command.
?system
I hope that this helps,
Andrew
On Thu, Oct 13, 2005 at 05:14:38PM -0400, Marco Grazzi wrote:
Hi,
How can I execute some scripts from within R.
I have a large data file which I process (for instance with gawk, but not
only) before performing some
In a package, i type a function name and got the following message:
...
tmp - .Fortran(master, x = as.double(x), y = as.double(y),
sort = as.logical(sort), rw = as.double(rw), npd = as.integer(npd),
ntot = as.integer(ntot), nadj = integer(tadj), madj =
as.integer(madj),
Dnia czwartek, 13 października 2005 23:25, Andrew Robinson napisał:
Marco,
use the system command.
?system
I hope that this helps,
system(ls, intern = FALSE, ignore.stderr = TRUE)
Error in as.character(args[[i]]) : cannot coerce to vector
In fact it ignores the stderr !
--
pozdrawiam,
Dear all
I am using R to produce ordinations library(vegan) and the plot function
produced looks great on the screen but when I send it to jpg or pdf or eps
the resolution is not so good. Can you tell me how to get high resolution
images out of R for publication?
Cheers
Chris
Benedykt P. Barszcz wrote:
Dnia czwartek, 13 października 2005 23:25, Andrew Robinson napisał:
Marco,
use the system command.
?system
I hope that this helps,
system(ls, intern = FALSE, ignore.stderr = TRUE)
Error in as.character(args[[i]]) : cannot coerce to vector
In fact it
On Fri, 2005-10-14 at 00:04 +0200, Benedykt P. Barszcz wrote:
Dnia czwartek, 13 października 2005 23:25, Andrew Robinson napisał:
Marco,
use the system command.
?system
I hope that this helps,
system(ls, intern = FALSE, ignore.stderr = TRUE)
Error in as.character(args[[i]]) :
Benedykt,
No, there is no ls object in your workspace. Try this:
system(ls, intern = FALSE, ignore.stderr = TRUE)
Andrew
On Fri, Oct 14, 2005 at 12:04:52AM +0200, Benedykt P. Barszcz wrote:
Dnia czwartek, 13 pa?dziernika 2005 23:25, Andrew Robinson napisa?:
Marco,
use the system
Dave,
that's an interesting start for a comparison. Let me point out some
ways that you might construct a compelling argument. Of course, these
aren't exhaustive, and others may well provide further depth.
1) If I understand correctly, you're trying to estimate parameters
from a real
On Thu, 2005-10-13 at 15:20 -0600, Chris Buddenhagen wrote:
Dear all
I am using R to produce ordinations library(vegan) and the plot function
produced looks great on the screen but when I send it to jpg or pdf or eps
the resolution is not so good. Can you tell me how to get high resolution
Dnia piątek, 14 października 2005 00:13, Sundar Dorai-Raj napisał:
Don't you mean system(ls)? See ?system.
Arguments:
command: the system command to be invoked, as a string.
This is the kind of obstacles a newbie has to overcome. Whoeve is writing the
documentation for R, please do not
I am having some trouble getting the colors correct on county maps using
the maps package. I have a data.frame that contains coloration data for
every county -- it also contains a variable 'mapm' which fits the
'state,county' format used in the mapping package.
I use this to define colors
My thanks to Marc Schwartz who provided the solution - put unlist() around
the tapply() statement. Looks like it works.
Brian
Brian S. Cade
U. S. Geological Survey
Fort Collins Science Center
2150 Centre Ave., Bldg. C
Fort Collins, CO 80526-8818
email: [EMAIL PROTECTED]
tel: 970 226-9326
Actually one can download a working verssion of our software for free.
So anyonme can verify these results.
It is restricted by needing a network connection to get permission to
operate and is for evaluation only. The current version does not have
the Gauss Hermite
integration which I put in to
Le 13 Octobre 2005 18:39, Benedict P. Barszcz a écrit :
Dnia piątek, 14 października 2005 00:13, Sundar Dorai-Raj napisał:
Don't you mean system(ls)? See ?system.
Arguments:
command: the system command to be invoked, as a string.
This is the kind of obstacles a newbie has to overcome.
On Fri, 2005-10-14 at 00:39 +0200, Benedict P. Barszcz wrote:
Dnia piątek, 14 października 2005 00:13, Sundar Dorai-Raj napisał:
Don't you mean system(ls)? See ?system.
Arguments:
command: the system command to be invoked, as a string.
This is the kind of obstacles a newbie has to
First I have to correct myself; there probably is an ls object in your
workspace, but it's a function.
Second I have to suggest you should try to run the example code in the
help file. Surely this would clarify that the string has to be
quoted?
t1 - system(who, TRUE)
I'm not so sure that it's
Convert dat to a matrix and see if working with the
matrix instead of a data frame speeds things up
enough.
On 10/13/05, Hans-Peter [EMAIL PROTECTED] wrote:
Hi,
I use the code below to aggregate / cnt my test data. It works fine,
but the problem is with my real data (33'000 rows) where the
Gabor Grothendieck wrote:
Convert dat to a matrix and see if working with the
matrix instead of a data frame speeds things up
enough.
In the Hmisc package the asNumericMatrix and matrix2dataFrame functions
facilite this.
Also look at the summarize and mApply functions in Hmisc, which can be
Vincent,
On 13 October 2005 at 13:06, Vincent Goulet wrote:
| Le 12 Octobre 2005 18:03, Dirk Eddelbuettel a écrit :
| i) In general, and especially between 'testing' and 'unstable', use
| apt-pinning, explained in the apt-howto packages, esp apt-howto-en for
| English; and on various places
75 matches
Mail list logo