On Wed, 9 Jun 2004, Zhu Wang wrote:
I have a problem to build a package on Windows XP while there is no
problem on Linux. The Makefile is something like:
There would be a problem on Linux, if that Makefile were used. I suspect
it is not used.
###
LIBNAME=cts
PKG_LIBS =
Rmetrics availalble from www.rmetrics.org has implemented several R
functions
which deal with column statistics, row statistics and rolling analysis.
Diethelm Wuertz
David Brahm wrote:
Mark Leeds [EMAIL PROTECTED] wrote (to S-News):
does anyone know the R equivalent of the SPlus rowVars
I am trying to compile R on an Itanium machine running Red Hat 7.2 and
gcc version 2.96.
The build fails with the following symptoms:
g77 -fPIC -g -O2 -c cmplx.f -o cmplx.lo
g77 -fPIC -g -O2 -c cmplxblas.f -o cmplxblas.lo
gcc -shared -L/usr/local/lib -o libRlapack.so dlapack0.lo dlapack1.lo
my major is bioinformatics, and i'm trying to cluster ( agglomerate
the closest pari of observations ) in R.
i have already got my own similarities metric, but do not know how to
clust it based on similarities instead of dissimilarities.
since the help document of hierclust mentions the parameter
Deepayan == Deepayan Sarkar [EMAIL PROTECTED]
on Wed, 9 Jun 2004 18:04:38 -0500 writes:
Deepayan On Wednesday 09 June 2004 01:58, Tom Mulholland wrote:
While using Lattice I received the following error.
Error in if (xx != 0) xx/10 else z/10 : argument is of length zero
On Thu, 10 Jun 2004, F. Tusell wrote:
I am trying to compile R on an Itanium machine running Red Hat 7.2 and
gcc version 2.96.
The build fails with the following symptoms:
g77 -fPIC -g -O2 -c cmplx.f -o cmplx.lo
g77 -fPIC -g -O2 -c cmplxblas.f -o cmplxblas.lo
gcc -shared
David == David J Netherway [EMAIL PROTECTED]
on Thu, 10 Jun 2004 11:19:05 +0930 writes:
David Thanks for all the relies.
David I recently discovered names and applied it to lm
David objects but did not think to apply it to the
David summary object.
ok, and now you should
Hi!
I would convert the simmilarities into dissimilarities by myself.
Its quite easy.
For example to make a dissimilarity from correlation which is a simmilarity measures
you can:
dcorr = 1 - cor
and dcorr will be a dissimilarity.
More general. Make the smallest value the largest and the
I apologize for posting this in essence the second time (no light at the
end of the tunnel yet..):
is there a way to enforce that nls takes both, the data *and* the
model definition from the parent environment? the following fragment
shows the problem.
# cut here==
wrapper -
As I've already said on the R-devel list,
MM why on earth are you using hierclust() from the ORPHANED package
MM 'multiv', when there's hclust() in the core 'stats' package
MM and 'agnes' in the recommended 'cluster' package ?
{and your question is not about hierclust but about
On Thu, 10 Jun 2004, joerg van den hoff wrote:
I apologize for posting this in essence the second time (no light at the
end of the tunnel yet..):
is there a way to enforce that nls takes both, the data *and* the
model definition from the parent environment? the following fragment
shows
I want to start R processes on multiple processors from single shell
script
and I want all of them to have different random seeds.
The usual way would be using the process ID for the initialization.
Or the shell script could just hand out consecutive numbers 1, 2, 3, ...
Lutz
Hi,
I would like to conduct a MANOVA. I know that there 's the manova() funciton and the
summary.manova() function to get the appropriate summary of test statistics.
I just don't manage to specify my model in the manova() call. How to specify a model
with multiple responses and one explanatory
On Thu, 10 Jun 2004, Prof Brian Ripley wrote:
Around R 1.2.x the notion was introduced that variables should be looked
for in the environment of a formula. Functions using model.frame got
converted to do that, but nls did not. I guess that the best way forward
is to ensure that nls (and
[Forwarding on behalf of a colleague]
She's got a list with several tables:
tab - list()
for(i in 1:6) {
+ tab[[i]] - table(freq[i])
+ }
tab
[[1]]
0 0.17 0.3 0.51 2.534
1962515252
[[2]]
0 0.17 0.31 2.534
199136
Hallo
On 23 May 2004 at 13:44, Stefanie von Felten wrote:
Hi,
I would like to conduct a MANOVA. I know that there 's the manova()
funciton and the summary.manova() function to get the appropriate
summary of test statistics.
I just don't manage to specify my model in the manova() call.
Prof Brian Ripley wrote:
On Thu, 10 Jun 2004, Prof Brian Ripley wrote:
Around R 1.2.x the notion was introduced that variables should be looked
for in the environment of a formula. Functions using model.frame got
converted to do that, but nls did not. I guess that the best way forward
is
On Thu, 10 Jun 2004, joerg van den hoff wrote:
Prof Brian Ripley wrote:
On Thu, 10 Jun 2004, Prof Brian Ripley wrote:
Around R 1.2.x the notion was introduced that variables should be looked
for in the environment of a formula. Functions using model.frame got
converted to do
Ko-Kang Kevin Wang wrote:
[Forwarding on behalf of a colleague]
She's got a list with several tables:
tab - list()
for(i in 1:6) {
Most easy solution:
Make freq[i] a factor with all the levels that may appear. Then all
tables have same dimension, e.g.:
tab - lapply(freq, function(x) table(
Hi!
First as Uwe Ligges explained and finally.
do.call(rbind,tab)
Sincerely
Eryk
*** REPLY SEPARATOR ***
On 6/10/2004 at 10:17 PM Ko-Kang Kevin Wang wrote:
[Forwarding on behalf of a colleague]
She's got a list with several tables:
tab - list()
for(i in 1:6) {
+
Hi there fellow R users,
R has many different clustering packages (e.g. mclust,cluster,e1071).
However, can anyone recommend a method to deal with data sets that contain
categorial and continuous variables?
Regards
Wayne
KSS Ltd
Seventh Floor St James's Buildings 79 Oxford Street
# Generate data
universe - c(0, 0.17, 0.5, 1, 2.5, 3, 4)
prob - c(94, 1, 1, 1, 1, 1, 1)/1000
tab - list(NULL)
for(i in 1:6){ tab[[i]] - table(sample(universe, 218, prob=prob,
rep=TRUE)) }
# Re-table
all - unique(unlist(sapply(tab, function(x) names(x
retab - t(sapply( tab, function(x)
Hi!
You need a apropriate dissimilarity measure.
look for daisy in package cluster
help(daisy,package=cluster)
x: numeric matrix or data frame. Dissimilarities will be
computed between the rows of 'x'. Columns of mode 'numeric'
(i.e. all columns when 'x' is a matrix) will
There must be a simple answer. I want to plot an expression,
where the expression is held in a string variable. The obvious
solution along the lines of
ex-x^2
plot( c(0,1), c(0,1), main=as.expression(ex) )
gives me the a title x^2 - ie doesn't treat it like an expression. I
suspect I don't
You *parse* strings to form expressions.
plot( c(0,1), c(0,1), main=parse(text=ex) )
On Thu, 10 Jun 2004, Steve Roberts wrote:
There must be a simple answer. I want to plot an expression,
where the expression is held in a string variable. The obvious
solution along the lines of
ex-x^2
Am 09.06.2004 22:23:51 schrieb(en) Peter Alspach:
2. How can I make outputs of the grafics (plot, hist, ...) into a
file?
what platform (Windows or Unix)?
Unix. (Gentoo Linux in special)
Martin
__
[EMAIL PROTECTED] mailing list
Hi
On 10 Jun 2004 at 12:28, Steve Roberts wrote:
There must be a simple answer. I want to plot an expression,
where the expression is held in a string variable. The obvious
solution along the lines of
ex-x^2
plot( c(0,1), c(0,1), main=as.expression(ex) )
Plotmath example led me to
Hi Rick,
learning R. A question to this end, though.
Would all of my existing S-Plus arrays, functions,
and so on have to be re-created from scratch in R,
or is there a way to copy them into the .RData
directory? The answer to this question has major
I am not sure if you have had this
I make a study in health econometrics and have a categorical dependent variable (take
value 1-5). I would like to fit an ordered probit or ordered logit but i didn't find a
command or package who make that. Does anyone know if it's exists ?
other question : i would like to make a least square
On Thu, 10 Jun 2004, Pikounis, Bill wrote:
Hi Rick,
learning R. A question to this end, though.
Would all of my existing S-Plus arrays, functions,
and so on have to be re-created from scratch in R,
or is there a way to copy them into the .RData
directory? The answer to this
On Thu, 10 Jun 2004, VUILLEUMIER Mathieu wrote:
I make a study in health econometrics and have a categorical dependent
variable (take value 1-5). I would like to fit an ordered probit or
ordered logit but i didn't find a command or package who make that. Does
anyone know if it's exists ?
On Thu, 10 Jun 2004, Prof Brian Ripley wrote:
On Thu, 10 Jun 2004, Prof Brian Ripley wrote:
Around R 1.2.x the notion was introduced that variables should be looked
for in the environment of a formula. Functions using model.frame got
converted to do that, but nls did not. I guess
Here is my best solution so far using $RANDOM in bash
or maybe somebody has pure R solution ?
#!/bin/bash
for i in 1 2; do
p=tmp$i
cat $p.R EOF
set.seed($RANDOM)
cat('rnorm(3) =', rnorm(3), '\n')
EOF
nice R CMD BATCH --no-save --no-restore $p.R $p.log
done
sleep 3
grep rnorm tmp[12].log
Best
No that isn't it I'm afraid - my expression is pasted together from
various bits and pieces and is at some stage a string variable. I
needed parse as Brian pointed out
Thanks anyway,
Steve.
From: Petr Pikal [EMAIL PROTECTED]
To: Steve Roberts [EMAIL
Dear Hank,
Last question first: really, only you can say for sure if 4e-281 and
5e-11 are small enough; it depends on the units you measure your state
variables in. However, this strategy cannot get the state variables to
exactly 0. Obviously, you could get closer to 0.0 faster by setting
perldoc -f srand
has a quick and good advice on initializing a random sequence.
I second the recommendation to prepare the list of random numbers
in advance for sake of reproducibility.
Itay
--
[EMAIL PROTECTED]
[EMAIL PROTECTED] writes:
Here is my best solution so far using $RANDOM in bash
or maybe somebody has pure R solution ?
Sure -- use the rsprng or rlecuyer packages from CRAN, which support
(independent?) streams of random numbers from a common set of seeds.
best,
-tony
--
[EMAIL PROTECTED]
Hi,
I would like to repsond to a few questions raised in a reply to a question
posted by another user on March 17, 2004. The entire message is copied at the
end of this email.
The relevant questions and statements are as follows:
What did you not understand about help(memory.size)?
This is
Does the 'groups' option on qqmath just color the points differently in the main
distribution or does it actually overlay separate quantile plots for each subset? I
would like to be able to do the latter.
--Rich
Richard Kittler
AMD TDG
408-749-4099
On Thursday 10 June 2004 09:21, [EMAIL PROTECTED] wrote:
Does the 'groups' option on qqmath just color the points differently
in the main distribution or does it actually overlay separate
quantile plots for each subset? I would like to be able to do the
latter.
The 'groups' option doesn't
I've started a tribe for discussing R and
sharing scripts. Tribe.net is one of the popular
on-line social communities, like Friendster.
Visit and see if it is a forum that you find useful.
To join the tribe you will need to register with
Tribe.net.
I hope it will be of help to newbies,
On Thu, 2004-06-10 at 06:40, Prof Brian Ripley wrote:
On Wed, 9 Jun 2004, Zhu Wang wrote:
I have a problem to build a package on Windows XP while there is no
problem on Linux. The Makefile is something like:
There would be a problem on Linux, if that Makefile were used. I suspect
it is
Hi folks,
I tried to compile a FORTRAN routine using the command g77 o2 c
all.f, it seems everything is fine, there is no error message at all,
but finally I can not get the file all.o. Therefore, I can not move on
to the next step using g77 shared c all.dll all.o to get the file
Thanks for the info. Given that I rule out qqmath would the best method be to make
repeated calls to qqnorm without plotting and then overlay the results or is there a
more elegant method?
--Rich
Richard Kittler
AMD TDG
408-749-4099
-Original Message-
From: Deepayan Sarkar
I am afraid your message is in an unreadable character set and none was
specified by your mailer. So I don't really know what you used, but I
have ^V where - should be, and so on.
What you would need is
g77 -O2 -c all.f
^ note
You appear to have asked for a file '2' here via -o2.
g77
I'm trying to use odesolve for integrating various series of coupled 1st
order differential equations (derived from a system of enzymatic
catalysis and copied below, apologies for the excessively long set of
parameters).
The thing that confuses me is that, whilst I can run the function rk4:
out
Hi folks,
I tried to use Mirosoft Fortran Powerstation 4.0 to create a dll file.
However, when I used the command dyn.load(test.dll), I got the
following message:
NULL
Warning message:
DLL attempted to change FPU control word from 9001f to 90003
I read the instruction on Duncan Murdochs
Rui wrote:
Hi folks,
I tried to use Mirosoft Fortran Powerstation 4.0 to create a dll file.
And what is the reason not to use the recommended gcc?
Uwe Ligges
However, when I used the command dyn.load(test.dll), I got the
following message:
NULL
Warning message:
DLL attempted to change FPU
I have a data frame called totaldata that is 10,000 rows by about 9 columns. The
data frame contains many zero entries which are important. If I type totaldata
it only prints out the first two columns. I expect (and want) it to print out all
9 columns. (I am actually sinking this to a text
On Thu, 10 Jun 2004, Rui wrote:
Hi folks,
I tried to use Mirosoft Fortran Powerstation 4.0 to create a dll file.
However, when I used the command dyn.load(test.dll), I got the
following message:
NULL
Warning message:
DLL attempted to change FPU control word from 9001f to 90003
I read
Rishi Ganti wrote:
I have a data frame called totaldata that is 10,000 rows by about 9 columns.
If about 9 equals 2, the behaviour reported below is expected.
The
data frame contains many zero entries which are important. If I type totaldata
it only prints out the first two columns. I expect
On Thu, 2004-06-10 at 11:26, Uwe Ligges wrote:
Rishi Ganti wrote:
I have a data frame called totaldata that is 10,000 rows by about 9 columns.
If about 9 equals 2, the behaviour reported below is expected.
That is, of course, for sufficiently large values of about...
;-)
Marc
On Thursday 10 June 2004 10:38, [EMAIL PROTECTED] wrote:
Thanks for the info. Given that I rule out qqmath would the best
method be to make repeated calls to qqnorm without plotting and then
overlay the results or is there a more elegant method?
I don't think there's any good way currently.
Chris Knight [EMAIL PROTECTED] writes:
I'm trying to use odesolve for integrating various series of coupled 1st
order differential equations (derived from a system of enzymatic
catalysis and copied below, apologies for the excessively long set of
parameters).
The thing that confuses me is
On Thu, 10 Jun 2004 09:12:05 -0700 (PDT), Thomas Lumley
[EMAIL PROTECTED] wrote :
DLL attempted to change FPU control word from 9001f to 90003
I read the instruction on Duncan Murdochs website about preserving
registers, but I still dont understand it. For example,
This code is for Delphi.
I always thought that the use of we was a reference to the
author(s) and the reader.
-roger
Erin Hodgess wrote:
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
Thank you very much to those who contributed to this rather interesting
discussion/debate. I was very surprised (and almost overwhelmed) by the
volume of replies based on this topic!
I have prepared some slides and put the draft version on
www.stat.auckland.ac.nz/~kwan022/tmp/RFin.ppt. Feel
Hi,
I have some code that looks like:
dftc - df[sets$tcset,]
pt - numeric(nrow(dftc))
sub - 1:nrow(dftc)
for (i in 1:nrow(dftc)) {
n - nnet( fmla, data=dftc, weights=wts, subset=sub[-i], size=4,
decay=0.01)
pt[i] - predict( n, dftc[ i, ], type='class' )
}
With tryCatch() it is possible to catch interrupts with tryCatch(). Then you
can use a 'finally' statement to clean up, release resources etc. However,
how can I protect against additional interrupts? This is a concern when
the hold down Ctrl+C and generates a sequence of interrupts. Example:
This looks like it probably is a scope problem with non-standard evaluation
rules for the argument subset= of nnet.
Instead of subset=sub[-i], try data=dftc[-i,] (I've not tested this since
I don't have the data objects you used.)
hope this helps,
Tony Plate
At Thursday 04:38 PM 6/10/2004,
Hi to you all
My question is:
there is a package written in UNIX for which there is no Mac version.
I would like to know if it's possible to install the R UNIX version on the
MacOSX and run that UNIX package on my Mac (through this UNIX R Vresion on
a Mac)
I have seen a porfile for r version
Hi:
Please tell us: which Mac OS are you using?
If you are using any of the 10.x, then you should download R-1.9.0 at:
http://cran.us.r-project.org/bin/macosx
If you have anything older (MacOS 9 or lower), then you should visit:
http://cran.us.r-project.org/bin/macos
Regards...
On Thu, 10 Jun
More complete support for interrupt management is still on the to do
list, unfortunately. onexit actions have similar issues (and have for
a long time). What I was hoping to do is provide a user level
mechanism for turning off interrupts around a piece of code and then
optionally enabling them
Hola ulisses
thank you very much for replying to my posting
I am using the last Mac Version:
OS X 10.3 Panther
I use R in the Mac and it works fine
The thing is that I want to use a package that was made in UNIX and I
wanted to know how easy it would be to run the UNIX version of R in the
Mac
Hola!
If you are already using Panther with R.dmg, then is pretty straight
forward to use a Unix package (panther it's a kind of Unix). If the
package is already on CRAN, all you have to do is to search the package
(either from the command line or from the menus) and then install. If
what
65 matches
Mail list logo