Angelo Secchi wrote:
Hi,
I'm using the ineq package to calculate some concentration measures
(Gini, Herfindal, ...) and I was wondering if there's around also a
function to calculate standard error on these measures. If not, is anybody
aware of where I can find a reference on this point?
Hi
I have a two-column integer matrix like this:
R jj
[,1] [,2]
[1,] -11
[2,] -22
[3,] -76
[4,] -87
[5,] -65
[6,] -98
[7,] -54
[8,]3 -3
[9,] -109
[10,] -43
I want a diagnostic that detects whether a row is a multiple of
the
On Mon, 7 Feb 2005, John Fox wrote:
Dear Shlomo,
Just last week I experienced a similar problem with RODBC and MySQL under
Windows, though the problem extended to both Rterm *and* Rgui. After a fair
amount of detective work, I discovered that the source of the difficulty was
a newly downloaded
Gabor == Gabor Grothendieck [EMAIL PROTECTED]
on Mon, 7 Feb 2005 04:39:40 + (UTC) writes:
Adrian Baddeley adrian at maths.uwa.edu.au writes:
: I want to equip a data frame with an attribute
: which specifies how to plot some of the columns.
:
: Up to now we
On Mon, 7 Feb 2005 04:39:40 + (UTC), Gabor Grothendieck
[EMAIL PROTECTED] wrote :
Adrian Baddeley adrian at maths.uwa.edu.au writes:
: I want to equip a data frame with an attribute
: which specifies how to plot some of the columns.
:
: Up to now we have been doing this by giving the data
Dear all,
I came across recently to the Rmatlab package
(http://www.omegahat.org/RMatlab/) that allows R and Matlab to talk
to each other. I made several attempts to install it on my Mac but
without any success as some headers seem to miss. I am just wondering
if anyone has already tried and
bonjour,
Je suis actuellement en derniere annee d'ecole d'ingenieur en informatique et
statistiques et je dois réaliser mon projet de fin d'études sur le logiciel R.
En fait, je dois réaliser un scoring sous R puis le meme sous SAS et comparer
les resultats. Mon fichier se prete à une
Help needed with lm function:
Dear R's,
Could anyone tell me how to replace the link function (probit logit,
loglog etc.) in lm
with an abitrary user-defined function? The task is to perform ML
Estimation of betas
for a dichotome target variable.
Maybe there is already a package for this (I did
nathalie bouez wrote:
bonjour,
Je suis actuellement en derniere annee d'ecole d'ingenieur en informatique et statistiques et je dois réaliser mon projet de fin d'études sur le logiciel R.
En fait, je dois réaliser un scoring sous R puis le meme sous SAS et comparer les resultats. Mon fichier se
see at ?glm and ?family
and use, e.g., `family=binomial(link=probit)'.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
Web:
I am not sure that fully answers Jeff's question. If the available link
functions (even with the quasi family) are not sufficient for your
needs, then you need to make your own constructor for a family object
- which provides the necessary information to the glm engine - and use
this as the
Hi,
you can see this link:
http://www.kb.u-psud.fr/kb/acces-etudiant/cours/biostat/html/Syntaxe%20et%20listings%20CUSSV/html/glm.htm
The R functions (ex. glm, polr) were compared to the SAS procédures (ex.
proc logistic, genmod).
perhaps the solution of your problem is here.
P.BADY
At
I'm learning to use Sweave (it's really great!) and everything works fine,
except that occasionally I get, after having sweaved many times,
--
Sweave(report.Snw)
Writing to file report.tex
Processing code chunks ...
1 : term hide
Error:
From: Martin Maechler [EMAIL PROTECTED]
maybe Adrian save()s that data.frame susequently?
Then, I assume the environment will copied.
In all(?) other circumstances that should only be a pointer and
not really use much memory.
Yes, sorry for the garbled message,
it is only when the object
On Mon, 7 Feb 2005 12:32:13 +0100,
Göran Broström (GB) wrote:
I'm learning to use Sweave (it's really great!) and everything works fine,
except that occasionally I get, after having sweaved many times,
--
Sweave(report.Snw)
Robin Hankin r.hankin at soc.soton.ac.uk writes:
:
: Hi
:
: I have a two-column integer matrix like this:
:
: R jj
:
:[,1] [,2]
: [1,] -11
: [2,] -22
: [3,] -76
: [4,] -87
: [5,] -65
: [6,] -98
: [7,] -54
: [8,]3 -3
:
On 07-Feb-05 [EMAIL PROTECTED] wrote:
Help needed with lm function:
Dear R's,
Could anyone tell me how to replace the link function (probit logit,
loglog etc.) in lm
with an abitrary user-defined function? The task is to perform ML
Estimation of betas
for a dichotome target variable.
Just for the record -- NEWS for 2.1.0 includes:
o binomial() has a new cauchit link (suggested by Roger Koenker).
the MASS polr for ordered response is also now adapted for the Cauchit
case.
url:www.econ.uiuc.edu/~rogerRoger Koenker
email [EMAIL PROTECTED]
Urs Wagner wrote:
Hello
I built a dll with cggwin.
When I want to load it with RGUI the RGUI
hangs. How can I troubleshoot it?
Not using cygwin, which is unsupported.
Please read the R for Windows FAQ, in particular Chapter 7 Building
from Source and follow the advices given therein - and given
Petr,
Thank you! Yes, rowSums appears to be even a little bit faster than
unique(which()), and it also maintains the original order. I do want
original order maintained, but I first apply a function to one of my
data.frames (that without any -s ... yes, these do represent nulls,
as someone
roger koenker rkoenker at uiuc.edu writes:
:
: Just for the record -- NEWS for 2.1.0 includes:
:
: o binomial() has a new cauchit link (suggested by Roger Koenker).
:
: the MASS polr for ordered response is also now adapted for the Cauchit
: case.
Do any of the data sets that come with
This works for your first example:
f - function(m) {
k - t(t(m) %/% m[1,])
which(rowSums(k - k[, 1]) == 0)
}
f(jj)
[1] 1 2 8
but not the second one. The 0s are problematic...
Andy
From: Robin Hankin
Hi
I have a two-column integer matrix like this:
R jj
[,1]
hi, i got just one question concerning R-Project. How can we simply
include in C++ program , the r library? can we find sources of examples.
Thanks for your help. Bests regards.
BILLY Frédéric
RESPONSABLE TECHNO ADI
Tél : 01.56.88.84.44
Fax : 01.56.88.84.05
[EMAIL PROTECTED]
On Mon, 07 Feb 2005 14:43:45 +0100, Uwe Ligges
[EMAIL PROTECTED] wrote :
Urs Wagner wrote:
Hello
I built a dll with cggwin.
When I want to load it with RGUI the RGUI
hangs. How can I troubleshoot it?
Not using cygwin, which is unsupported.
Please read the R for Windows FAQ, in
I would like to use pipes. The mkfifo call is not in the mingw
distribution that I am using.
Urs
Uwe Ligges wrote:
Urs Wagner wrote:
Hello
I built a dll with cggwin.
When I want to load it with RGUI the RGUI
hangs. How can I troubleshoot it?
Not using cygwin, which is unsupported.
Please read
Here is one way. It is the custom to return a value that will be assigned
to the variable, so I changed your 'macro' to a function that returns the
value and then assigns it to your variable:
test - function(name, value){
+ .result - NULL # initialize to NULL
+ .result[name] - value
How can I download the R1.5.O version for windows with its extend packages
Splancs 1.5.0 and Spatial 1.5.0?
-
mails !
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing
Hi all:
I have a question on how to go about creating a correlation matrix. I have
a huge amount of data21 variables for 3471 times. I want to see how
each of the variables correlate to each other. Any help would be
appreciated, including which package and which functions I should use to
Hi,
we try to do a logistic regression with the function glm.
But we notice that this function don't give the same results as the SAS proc
catmod (differents estimate given).
We try to change the contrast on R system with:
options(contrasts=c(unordered=contr.SAS,ordered=contr.poly))
We also try
On Mon, 7 Feb 2005, David Enot wrote:
I came across recently to the Rmatlab package
(http://www.omegahat.org/RMatlab/) that allows R and Matlab
to talk to each other. I made several attempts to install
it on my Mac but without any success as some headers seem
to miss. I am just wondering if
Have you considered cor? The command
'help.search(correlation)' suggests among other functions
var(stats), the documentation for which also includes cor.
If this is not adequate, PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html . It may help you get
Hi,
see ?cor in base package to get correlation matrix for
your data. Maybe it could be usefull getting principal
components (give a look to: ? princomp (base)) to
reduce the number of variables.
Hoping I helped you.
Best regards,
Vito
You wrote:
Hi all:
I have a question on how to go about
From the help for cor (from the stats package):
If x and y are matrices then the covariances (or correlations) between the
columns of x and the columns of y are computed.
So if you make a matrix with each column corresponding to one of your
variables, you can get what you're after.
For future
Gabor Grothendieck [EMAIL PROTECTED] writes:
Robin Hankin r.hankin at soc.soton.ac.uk writes:
:
: Hi
:
: I have a two-column integer matrix like this:
:
: R jj
:
:[,1] [,2]
: [1,] -11
: [2,] -22
: [3,] -76
: [4,] -87
: [5,] -65
CATMOD in SAS builds log-linear glms. If you are truly trying to create
logistic
regressions in SAS, use PROC LOGISTIC instead. If you really meant that you
are trying
to create log-linear models in R - then look up the usage of the function
loglin (e.g.
loglin())
Cheers,
Joe
Quoting [EMAIL
Le 07.02.2005 16:58, Jessica Higgs a écrit :
Hi all:
I have a question on how to go about creating a correlation matrix. I
have a huge amount of data21 variables for 3471 times. I want to
see how each of the variables correlate to each other. Any help would
be appreciated, including which
Dear All,
I' ve used ANCOVA (through lm(Y~factor*x)) to study the influence/differences
between levels of factor upon the Y-x relationship. I found that both
intercepts and slopes differ among levels of the factor. I understand from the
results of summary(lm(Y~factor*x)) which of the
Robin
The attached script works but it ain't pretty/
Martin
--
No virus found in this outgoing message.
Checked by AVG Anti-Virus.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
I've tried using cor() by the following sequence:
C - cor(x, y = NULL, use = all.obs, method = c(pearson))
where x is my matrix of 21 columns and 3471 rows.
and I get this error:
Error in cor(x, y = NULL, use = all.obs, method = c(pearson)) :
missing observations in cov/cor
any
Jessica Higgs wrote:
I've tried using cor() by the following sequence:
C - cor(x, y = NULL, use = all.obs, method = c(pearson))
where x is my matrix of 21 columns and 3471 rows.
and I get this error:
Error in cor(x, y = NULL, use = all.obs, method = c(pearson)) :
missing observations in
Le 07.02.2005 18:53, Jessica Higgs a écrit :
I've tried using cor() by the following sequence:
C - cor(x, y = NULL, use = all.obs, method = c(pearson))
In ?cor you find that :
If 'use' is 'all.obs', then the presence of missing observations will
produce an error.
try use=complete instead.
It looks like you have missing observations. With the use argument, you
can specify complete observations or pairwise-complete observations.
Hope this helps,
Matt Wiener
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Jessica Higgs
Sent: Monday,
Helene,
The problem is not with R, but with the specification of the
SAS procedure you used. It would be simpler and more
comparable if you specify your model in SAS with Proc Logistic
or Proc Genmod.
See if they work out for you,
Tim
Original message
Date: Mon, 7 Feb 2005 17:19:54
Sorry! Something went wrong with the attachment in my previous post.
Here it is.
---
intmult - function(X)
{
# extract columns with non-zero top cell
X1 - X[,which(X[1,]!=0)]
# look for integer multiples in reduced table
m - nrow(X1);n - ncol(X1)
D -
Gorjanc Gregor wrote:
Hello to Rusers!
I am puzzled with R and I really do not know where to look
in for my problem. I am moving from SAS and I have difficulties
in translating SAS to R world. I hope I will get some hints
or pointers so I can study from there on.
I would like to do something
Hi,
I've got a Dual G5 running 10.3.7 and I'm trying to install RMySQL.
I've already got R up and running. When I try the command
install.packages(RMySQL)
It downloads a few things and then produces the following error:
---
Configuration error:
Could not locate the
R comes with it's own copy of zlib so even if it is compiled into R,
there isn't necessarily a system-wide installation of the library.
You may still need to install it.
-roger
Xander Meadow wrote:
Hi,
I've got a Dual G5 running 10.3.7 and I'm trying to install RMySQL.
I've already got R up
Hi,
Thank you for responding. One thing I forgot to mention (although I'm
guessing you figured it out anyway) is that I'm running Mac OS X.
I checked and I've got the file /usr/include/zlib.h installed on my
machine. Is this enough for R or does it need a more robust
installation of zlib? If
Dear R-help,
Has anyone tried running R on a Beowulf-type cluster? I can get R to run in
batch (using R CMD BATCH) on a cluster, but am wondering if it is possible
to get an interactive R session on a compute node. Right now, if I run:
beorun --nolocal R
I just get the R start-up message
look into the following knot of CRAN packages:
SNOW
Rmpi
Rpvm
rlecuyer
rsprng
SNOW-FT
RScaLAPACK
et cetera.
these provide some pretty interesting bits - like a suite of parallelized
apply() functions.
--elijah
Has anyone tried running R on a Beowulf-type cluster? I can get R to run in
batch
Hi Elijah,
Thanks for the reply. I'm aware of almost all of those things you
mentioned, and have played with some of them. I'm not looking to do
distributed computing within R at this point. What I'm after, though, is
not really addressed by any of them. I want a user to be able to start an
Hi,
I am using R to plot the graph and the problem I am facing with my graph is
that I have lots of points concentrated in one area and It is creating a
visualization challenge.
Is there any commands in R I could use to solve this problem. Even if there is
no command, do you know how I could
R-DCOM server problem
My program (VB program) runs R-program through R DCOM
server. No problem to run statistical functions and my
own functions. However When I tried to load data using
ODBC I got an error message.
I need to set up any parameter for R-DCOM server for
ODBC ?
library(DBI);
ROracle ploblem
I can't install ROracle in the window XP. Is there any
way to install ROracle package on the windows OS ? I
think I can enhance data acquisition performance using
oracle native dirver instead of using ODBC.
Tae Sik Han.
North Carolina State University.
On Mon, 7 Feb 2005, Latha Raja wrote:
I am using R to plot the graph and the problem I am facing with my graph
is that I have lots of points concentrated in one area and It is
creating a visualization challenge.
How about 'plot( ,pch='.', )?
--
Robert Burrows, PhD
New England Biometrics
Hi, guys:
These days I keep using R and Weka to do data mining. I think my next
step is open the source codes so that I can customrize them and make
them better server my purpose. But now I kinda hesitate to do so b/c I
am really not sure which is better to start with. You know, both
require some
Latha Raja wrote:
Hi,
I am using R to plot the graph and the problem I am facing with my graph is
that I have lots of points concentrated in one area and It is creating a
visualization challenge.
Is there any commands in R I could use to solve this problem. Even if there is
no command, do you
Thank you all for your help :-) Much Much Aprreciated
I tried with the jitter command and it worked :-)
Cheers,
Latha
-Original Message-
From: Latha Raja
Sent: Monday, February 07, 2005 3:21 PM
To: 'r-help@stat.math.ethz.ch'
Subject: Need your help with
Glad that jitter() solved your problem. However, this would not be the case
for largeish bivariate displays in general. So may I add another
suggestion for those who might face similar problems: hexbin() in the hexbin
package on Bioconductor, which was specifically developed for the problem
you
Oi Bernardo
Maybe it would also be useful to show them some of the libraries available
for medical research. That list would also be useful for other people on
this forum including myself ;)
In the past I have made demos running R and other commercial software side
by side, to demonstrate
Hi R users!
I have a likelihood ratio statistic that depends on a parameter delta and I am
trying to get confidence intervals for this delta using the fact that the
likelihood ratio statistic is approx. chi-squared distributed.
For this I need to maximize the two likelihoods (for the ratio
Hello,
I have several data sets involving repeated measures over time of annual
grass survival in plots following application of several levels of
herbicide treatments. Since the data set contains quite a number of
zeros, it doesn't fit the assumptions for parametric analysis.
I read in
Shah,
[EMAIL PROTECTED] writes:
Hi R users!
I have a likelihood ratio statistic that depends on a parameter delta and I am
trying to get confidence intervals for this delta using the fact that the
likelihood ratio statistic is approx. chi-squared distributed.
For this I need to maximize the
Hello,
I would like to make a contourplot of the following data;
x - 1:10
y - 1:10
z - 100:110
By doing contour(x,y,z) I get the following error;
Error in contour.default(x, y, z) : no proper `z' matrix specified
How do I fix this??
Kind regards, Datius Blaszijk
[[alternative
Hi R folks,
I'm using Manova to evlaluate the effects of a prescribed fire treatment on
some forest stand structure data MODEL - manova(Y ~ FIRE, test =
Hotelling-Lawley) where Y is a 2 parameters describing tree distributions and
fire has 4 levels, 0,1,2,3.
Can anybody tell me how I can
On Tue, 8 Feb 2005 01:15:06 +0100 [EMAIL PROTECTED] wrote:
Hello,
I would like to make a contourplot of the following data;
x - 1:10
y - 1:10
z - 100:110
By doing contour(x,y,z) I get the following error;
Error in contour.default(x, y, z) : no proper `z' matrix specified
How
Hi all. Thanks for all of your help/suggestions. I found an old email in
the R-help archives, pieced together a couple things and arrived at the
solution below. As an additional followup, I thought I would go ahead and
post it should other readers come across this same situation. Here goes..
Dear R users,
can anybody explain the reason, why the first piece of code
below gives a parsing error, while the other two variations
work?
# Gives a parsing error
x - 1
if (x 0)
{
y - 1
}
else # Error occurs at this line
{
y - -1
}
# This works
x - 1
{
if (x 0)
{
On Tue, 8 Feb 2005, Christian Lederer wrote:
can anybody explain the reason, why the first piece of code
below gives a parsing error, while the other two variations
work?
# Gives a parsing error
x - 1
if (x 0)
{
y - 1
}
else # Error occurs at this line
{
y - -1
}
My first thought was that all it looked a bit complicated for something that
should be straightforward.
I created a file called t.txt. I worked out the way I would have done it and
then I tested to see which was fastest. One little hiccup is that the two
objects are not identical and I though
That's a shell/R script issue. I've come across it before, but I've
no longer got access to any clusters of that form (yet :-).
Depending on the jobs, it's not clear to me that you want a truly
interactive process on a subnode, but the idea is the same as the
chicanery we used to have a remote R
Hello!
Sorry, may be I will just say something trivial, but:
I can get R started (but the R prompt does not appear) and do computations,
but no access to an x11() device.
looks like you have started R with --slave switch and remote console has
no DISPLAY environment variable set (or your
Content-Type: text/plain; charset=iso-8859-1
Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate
permitted sender hosts)
X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch
Content-Transfer-Encoding: 8bit
X-MIME-Autoconverted: from quoted-printable to 8bit by
Tae Sik Han wrote:
ROracle ploblem
I can't install ROracle in the window XP. Is there any
way to install ROracle package on the windows OS ?
You can both, install from sources as descibed in the docs, as well as
downloading a binary version from the maintainer's side:
Helmut Kudrnovsky wrote:
Content-Type: text/plain; charset=iso-8859-1
Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted sender hosts)
X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch
Content-Transfer-Encoding: 8bit
X-MIME-Autoconverted: from quoted-printable
75 matches
Mail list logo