Not an expert in programming either, but to me it seems like you've forgotten
to initialize the variable tr. It just picks up garbage from allocated memory
previously initialized by other processes.
Med venlig hilsen
Frede Aakmann Tøgersen
-Oprindelig meddelelse-
Fra: [EMAIL
You forgot to tell us you are running MacOS X. Looks like a problem with
your compiler installation (and are you using a pre-compiled R from CRAN
and an incompatible compiler?).
The installation procedure is complicated because the system dependencies
need to be checked. It works well under
Is your OS set to allow 3GB users address space for processes?
Since you are running a 32-bit executable under `Windows XP-64' I don't
know if that is possible or how to do it, so please ask whoever advised
you to buy that OS.
You might want to send some sample code to the mgcv maintainer to
Hello,
Check under par() the option srt and las.
Although the way to rotate the labels is via manipulating the axis() of
your plot separately.
Regards,
Carlos Ortega.
On 3/8/06, Lisa Wang [EMAIL PROTECTED] wrote:
Hello there,
I would like to plot a graph with the x axis's label displayed
On Wed, 8 Mar 2006, Liaw, Andy wrote:
[...]
Now, Roger Peng has kindly made available a note on debugging R code:
http://www.biostat.jhsph.edu/~rpeng/docs/R-debug-tools.pdf
There is also a new chapter in the `Writing R Extensions' manual in the
R-devel version of R (to be R 2.3.0). So now is
Please look at the help page for arima.sim(). You first argument is being
matched to 'rand.gen', and is not a function.
I suspect you meant to call arima().
On Tue, 7 Mar 2006, (s) Richard Nuttall wrote:
hi,
I am trying to fit an ARIMA model to some time series data, I have used
On Tue, 7 Mar 2006, hadley wickham wrote:
I tend to have to use trial and error myself. Here is another
possibility.
That's got the subsetting solved, so here's the next challenge
lm(x ~ y, z)
Error in as.data.frame.default(data) : cannot coerce class myobj
into a data.frame
you can also use str(data.frame) which allows you to access the structure of
your object.
Florence.
On 3/8/06, Federico Calboli [EMAIL PROTECTED] wrote:
Hi All,
is there an equivalent of the Unix command 'less' (or 'more'), so I can
look at what's inside a data.frame or a matrix without
Hi,
I am using the mAr package to calibrate an Multivariate model (size 3,
order 12). I am trying to do the two following things:
1. I would like to calibrate the model using not a single time
series, but several of them: each time series should be seen as one
independent realisation
Dear R-users,
I am trying to make my plot nicer. Does anyboby know how to make some
characters in the xlab of a plot to be superscripted or subscript? ?
for example m3/s in the x-axis label, how to make 3 supersripted?
Best,Jing
__
Dear all,
when making a DLL via Rcmd SHLIB is there a way to link against a
library such as Rblas (I am on a Windows platform) on a case to case
basis?
I played a bit around with some self-written C-code which should call
the function 'dasum' defined in Blas.h.
I encountered the following
Dear community,
if I install R locally as some user on a server without root rights, will R
be cabable of performing calculations on a cluster spreading the job on
several cpu's and nodes? Or will R just simulate that it is doing so while
just using one cpu in reality?
regards
Benjamin
Yes, lots of packages link to Rblas. The help page for SHLIB tells you
Please consult section 'Creating shared objects' in the manual
'Writing R Extensions' for how to customize it (for example to add
'cpp' flags and to add libraries to the link step) and for details
of
On Tue, 7 Mar 2006, Jason Horn wrote:
Thank you again Gabor, that did the trick. Any thoughts on where I
go go for a reference for these time codes? Where did you get
CDT6CST from? Or is this just one of those things that is common
knowledge in UNIX circles.
To the R developers: I
On 3/9/06 5:58 AM, Benjamin Otto [EMAIL PROTECTED] wrote:
Dear community,
if I install R locally as some user on a server without root rights, will R
be cabable of performing calculations on a cluster spreading the job on
several cpu's and nodes? Or will R just simulate that it is doing
Dear Prof. Ripley,
thank you very much. It works now absolutely fine using a Makevars
file with the contents you suggested.
Thanks again,
Roland
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Thursday, March 09, 2006 12:09 PM
To: Rau, Roland
Cc:
Jing Yang wrote:
Dear R-users,
I am trying to make my plot nicer. Does anyboby know how to make some
characters in the xlab of a plot to be superscripted or subscript? ?
for example m3/s in the x-axis label, how to make 3 supersripted?
Best,Jing
See ?plotmath.
plot(1:10, xlab =
Hi
see
?plotmath
# not tested
expression(paste(m^3, /s))
but I do not know how / is printed.
HTH
Petr
On 9 Mar 2006 at 11:48, Jing Yang wrote:
Date sent: Thu, 9 Mar 2006 11:48:45 +0100
From: Jing Yang [EMAIL PROTECTED]
To: r-help
Hello
I realized
Anova repeated mesur
And I have result following:
summary(a)
Error: indiv
Df Sum Sq Mean Sq F value Pr(F)
Alt 2 28.377 14.188 63.5236 2e-16 ***
Site 9 75.238 8.360 37.4279 2e-16 ***
Espèce1 1.002 1.002 4.4867
The problem with the simulation is with the second group where there
is a high probability of obtaining all zeroes for the sample and this
is causing problems with the Wald statistics you are using to check
for a difference. Here is an example.
Browse[1] summary(res)
Call:
glm(formula =
Michael -
I recall reading something Breiman wrote that said essentially don't skimp on
the number of trees - they are cheap to build and it makes for a better model.
Also, look at your error rates (using plot), and make sure you run enough
trees so that the error settles down. You'll likely
Hi
Sorry if this is in the help :-S
I've looked at example(dendrogram) and though it gives some indication of what
I want, it doesn't do all.
OK, so here is what I want to do: draw a tree, and then have an action, on
user-click, to either draw a sub tree or a plot of the data. I also want
Dear all,
I found a link to the Numerical Mathematics Consortium
http://www.nmconsortium.org/index.aspx
Their rationale can be summarized (as far as I understood) by the
paragraphs which I copied from their homepage and pasted below.
Maybe this could be of interest also for the development of R?
On 3/9/06 7:47 AM, michael watson (IAH-C) [EMAIL PROTECTED]
wrote:
Hi
Sorry if this is in the help :-S
I've looked at example(dendrogram) and though it gives some indication of what
I want, it doesn't do all.
OK, so here is what I want to do: draw a tree, and then have an action,
I'm guessing this is pretty much impossible to get around, because
there is no way to tell eval how to deal with myobj type objects, and
lm only dispatches based on the type of the first argument.
Did you write an as.data.frame method? From ?model.frame
My as.data.frame is :
dear list,
i have a matrix with missing values like
1234x
x2222
222x2
the x stands for the missing value.
i have to substitute it to NA
substitute or replace didn´t work out, as they are for vectors only
however matrix[,i] also didn´t work.
X - matrix(c( 1, 2, 3, 4, x,
x, 2, 2, 2, 2,
2, 2, 2, x, 2),
ncol=5, byrow=TRUE)
X.new - matrix(as.numeric(X), ncol=ncol(X))
or
X.new - apply(X, 2, function(x){replace(x, x == x, NA)})
Stefan Semmeling wrote:
dear list,
i have a matrix
One way:
A - matrix(round(rnorm(20)), ncol=4)
A[2,1] - x
A[which(A == x, arr.ind=TRUE)] - NA
A - matrix(as.numeric(A), ncol=dim(A)[2])
Best,
Matthias
-Ursprüngliche Nachricht-
Von: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Im Auftrag von
Stefan Semmeling
Gesendet: Donnerstag,
On Wed, 8 Mar 2006, Don MacQueen wrote:
I would have suggested US/Central, by analogy with US/Pacific
which I use on my unix-like Mac OS X 10.3.9 system. I don't know what
might be common knowledge in UNIX circles, but here's a bit of
information from my system.
It has a directory named
On Thu, Mar 09, 2006 at 08:44:27AM -0500, Chuck Cleland wrote:
ChuckX - matrix(c( 1, 2, 3, 4, x,
Chuck x, 2, 2, 2, 2,
Chuck 2, 2, 2, x, 2),
Chuck ncol=5, byrow=TRUE)
Why not:
X[X==x] - NA
mode(X) - numeric ## to get a numeric matrix
Btw if you
Dear r-helpers,
I can't get lsa to run because:
library(lsa)
Loading required package: Rstem
Error in library(pkg, character.only = TRUE, logical = TRUE, lib.loc
= lib.loc) :
'Rstem' is not a valid package -- installed 2.0.0?
In addition: Warning message:
cannot create HTML package
Thomas,
This is what I use for writing R programs. I use it with Linux but
adapting it for Windows shouldn't be that much of a problem for us MIT
guys.
In ~/.Rprofile (so it gets loaded every time I start R) I have
myedit - function(object) {
system(if [ ! -d $HOME/stat-misc/Rsrc ]; then
Michael--
I had emailed Duncan Lang about this a while back (the maintainer of Rstem).
He
gave me the following suggestion (which worked for me with R v2.2.1 and Windows
XP):
install.packages(Rstem, repos = http://www.omegahat.org/R;)
cheers, Dave
Dear r-helpers,
I can't get lsa to run
Can you create a small self contained reproducible example
that does not work? The reproducible example I provided earlier on
this thread worked fine.
One idea is to check what the class is of the output of
your .GGobiCall. If it were of class ggobiDataset then
it would in turn be calling
Hi,
I am very confused about those note that I found in
this web page about TYPE III SS in ANOVA. In my mind
Type III is not useful when there are interaction.
Can someone help me ?
web site :
http://www.xlstat.com/fr/support/tutorials/ano2.htm
__
I suppose you're referring to the sentence (or whatever it corresponds to in
the French version you cited):
This means that the order in which the variables are selected will not have
any effect on the values in the Type III SS. The Type III SS is generally
the best method to use to interpret
Hi Sean
Thanks for the help, but I really wanted to do this in R :)
Any suggestions?
Mick
From: Sean Davis [mailto:[EMAIL PROTECTED]
Sent: Thu 09/03/2006 1:01 PM
To: michael watson (IAH-C); r-help
Subject: Re: [R] Identifying or searching for labels in a
Prof Brian Ripley [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
On Wed, 8 Mar 2006, Liaw, Andy wrote:
So now is a very good time for
people to offer further information/suggestions for that chapter.
This may be useful:
Mark Bravington, Debugging Without (Too Many) Tears, R News,
Shown below is most of the FORTRAN subroutine named HCLUST.
My question concerns the DO loop labeled as '10'. What happened to its
CONTINUE statement? I will assume that after FLAG(I)=.TRUE. is executed that
control returns to DO 10 I=1,N. Am I correct?
Dave
C
David Emmith [EMAIL PROTECTED] writes:
Shown below is most of the FORTRAN subroutine named HCLUST.
My question concerns the DO loop labeled as '10'. What happened to its
CONTINUE statement? I will assume that after FLAG(I)=.TRUE. is executed that
control returns to DO 10 I=1,N. Am I
Can you create a small self contained reproducible example
that does not work? The reproducible example I provided earlier on
this thread worked fine.
I wish I could, and I'm very grateful for the help, but because the
data is an external pointer it's not easy to make a self-contained
example.
I have written some R functions and want to call them in C. Is there any way to
do it? I have read through the Writing R Extensionsdocument. It seems that
the only R functions that we can call in C is wrapped in Rmath.h or by
embedding R.dll. Thanks.
[[alternative HTML version
Dear all,
I try to figure out how to use R to count the number of pixels of the
same color in some gray-level picture. I managed to read it in either
tiff or jpeg format, but the returned pixmap object keeps its
information out of (my) reach. Is there an easy way to tabulate the
different
On Thu, 9 Mar 2006, Christian Jost wrote:
Dear all,
I try to figure out how to use R to count the number of pixels of the
same color in some gray-level picture. I managed to read it in either
tiff or jpeg format, but the returned pixmap object keeps its
information out of (my) reach. Is
I'm having this error when I try to output out a gene expression data frame
Error in file(file, ifelse(append, a, w)) :
unable to open connection
In addition: Warning message:
cannot open file 'SAMPLES.txt'
when issuing the command:
write.table(SAMPLES, file = SAMPLES.txt)
mark salsburg [EMAIL PROTECTED] writes:
I'm having this error when I try to output out a gene expression data frame
Error in file(file, ifelse(append, a, w)) :
unable to open connection
In addition: Warning message:
cannot open file 'SAMPLES.txt'
when issuing the command:
I think the problem is that your ggobiDataset objects are also
data.frame objects. They must NOT be. For example,
note how the following example fails once we add data.frame
to the class vector of x: The reason is that ordinary inheritance
is not used by model.frame; rather, it uses
Dear All,
I would be grateful if you can tell how can I suppress the legend
(automatically created) when I plot a groupedData.
suppose that I have data farme df_0 which contains the following columns:
ID represents the clusters, t represents time of observation, Y: the variable
Hi,
Having trouble loading tcltk in R 2.2.1 built from source.
./configure, make, make check, and make install run ok.
library(tcltk)
Error in firstlib(which.lib.loc, package) :
Tcl/Tk support is not available on this system
Error in library(tcltk) : .First.lib failed for 'tcltk'
On Wed, 8 Mar 2006, Paulo Brando wrote:
summary(model.fit) # just one species from one treatment shown below
Call: survfit(formula = Surv(time, censo) ~ treatment + species, data =
wsuv)
treatment=0, species=1
time n.risk n.event survival std.err lower 95% CI upper 95% CI
Patrice -
I had a very similar problem using TCL/TK 8.3.
Below is the email I sent to my computing group at work about how I
fixed it. Note that since my TCL/TK header (.h) files were in an odd
location, the first step probably isn't relevant for you. But I bet the
second step is.
Hi all,
I have a dataset which includes 84 rows and 4313 columns. Starting from the
2nd row, each row represents
a patient. The 1st column is for arrayID
2nd column is for time
3rd column is for cancer
4th column is for patientID
Hi All,
While attempting to build R-patched using the Portland Group compiler suite
on our dual Opteron 250 box running Scyld 29cz5 (based on RH, kernel
2.4.29), I get:
make[5]: Leaving directory
`/home/andy/Rbuild/R221-pgi/src/library/tools/src'
make[5]: Entering directory
On 3/9/06, Marc Bernard [EMAIL PROTECTED] wrote:
Dear All,
I would be grateful if you can tell how can I suppress the legend
(automatically created) when I plot a groupedData.
suppose that I have data farme df_0 which contains the following columns:
ID represents the clusters, t
For some reason I had to see my own post to realize my stupidity...
Defining PG_HOME _before_ running the configure solved the problem. Sorry
for wasting the bandwidth.
Andy
From: Liaw, Andy
Hi All,
While attempting to build R-patched using the Portland Group
compiler suite on our dual
In addition to the multinom(nnet) function mentioned below there is
some literature on how one can divide such polytomous problems into an
set of dichotomous classifications and then aggregate the results,
e.g.:
1) one-vs-all
2) pairwise comparisons (aka [double] round-robin) (Führnkranz)
Patrice Seyed [EMAIL PROTECTED] writes:
Hi,
Having trouble loading tcltk in R 2.2.1 built from source.
./configure, make, make check, and make install run ok.
library(tcltk)
Error in firstlib(which.lib.loc, package) :
Tcl/Tk support is not available on this system
Error
Thank you, this did the trick.
In which case
$ apt-get install r-cran-car
is your friend.
Is there anyway to upgrade my R install through CRAN? Or is there
some Debian repository that has an upgraded version of R?
Thanks again.
Jeremy
How to do a shift/rotate os a list?
if
a = c(1,2,3)
what is the best way to make a equal
c(3,1,2)
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
On 9 March 2006 at 20:47, Jeremy Morris wrote:
| Thank you, this did the trick.
|
| In which case
|
| $ apt-get install r-cran-car
|
| is your friend.
|
| Is there anyway to upgrade my R install through CRAN? Or is there
| some Debian repository that has an upgraded version of R?
On Thu, 9 Mar 2006, Omar Lakkis wrote:
How to do a shift/rotate os a list?
if
a = c(1,2,3)
what is the best way to make a equal
c(3,1,2)
a - c(a[length(a)],a[-length(a)])
or
n - length(a)
a - c(a[n],a[-n])
David Scott
_
a = c(1,2,3)
a
[1] 1 2 3
rev(a)
[1] 3 2 1
PS: a in your example is not a list; i.e class(a)
From: Omar Lakkis [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject: [R] shift / rota
Date: Thu, 9 Mar 2006 15:51:51 -0500
How to do a shift/rotate os a list?
if
a = c(1,2,3)
what is the best way
I'm looking to buy a new desktop which will primarily be used for
analyses of large datasets (100s of MB). I've seen postings from several
years back re the 'optimal' platform for running R, but nothing more
recently.
Specifically, I want to know: 1) if I run R under Windows, does having a
On 3/9/2006 4:47 PM, [EMAIL PROTECTED] wrote:
I'm looking to buy a new desktop which will primarily be used for
analyses of large datasets (100s of MB). I've seen postings from several
years back re the 'optimal' platform for running R, but nothing more
recently.
Specifically, I want
Hi all,
I have a very simple question that I can't seem to find the answer to.
How do I extract rows that meet a certain criteria from a data frame
and group them into a new data frame? For example, if I want to make a
new data frame that only includes rows of data for which the p values
(given
Dear List-mates,
I have been trying to set up a 4x8 trellis plot (that is, 4 columns, 8
rows), and would like to have the axis limits set based on the data range of
rows (for ylim) and columns (for xlim). I've been using the call:
foo-xyplot(y~x|Epoch+Subject,
type=c(l,r),
Matt,
Have a look at subset specially the examples
at the end.
I use it a lot.
Hope it helps,
Augusto
Augusto Sanabria. MSc, PhD.
Mathematical Modeller
Risk Research Group
Geospatial Earth Monitoring Division
Geoscience Australia (www.ga.gov.au)
You can try:
new.dataframe - my.dataframe[my.dataframe$p.value 0.05, ]
This will select all columns. Alternatively, you can specify the columns
that you want after the ,.
-Christos
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Matthew Scholz
Sent:
?subset
Marc
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Matthew Scholz
Sent: Thursday, March 09, 2006 2:18 PM
To: r-help@stat.math.ethz.ch
Subject: [R] newbie question: grouping rows
Hi all,
I have a very simple question that I can't seem to find
I am using R and LISREL for simulation studies. R generates the data
that is analyzed with LISREL.
In LISREL I use PV in the LISREL output statement to request estimated
variances. LISREL writes these in a file that looks like this:
1 0 0
0.100331D+01 0.144845D+01 0.141009D+01
The ouput from LISREL is a matrix,so the read.matrix(tseries) will do
the job in this situation.
?read.matrix
read.matrix(tseries) R Documentation Read Matrix Data
Description
Reads a matrix data file.
Usage
read.matrix(file, header = FALSE, sep = , skip = 0)
2006/3/10,
hello all, i'm forwarding this question for a colleague.
Is it possible to open a 2nd R Console?
regards,
mark+
--
mark garey
ucsf
department of epidemiology and biostatistics
division of biostatistics
185 berry street, suite 5700
san francisco, ca. 94107-1739
415.514.8147
On Thu, 9 Mar 2006, mark garey wrote:
hello all, i'm forwarding this question for a colleague.
Is it possible to open a 2nd R Console?
The answer is probably No, but since it isn't clear what your colleague
means by a 2nd R Console or what OS this is, it's hard to be sure.
-thomas
Greetings,
I have perused the r-help mailing list archives for an answer to this
question, without avail.
I would like to color the leaves of a dendrogram plot based on a cutoff in
one of the variables involved in the initial clustering.
My input data is in the form of:
B
I suggest you go to www.r-project.org - CRAN - (select a local
mirror) - Packages. Among the offerings there, you will find the
following:
distr Object orientated implementation of distributions
distrEx Extensions of package distr
distrSimSimulation classes based on package
Does anyone know of a resource for learning the basics of how to manage and
manipulate dates and times in R? I have been reading Introductory
Statistics with R by Peter Dalgaard which is fantastic. But alas, I could
find no reference to date and time. I have looked at the reference manual
but it
I'm familiar with Box and Tiao (1975) intervention analysis; I
studied time series under Box and Tiao. I don't know how to do that in
R, but there must be a way. Have you looked at the 'dse' bundle? That
comes with vignettes that make it relatively easy to learn (or at least
to
To manipulate date/time I use packages zoo and
survival.
Hope it helps,
Augusto
Augusto Sanabria. MSc, PhD.
Mathematical Modeller
Risk Research Group
Geospatial Earth Monitoring Division
Geoscience Australia (www.ga.gov.au)
Cnr. Jerrabomberra Av.
On Thu, 2006-03-09 at 22:38 -0500, Farrel Buchinsky wrote:
Does anyone know of a resource for learning the basics of how to manage and
manipulate dates and times in R? I have been reading Introductory
Statistics with R by Peter Dalgaard which is fantastic. But alas, I could
find no reference
On Windows, I often have more than one R console (instance of Rgui.exe)
open at the same time. Unless by mistake, each instance is open on a
unique working directory.
MHP
mark garey wrote on 3/9/2006 8:05 PM:
hello all, i'm forwarding this question for a colleague.
Is it possible to open a
[EMAIL PROTECTED] wrote on 3/9/2006 4:47 PM:
I'm looking to buy a new desktop which will primarily be used for
analyses of large datasets (100s of MB). I've seen postings from several
years back re the 'optimal' platform for running R, but nothing more
recently.
Specifically, I want to
On 3/9/06, Keith Chamberlain [EMAIL PROTECTED] wrote:
Dear List-mates,
I have been trying to set up a 4x8 trellis plot (that is, 4 columns, 8
rows), and would like to have the axis limits set based on the data range of
rows (for ylim) and columns (for xlim). I've been using the call:
Hi Dave,
On Mar 9, 2006, at 9:25 AM, Dave Atkins wrote:
install.packages(Rstem, repos = http://www.omegahat.org/R;)
This works.
install.packages(lsa)
Also works.
But:
library(lsa)
Loading required package: Rstem
Error in dyn.load(x, as.logical(local), as.logical(now)) :
unable
Dear Deepayan,
My deepest thanks! Your example code worked perfectly. I understood the
caveats you detailed about figuring out the axis limits before resolving the
layout, which depend on a bunch of other stuff. I'm glad this is possible.
Yes, this particular layout does only make sense in this
On Thu, 9 Mar 2006, Michael Kubovy wrote:
Does anyone have further advice?
Build Rstem from the sources on your own machine.
install.packages(Rstem, repos = http://www.omegahat.org/R;,
type = source)
should work if you have your compilers etc set up to build packages from
On Thu, 9 Mar 2006, Thomas Lumley wrote:
On Thu, 9 Mar 2006, mark garey wrote:
hello all, i'm forwarding this question for a colleague.
Is it possible to open a 2nd R Console?
The answer is probably No, but since it isn't clear what your colleague
means by a 2nd R Console or what OS this
Hi all,
Why cov(y, y) only gives one value, and cov(t(y), t(y)) gives 3x3 NA matrix?
Here my y is listed below and it is a 3x1 matrix.
I am expecting that if I have a random vector y=[y1 y2 y3]', here '
denotes a transposition so that
y is a column vector, where y1, y2, y3 are independent
On Thu, 9 Mar 2006, [EMAIL PROTECTED] wrote:
I'm looking to buy a new desktop which will primarily be used for
analyses of large datasets (100s of MB). I've seen postings from several
years back re the 'optimal' platform for running R, but nothing more
recently.
It is a subject which comes
Another solution is to use external tools. It really depends on what you
have to do: if it is to count pixels of one, or a few gray levels, it
could be fine to do it in R,... but if you want to count *all* pixels,
this could be more efficient using a C program, especially if you work
with
89 matches
Mail list logo
