I have used gauss for a number of years. It is my experience that R is
slower than gauss, but on the other hand R has much more functionality
available than gauss and that saves programming time (and time to debug).
The trick is of course to program efficiently in R, just as is the case
when you
Bob,
on a remark about Brian Ripley's remark
if you're actually using indicies explicitly,
you probably haven't wrapped your head around
how powerful the indexing structure and
whole object approach is in S and in R,
Gabor Grothendieck wrote:
I believe the point was not so much that R
I think that is an excellent illustration of my point:
If you are writing code that works with single elements, you are
probably a lot better off writing C code to link into R (and C is
0-based ...).
but even in R it is not following
However, the R thinking is to work with whole objects
At 01:39 pm -0500 31/03/04, Raubertas, Richard wrote:
Another alternative is to use the underappreciated function
'sweep()':
sweep(A, 1:2, a, +)
I find the following helpful (it's not due to me but I cannot find
the original poster):
%.+% - function(a,x){sweep(a , 2:1 , x ,+ )}
%+.% -
Gabor Grothendieck wrote:
Now you can call f(x) and if x is a legacy variable f.legacy gets
used and if x is Oarray then f.Oarray gets used.
Thanks, Gabor. While I didn't understand the answer at all,
I expect that I will eventually and I appreciate the thought
that went into the affirmative.
Hallo all
I don't understand why the Search Engine Keywords
comand on the html help of R 1.8 doesn't work, is it
active? or the problem is in my machine?
thanks michele
__
The problem is with your machine (and some other people's). That search
engine does work on a properly configured machine.
Since haven't given us anything to go on (like which OS, browser, R
version (there is no such thing as R 1.8) ...) I can only suggest you look
at the R 1.9.0 beta which is
Thank you for your answers and suggestioins!
I just read a tutorial on the use of SAS with Discrete Choice Models by
Warren Kuhfeld (Marketing Research Methods in SAS, 2004, free for
download). He suggests using the Cox PH model. It necessitates reformatting
the data so that the alternatives
Hi,
Thanks for the pointer - the 'lm(y~x:z)' model does give the slopes
directly and hence confint gives the confidence intervals.
The thing that puzzles me is that my dummy data explicitly sets the
three levels of the factor to have different variances and yet the
standard error is the same
Bob Cain [EMAIL PROTECTED] writes:
Rolf Poalis wrote:
Dear R users,
Biostatistics Denmark would like to annouce the availability of the
new utility: sas2R --- a SAS to R parser.
Has this been done by anyone for Matlab?
I think there was a group at U.of Apraphulia working on it.
--
Dear mailing list,
Do you know a specific package in R where I can create an artificial survival
data with controlled censoring condition? (For example, I want to simulate a
data set with 5 variables and 20% of the observations are censored.)
Many thanks.
Jei
On Thu, 1 Apr 2004 07:19:49 +0200
Rolf Poalis [EMAIL PROTECTED] wrote:
Dear R users,
Biostatistics Denmark would like to annouce the availability of the
new utility: sas2R --- a SAS to R parser.
For almost 40 years SAS has been the primary tool for statisticians
worldwide and its
Have a look at function varclus in package Hmisc.
Best
Jens Oehlschlägel
--
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
for testing learning purposes I create X11 device and specify layout like
layout(c(1,2,3), 3, 1), so I could play with parameters and see
several plots at the same time. That works fine until I try to create 4-th
plot - all other plots erased.
That's expected behaviour: you asked
On Thu, 1 Apr 2004 15:38:54 +0400 (MSD), Oleg Bartunov
[EMAIL PROTECTED] wrote :
Hello,
for testing learning purposes I create X11 device and specify layout like
layout(c(1,2,3), 3, 1), so I could play with parameters and see
several plots at the same time. That works fine until I try to create
Rolf Poalis wrote:
Biostatistics Denmark would like to annouce the availability of the
new utility: sas2R --- a SAS to R parser.
For almost 40 years SAS has been the primary tool for statisticians
worldwide and its easy-to-learn syntax, unsurpassed graphical system,
powerful macro
On Thu, 1 Apr 2004, Philipp Pagel wrote:
for testing learning purposes I create X11 device and specify layout like
layout(c(1,2,3), 3, 1), so I could play with parameters and see
several plots at the same time. That works fine until I try to create 4-th
plot - all other plots erased.
On Thu, 1 Apr 2004, Duncan Murdoch wrote:
On Thu, 1 Apr 2004 15:38:54 +0400 (MSD), Oleg Bartunov
[EMAIL PROTECTED] wrote :
Hello,
for testing learning purposes I create X11 device and specify layout like
layout(c(1,2,3), 3, 1), so I could play with parameters and see
several plots at
On Thu, Apr 01, 2004 at 02:13:18PM +0200, Peter Dalgaard wrote:
Bob Cain [EMAIL PROTECTED] writes:
Rolf Poalis wrote:
Dear R users,
Biostatistics Denmark would like to annouce the availability of the
new utility: sas2R --- a SAS to R parser.
Has this been done by anyone for
The command:
layout(c(1,2,3), 3, 1) specifies 3 plots
Try
layout(1:4,2,2,byrow=T)
Regards,
Alex
-Original Message-
From: Oleg Bartunov [mailto:[EMAIL PROTECTED]
Sent: April 1, 2004 7:39 AM
To: R-help
Subject: [R] multiple plots problem
Hello,
for testing learning purposes I create
Have you considered using
?split.screen
and
?erase.screen
with a non-transparent background
instead of layout.
Example:
par(bg = white) # default is likely to be transparent
split.screen(c(3,1))# split display into three screens
screen(1)
plot(1:10)
screen(2)
plot(2:20)
screen(3)
interesting, it works, but still need to explicitly specify screen number.
It'b be nice to have tabbed canvas, but I suspect it's wrong list.
Oleg
On Thu, 1 Apr 2004, Angel Lopez wrote:
Have you considered using
?split.screen
and
?erase.screen
with a non-transparent background
Depends on what you understand by _explicitly_ specify screen number:
par(bg = white) # default is likely to be transparent
screenlist-split.screen(c(3,1))# split display into three
i-screenlist[1]
# Create a function to simulate the tabbing
nextplot-function(i)
{
i-i+1
if
Hello,
-Original Message-
From: Jingky P. Lozano [SMTP:[EMAIL PROTECTED]
Sent: Thursday, April 01, 2004 2:15 PM
To: [EMAIL PROTECTED]
Subject: [R] Question on Data Simulation
Dear mailing list,
Do you know a specific package in R where I can create an artificial
Angel,
thanks for example, but I think you don't understand me - I was looking
native way to say R, please, don't erase canvas. Any emulations required
additional typing (so many in R! - readline is great, but no support for
prefix search as in bash)) and some attention
(don't forget to type
I am trying to analyze a general mixed model, but am confused about
the model statement to use. The structure of the problem is the
following:
y = X b + Z u + e
where X and Z are known incidence matrices for fixed and random
effects respectively, b is a vector of fixed effects to be
This would be accomplished via the nlme library for mixed linear models. See Pinhiero
and Bates (2000) for all relevant documentation and data applications.
Harold
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
Sent: Thursday, April 01, 2004 10:21 AM
To: [EMAIL
On 31 Mar 2004 at 18:31, Bob Cain wrote:
Peter Dalgaard wrote:
I think it was said quite early in the thread, but since noone
apparently listened, let me reiterate: One of the powerful indexing
features in R is the negative index (all, except) and x[-0] would
lead to some
Correction:
I should have wrote
layout(matrix(c(1,2,3,4), 2, 2, byrow = TRUE))
Sorry
Alex
-Original Message-
From: Hanke, Alex [mailto:[EMAIL PROTECTED]
Sent: April 1, 2004 9:25 AM
To: 'Oleg Bartunov'; '[EMAIL PROTECTED]'
Subject: RE: [R] multiple plots problem
The command:
Ed L Cashin [EMAIL PROTECTED] writes:
...
Is there a way to tell aggregate just do perform median on column runtime to
select the whole row?
Some helpful folks have emailed me requesting more info about what I'm
trying to do. Here's a simple R function to produce a data frame like
the one I
Hello. I am trying to fit a non-rectangular hyperbola function to data
of photosynthetic rate vs. light intensity. There are 4 parameters that
have to be estimated. I find the nls function very difficult to use
because it often fails to converge and then gives out cryptic error
messages. I
As far as I can see, the problems with legacy functions kill the
class-based idea. As others have noted, if you define a class for which
the indexing function [ is zero-based, and then pass an object x of that
class to a legacy function, anything that function does that assumes 1 is
the first
When I try writing the current device to a jpeg file I keep getting a gray background
for the jpeg images using dev.print.
I have tried using the option bg=white and bg=#FF
but neither seem to have any effect.
If instead of writing to a screen device I open a jpeg device I can get any
Hi Bill,
I've just spent a few months trying to fit a model to a dataset, and it's not
easy. However, in my case, what appears to be recalcitrance on R's part
actually turns out to be well-founded warnings that the structure of the
model and the data are not permitting a clean, unambiguous
You have not shown us how you call dev.print() or the commands preceding
that. Did you specify the bg=white in the call to dev.print()?
Andy
From: John Muller
When I try writing the current device to a jpeg file I keep
getting a gray background
for the jpeg images using dev.print.
I
Hi there
Please excuse this elementary question, but I have been fumbling with this for
hours and can't seem to get it right.
I have a nested anova, with random factor lakefac nested within
factor fishfac (fixed), with an additional fixed factor Habfac. If I
consider everything as fixed
I get the gray background when I use plot to write to the screen
and then try to write the screen contents to a file using
dev.print(jpeg,
file=H:/Secondary/cycle_time/test.jpg,width=800,quality=100,bg=white) ;
I have also tried using bg = #FF instead of bg=white
and have also tried
Hi,
I am a brand new user of R and I have a really stupid problem: I am having troubles
applying fft on 2D matrices.
I import a data file of 256*256 elements with read.table (it is actually a grey scale
image, the corresponding data file being generated with Matlab), and I find that the
fft
Anne Gosset wrote:
Hi,
I am a brand new user of R and I have a really stupid problem: I am having troubles applying fft on 2D matrices.
I import a data file of 256*256 elements with read.table (it is actually a grey scale image, the corresponding data file being generated with Matlab), and I
On Thu, 1 Apr 2004, Anne Gosset wrote:
Hi,
I am a brand new user of R and I have a really stupid problem: I am
having troubles applying fft on 2D matrices. I import a data file of
256*256 elements with read.table (it is actually a grey scale image, the
corresponding data file being
On Thu, 1 Apr 2004 17:48:11 +0400 (MSD), Oleg Bartunov
[EMAIL PROTECTED] wrote :
On Thu, 1 Apr 2004, Duncan Murdoch wrote:
In Windows you can record all plots by using the record=TRUE option
when you open the graphics display, or from the menu on the display
window. You cycle back through
On Thu, 1 Apr 2004, John Muller wrote:
When I try writing the current device to a jpeg file I keep getting a
gray background for the jpeg images using dev.print.
That's because you are not using the correct tool. dev.print is designed
for printing, as it says.
I have tried using the option
par(bg=white)
plot(1:10)
dev.print(jpeg, width=600, height=600)
gives a white background on the jpeg copy. Without the first line it is a
transparent background, and the default matte for JPEGS is grey74, so the
viewer displays it as transparent over a grey matte.
BTW, John Muller has given
Hallo all
can someone explain me how the exogenus variables work
in the arimax models is not clear for me...
Thanks Michele
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
I read a data set into a data.frame. I then operate on the set, add
columns to the frame etc. and now want to output the frame into a file.
Trouble is, I want the file to have a name based upon the original data
set name.
I cannot figure out how to cat strings together and use that new string
as
On Thu, 1 Apr 2004, Mark O. Kimball wrote:
I read a data set into a data.frame. I then operate on the set, add
columns to the frame etc. and now want to output the frame into a file.
Trouble is, I want the file to have a name based upon the original data
set name.
I cannot figure out how to
Peter Wolf wrote:
sort.6-function(a){
n-length(a)
adapt-function(i){i+1} # local function to perform the index correction
a-c(0,a)
for(i in 2:n){
j-i-1
a[adapt(0)]-a[adapt(i)]
while(a[adapt(j)]a[adapt(0)]){
a[adapt(j+1)]-a[adapt(j)]
j-j-1
}
On Thursday 01 April 2004 01:40 pm, Mark O. Kimball wrote:
I read a data set into a data.frame. I then operate on the set, add
columns to the frame etc. and now want to output the frame into a file.
Trouble is, I want the file to have a name based upon the original data
set name.
I cannot
Sorry about the email issue. I have a spam filter that only let's things on
my pass list into my Inbox, all else gets put into a suspect folder.
I have added the R help list to my pass list.
Prof. Ripley asked
why if you are using Windows *are* you using dev.print
I actually started by using
I want to use R function Matrix inverse in my c code, please tell me how I can.
If there is a sample which can tell me how it works. It will be fantastic.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do
What in the world am I missing??
x-rnorm(20)
mean(x)
[1] -0.2272851
results-boot(x,mean,R=5)
results[2]
$t
[,1]
[1,] -0.2294562
[2,] -0.2294562
[3,] -0.2294562
[4,] -0.2294562
[5,] -0.2294562
Jeff Morris
Ortho-Clinical Diagnostics
A Johnson Johnson Co.
Rochester, NY
Tel: (585)
On Wed, Mar 31, 2004 at 08:32:49AM -0300, Luiz Rodrigo Tozzi wrote:
My question is: can I generate graphics and tables in gif ou any graphical
format through shell script?? can I call R, run a package in my ascii data e
then export the results to a gif, png or whatever?
On Wed, Mar 31,
On Thu, 1 Apr 2004 17:19:55 -0500 Morris, Jeffrey [OCDUS] wrote:
What in the world am I missing??
The help page?
As help(boot) tells you:
statistic: snip
In all other cases 'statistic'
must take at least two arguments. The first argument passed
will always be
You have to let boot pass the index for selection to the function as well.
So, for example, try:
boot.mean - function(data, index)
mean(data[index])
results - boot(x, boot.mean, R=5)
Andrew.
On Thursday 01 April 2004 14:19, Morris, Jeffrey [OCDUS] wrote:
What in the world am I
What you are missing is reading the helpfile!
The function mean is not a valid statistic to be passed to boot
function(x, i) mean(x[i])
is.
Please do read the documentation. That's why it is provided.
--
| Angelo J. Canty
What in the world am I missing??
stype. As in (from help(boot)
stype: A character string indicating what the second argument of
statistic represents. Possible values of stype are 'i'
(indices - the default), 'f' (frequencies), or 'w'
(weights).
This works.
I asked where index origin 0 would help.
(I am familiar with Dijkstra's article on why counting should start
at 0 and agree whole-heartedly. I am also convinced that kicking
against the goads is not helpful.)
Peter Wolf [EMAIL PROTECTED]
has provided two examples.
(1) Converting pseudo-code
tmp
row.labelsa b c
1 1 deadlift 7 13
2 2squat 7 24
3 3clean 7 10
4 4 deadlift 8 8
5 5squat 8 20
6 6clean 8 2
7 7 deadlift 9 5
8 8squat 9 32
9 9clean 9 19
Is there any effective way to get distinct geometric plotting symbols and
colors for plots involving more than 25 groups?
Thanks.
Dr. Marc R. Feldesman
Professor and Chairman Emeritus
Anthropology Department - Portland State University
email: [EMAIL PROTECTED]
email: [EMAIL PROTECTED]
fax:
Hi Michael,
I have a nested anova, with random factor lakefac nested within
factor fishfac (fixed), with an additional fixed factor
Habfac. If I consider everything as fixed effects, it's addmittedly
not the correct model, but I can at least get this to work:
...
So now I try to run it
Dear R People:
This is off topic, but this group is amazing:
There is a state object, with all kinds of variables about
US states. However, the numeric data is from 1977.
Does anyone know of a similar data set with more recent data, please?
Thanks in advance,
Sincerely,
Erin Hodgess
Associate
The answer is yes.
Have a look to pch and col plot parameters on
R-help.
A.S.
Alessandro Semeria
Models and Simulations Laboratory
Montecatini Environmental Research Center (Edison Group),
Via Ciro Menotti 48,
48023 Marina di Ravenna (RA), Italy
Tel. +39 544
62 matches
Mail list logo