You're right!!!
Thank you very much
Selon Peter Dalgaard [EMAIL PROTECTED]:
[EMAIL PROTECTED] writes:
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).
Hi Peter, Gabor
On Feb 7, 2005, at 04:42 pm, Peter Dalgaard wrote:
Gabor Grothendieck [EMAIL PROTECTED] writes:
Robin Hankin r.hankin at soc.soton.ac.uk writes:
: I want a diagnostic that detects whether a row is a multiple of
: the first row or not.
..and in that case it might be easier just to
hello!
this is a question, how can i specify the random part in the GLMM-call
(of the lme4 library) for compound matrices just in the the same way as
they defined in the lme-Call (of the nlme library). For example
i would just need
random=list(my.Subject=pdBlocked(list(pdIdent(~...
Hi,
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED]
Sent: Friday, February 04, 2005 10:40 PM
To: Roger Bivand
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] simple example of C interface to R
On Fri, Feb 04, 2005 at 09:09:37PM +0100,
Liaw, Andy wrote:
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
Hi,
I'm using glm function to do logistic regression and now I want to know if it
exists a kind of R-squared with this function in order to check the model.
Thank you.
__
R-help@stat.math.ethz.ch mailing list
Dear list,
I would like to rename an object as follows:
SimLUall - matrix(c(0,1,0,0, 1,0,0,0, 0,0,1,0, 0,0,0,1),nrow=4, ncol=4)
j - 2
SimLUall2 - SimLUall and j
The value of j should be assigned automatically to SimLUall.
How can I achieve this?
Regards,
Ulrich
Ulrich Leopold wrote:
Dear list,
I would like to rename an object as follows:
SimLUall - matrix(c(0,1,0,0, 1,0,0,0, 0,0,1,0, 0,0,0,1),nrow=4, ncol=4)
j - 2
SimLUall2 - SimLUall and j
Not absolutely sure what you are going to do, I guess either
SimLUall2 - SimLUall
SimLUall - j
or
Hi Ulrich
See:
? assign and
? paste
with something like:
assign(paste(SimLUall,j,sep=),SimLUall)
Note though that if you are planning some simulations, many guRus would say
that you should use a list.
nsimul - 100
simOut - vector(mode=list,length=nsimul)
names(simOut) -
Pavel == Pavel Khomski [EMAIL PROTECTED]
on Tue, 08 Feb 2005 10:20:03 +0100 writes:
Pavel this is a question, how can i specify the random part
Pavel in the GLMM-call (of the lme4 library) for compound
Pavel matrices just in the the same way as they defined in
Pavel the
I guess Martin is struggling to have the most entries in the next version
of fortunes ;-)
Well... I hope so! Or we should really do something for new users - letting
any R core team member getting convulsive is dangerous for the futuRe.
Best wishes,
Eric
Pavel this is a question, how can i
Hi,
I don't know if a pseudo squared R for glm exists in
any R package, but I find some interesting functions
in S mailing list:
http://www.math.yorku.ca/Who/Faculty/Monette/S-news/0422.html
Here are some functions for calculating (pseudo-)R^2,
which may be of use to some.
Rsquared -
From: Barry Rowlingson
Liaw, Andy wrote:
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
I am a student and for my R project I have to program a memory game, that is a
game for the children with cards and you have to find pairs.
Wouldnt you have this program in your computer for make an example for me.
Thanks for helping me!!
__
On Tue, 08 Feb 2005 11:36:29 +0100, Eric Lecoutre
[EMAIL PROTECTED] wrote :
Note though that if you are planning some simulations, many guRus would say
that you should use a list.
nsimul - 100
simOut - vector(mode=list,length=nsimul)
names(simOut) - paste(simu,1:nsimul,sep=)
for (i in
Martin Maechler wrote:
Pavel == Pavel Khomski [EMAIL PROTECTED]
on Tue, 08 Feb 2005 10:20:03 +0100 writes:
Pavel this is a question, how can i specify the random part
Pavel in the GLMM-call (of the lme4 library) for compound
Pavel matrices just in the the same way as they defined
Terry Elrod posted functions for (pseudo-)R^2 on S-news back in 1998:
http://www.biostat.wustl.edu/archives/html/s-news/1998-02/msg00267.html
I am not sure that it performs well as a check of the model, so you
might want to investigate that further.
[EMAIL PROTECTED] wrote:
I'm using glm
Helene,
you should read up about AIC, deviance, and deviance
residuals, then look at the summary() for your model.
Dan
Dr. Daniel P. Bebber
Department of Plant Sciences
University of Oxford
OX1 3RB
UK
Message: 78
Date: Tue, 8 Feb 2005 11:15:35 +0100
From:
On Tue, 8 Feb 2005 12:35:17 +0100 (CET) [EMAIL PROTECTED] wrote:
I am a student and for my R project I have to program a memory game,
that is a game for the children with cards and you have to find pairs.
Wouldnt you have this program in your computer for make an example
for me. Thanks for
[EMAIL PROTECTED] wrote:
I am a student and for my R project I have to program a memory game, that is a game for the children with cards and you have to find pairs.
Wouldnt you have this program in your computer for make an example for me.
Thanks for helping me!!
!!! Generally, we don't do the
For OBVIOUS reasons, is there any chance that we could introduce
package() and deprecate library()?
(well, I'll also ask if we could deprecate = for assignment, but
that's hopeless).
best,
-tony
On Tue, 8 Feb 2005 11:49:39 +0100, Martin Maechler
[EMAIL PROTECTED] wrote:
Pavel == Pavel
Vito Ricci wrote:
Hi,
I don't know if a pseudo squared R for glm exists in
any R package, but I find some interesting functions
in S mailing list:
It is included in lrm in the Design package. But note that this is not
for checking fit but rather for quantifying predictive discrimination.
.
I need advice or opinions for the following problem:
In a sample a part of the respondents has experienced victimizing
events. Only a part of these events have been reported to the police.
The rate of events reported to the police is
number of experienced events / number of reported events.
A.J. Rossini [EMAIL PROTECTED] writes:
For OBVIOUS reasons, is there any chance that we could introduce
package() and deprecate library()?
usepackage() or usePackage() has been suggested, but someone got
ambitious and wanted it to be different from library(), and it sort of
didn't get any
Hello
Has anybody implemented a copula correlation in R in order to correlate
between two non-normally distributed time series.
Thanks
Sean
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hi,
it's dependend of your task and the amount of data.
R-project have a lot more possibilities manipulate the data and many
different plot possibilities.
Weka have the nice possibilty use the experimenter which could test
several algorithms
on different dataset's for comparison, further you
Hello all,
Ive been playing with nnet (package 'nnet') and Ive come across this
problem. nnet doesnt seems to like to have more than 1000 weights. If I
do:
data(iris)
names(iris)[5] - species
net - nnet(species ~ ., data=iris, size=124, maxit=10)
# weights: 995
initial value 309.342009
Douglas Bates wrote:
The GLMM function in the lme4 package allows you to specify crossed
random effects within the random argument without the need for the
pdBlocked and pdIdent constructions. Simply ensure that your grouping
factors are defined in such a way that each distinct group has a
Sorry for the very basic, general question I have read over the
methods aspects of the R Language Reference manual and some of the help
pages (methods, setGeneric, etc.) and remain a bit confused yet about R
methods. Looking through the contributed docs and mail archives didn't
give me
Dear useRs,
I have a script (Python) that every once in a while appends data to a
MySQL table. Meanwhile, I have a running R session, and I want it to be
aware of such table updates. I could write a loop in R to periodically
check whether new data has become available; however, are you aware of
a
A.J. Rossini wrote:
For OBVIOUS reasons, is there any chance that we could introduce
package() and deprecate library()?
(well, I'll also ask if we could deprecate = for assignment, but
that's hopeless).
Amen to both suggestions.
cheers,
Sean Davis wrote (inter alia):
Also, just a general question, but are S4 methods now the
standard?
My God! I hope not!!!
cheers,
Rolf Turner
[EMAIL PROTECTED]
The package blighty allows you to draw the UK coastline.
Regards
Mike
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of michael watson
(IAH-C)
Sent: 08 February 2005 13:25
To: r-help@stat.math.ethz.ch
Subject: [R] Drawing maps of UK
Hi
I have downloaded
You are not reading the help.
from the help file
## Default S3 method:
nnet(x, y, weights, size, Wts, mask,
linout = FALSE, entropy = FALSE, softmax = FALSE,
censored = FALSE, skip = FALSE, rang = 0.7, decay = 0,
maxit = 100, Hess = FALSE, trace = TRUE,
Hi all,
I develop and print from both Windows and Linux, and am seeing some
printing inconsistencies first described about a year and a half ago by
Andy Liaw (see below). Specifically, the line widths on my windows plots
are about 5 times smaller than that on Linux, and my windows printouts
do not
Dear all,
Can anyone help me plotting a panel plot with log-scaled x axes?
My formula looks like this:
coplot(response~div|treatment+grass,
panel=function(x,y,...){panel.xyplot(x,y,...);panel.lmline(x,y,...)})
And I´d like to have div plotted in log scale.
Thanks very much for your help!
Best
VR's S PROGRAMMING is an excellent reference with a thorough discussion and
examples of S3 methods. It is somewhat less thorough about S4 methods,
perhaps because the implementation was still in process at the time of the
book's writing. Maybe BDR's online complements have more, but I haven't
On Tue, 8 Feb 2005, dr mike wrote:
The package blighty allows you to draw the UK coastline.
Yes, everyone outside the US gets good, free data from the US government,
like coastlines, but has to jump through hoops to get administrative area
boundary files, that in the US are simply free for
I am using GRASS 5.7 and R 2.0 for OS X. I have a simple xy location I
am doing some theoretical work with and would like to import maps into
R using the interface in order to do analyses. Yet when I attempt to
do this, I get the following error:
G-gmeta()
Error in gmeta() : region for
Hello,
applying a function to a list of variables I face the following problem:
Let's say I want to compute tables for several variables. I could write a
command for every single table, like
bravo-c(1,1,2,3,5,5,5,);charly-c(7,7,4,4,2,1)
table(bravo); table(charly)
table(bravo); table(charly)
I understand that I need to have a (in this case) square matrix with all the
data. But the question now is;
- can the contourplot not interpolate the missing values
or alternatively
- I have fit a model to the z data (z = 100 + 0.5x + 0.5y). How can I make
from this model a square matrix z to
Hi,
I do have it installed on 2 Macs as well (OS X 10.2.8 and 10.3.7) and
what I need does work, however if you do the command check routine some
problems will likely be revealed. At least there were problems for me.
Denis
Le 08 févr. 2005, à 12:23, [EMAIL PROTECTED] a écrit :
De: Don
On Tue, 8 Feb 2005, Steven T.Stoddard wrote:
I am using GRASS 5.7 and R 2.0 for OS X. I have a simple xy location I
am doing some theoretical work with and would like to import maps into
R using the interface in order to do analyses. Yet when I attempt to
do this, I get the following
Thanks for your patch, I've modified the code of bandplot appropriately.
This change will be in the next release of the gregmisc bundle.
-G
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of R user
Sent: Tuesday, January 25, 2005 1:02 PM
To:
Those printers AFAIK support postscript.
How are you printing to them on Linux? I suggest you use dev.print under
Windows (it needs some setup, see ?postscript). That makes more sense
than going via PDF as the support is all already in R and it is AFAIK the
printer's native mode.
We've seen
Ok. The work-around I can grow accustomed to, or I might fall back on
GRASS5.4 if that works.
I wonder if another problem that cropped up is somehow related to the
first. Using interp=TRUE I managed to import a raster map to R. But
when trying to plot it, I get the following:
On Tue, Feb 08, 2005 at 10:55:01AM +0100, Rau, Roland wrote:
do you know already the page of Roger D. Peng?
He has a document entitled An Introduction to the .C Interface to R.
It is located at:
http://www.biostat.jhsph.edu/~rpeng/docs/interface.pdf
thanks. that's a nice tutorial, but
R People
Is it possible to make financial calculus with are. anybody can help me
on it?
TU
Adrián
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
__
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
white is supposed to be the default background for BMPs (according to
?bmp) but it doesn't work
bmp(test.bmp, bg=white)
plot(0)
dev.off()
# results in grey background
# This seems to be a good enough workaround for now.
bmp(test.bmp)
par(bg=white)
plot(0)
dev.off()
#
This might be easier for your purpose:
lapply(list(bravo=bravo, charly=charly), table)
$bravo
1 2 3 5
2 1 1 3
$charly
1 2 4 7
1 1 2 2
Andy
From: Heinz Tuechler
Hello,
applying a function to a list of variables I face the
following problem:
Let's say I want to compute tables
On Tuesday 08 February 2005 10:19, Christoph Scherber wrote:
Dear all,
Can anyone help me plotting a panel plot with log-scaled x axes?
My formula looks like this:
coplot(response~div|treatment+grass,
panel=function(x,y,...){panel.xyplot(x,y,...);panel.lmline(x,y,...)})
Are you sure that
If you pass your function a NAMED list, then the following works:
demofn-function(varlist)
{
nm-names(varlist)
for (i in seq(along=varlist))
{cat('\n',nm[i])
print(table(varlist[[i]]))
}
}
demofn(list(bravo=bravo, charly=charly))
If you don't pass a named list,
Hello,
I was wondering, is there a way to plot estimated betas and standard
deviations that I've generated?
Specifically, I have the following:
A matrix of estimated betas (two different betas, 4 different simulations -
categories).
A matrix of estimated standard errors corresponding
Hello,
I'm hoping that there is an easier way to rename columns in a data frame
other than by using the names() assignment, which requires you to type in
all the column names at once for a data.frame, in the case that I simply
want to rename a single column in a data frame.
Also, is there an
You could use Dotplot() in the Hmisc package. For example:
mydata - data.frame(BETA = runif(4, min=.50, max=.70), SIM = c(Sim A,
Sim B, Sim C, Sim D))
mydata$CIL - mydata$BETA - runif(4, min=.05, max=.20)
mydata$CIU - mydata$BETA + runif(4, min=.05, max=.20)
Dotplot(SIM ~ Cbind(BETA, CIL, CIU),
On Feb 8, 2005, at 1:47 PM, Ken Termiso wrote:
Hello,
I'm hoping that there is an easier way to rename columns in a data
frame other than by using the names() assignment, which requires you
to type in all the column names at once for a data.frame, in the case
that I simply want to rename a
Ken Termiso wrote:
Hello,
I'm hoping that there is an easier way to rename columns in a data
frame other than by using the names() assignment, which requires you
to type in all the column names at once for a data.frame, in the case
that I simply want to rename a single column in a data frame.
Hi,
I am trying to compile R as a shared library (need to run RMAGEML,
which depends on SJava) on 64bit SUSE Linux 9.1 . I am using the
source code for R 2.0.1 (Nov 15 04)
./configure R_PAPERSIZE=LETTER
R_BROWSER=/opt/kde3/share/applications/kde/konqbrowser.desktop
--enable-R-shlib
Duncan Murdoch [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
On Tue, 8 Feb 2005 11:35:22 -0600, Earl F. Glynn
[EMAIL PROTECTED] wrote :
In what version of R? I just tried in 2.01 and R-patched, and it was
fine. Maybe the problem is your viewer?
Sorry.
I thought I was using R
Lic. Adrián Cecotto - FACEA wrote:
R People
Is it possible to make financial calculus with are. anybody can help me
on it?
Depends on what you mean with financial calculus.
Please read the posting guide mentioned below and - if you have still
questions - ask a precise question.
Uwe Ligges
TU
Rtest.o is just an object file, and the second call to gcc is really a call
to ld and creates Rtest.so, a shared library. Neither is executable.
Probably you want to compile Rtest.c into an executable; it does have the
main function after all. Perhaps
gcc -L/usr/local/lib -o Rtest Rtest.o
The manual is the source code: we would have to read it to answer your
questions so you may as well read it yourself.
That driver was written last century by someone no longer active in the R
project.
On Tue, 8 Feb 2005, Earl F. Glynn wrote:
Duncan Murdoch [EMAIL PROTECTED] wrote in message
Did none of you notice the Makefile in that directory?
On Tue, 8 Feb 2005, Huntsinger, Reid wrote:
Rtest.o is just an object file, and the second call to gcc is really a call
to ld and creates Rtest.so, a shared library. Neither is executable.
Probably you want to compile Rtest.c into an
Hi,there
I'm doing some R batch jobs in Unix.
Something like
R prog1 output1 --save
R prog2 output2 --save
prog1 and prog2 are running at the same time and they are essentially same
except it contains different parameter values.
I was wondering if two processes will affect each other?
thanks for all the help. i've tried everyone's suggestions, to no
avail...
On Tue, Feb 08, 2005 at 04:19:10PM -0500, Huntsinger, Reid wrote:
Rtest.o is just an object file, and the second call to gcc is really a call
to ld and creates Rtest.so, a shared library. Neither is executable.
I'm doing some R batch jobs in Unix.
Something like
R prog1 output1 --save
R prog2 output2 --save
prog1 and prog2 are running at the same time and they are essentially
same except it contains different parameter values. I was wondering
if two processes will affect each other?
On Tue, 8 Feb 2005 [EMAIL PROTECTED] wrote:
thanks for all the help. i've tried everyone's suggestions, to no
avail...
Not *my* suggestion, though. Mine was to use the Makefile, which has
Rtest: Rtest.o embeddedRCall.o
$(R_CMD_LINK) -o $@ Rtest.o embeddedRCall.o $(LIBR)
and just needs
Dear all,
I try to test behavioural data durations whether they follow an
exponential distribution. My data have the particularity that there
are no values below a certain threshold, and the only test adapted to
this case I could find is the Shapiro-Wilk test for exponentiality
(Shapiro Wilk
At 12:51 08.02.2005 -0500, Liaw, Andy wrote:
This might be easier for your purpose:
lapply(list(bravo=bravo, charly=charly), table)
$bravo
1 2 3 5
2 1 1 3
$charly
1 2 4 7
1 1 2 2
Andy
Thank you for your answer. In the case of the example it is sufficient, but
I see that I did not choose
At 10:32 08.02.2005 -0800, Berton Gunter wrote:
If you pass your function a NAMED list, then the following works:
demofn-function(varlist)
{
nm-names(varlist)
for (i in seq(along=varlist))
{cat('\n',nm[i])
print(table(varlist[[i]]))
}
}
demofn(list(bravo=bravo,
On Tue, 8 Feb 2005, Chao Zhu wrote:
I'm doing some R batch jobs in Unix.
Something like
R prog1 output1 --save
R prog2 output2 --save
Consider R CMD BATCH instead: it redirects the error/warnings too, which I
expect you would want to do.
prog1 and prog2 are running at the same time and they
Perhaps this is useful:
QuantLib: http://www.quantlib.org/
Extensions: http://www.quantlib.org/extensions.shtml
RQuantLib: http://dirk.eddelbuettel.com/code/rquantlib.html
Cheers,
Jim
On Tuesday 08 February 2005 17:23, Lic. Adrián Cecotto - FACEA wrote:
R People
Is it possible to make
Joe Nocera wrote:
Helene -
In addition to some of the excellent suggestions already posited (e.g.
examining AIC,
pseudo R^2 in the Design package), you might want to consider another tool to
assess
logistic regression model accuracy: the area-under-curve (AUC) from a
receiver-operating
Here's a (rather ugly) function that takes its arguments and turn them into
a list, using the names of the arguments and the names of the components:
listWithName - function(...) {
aname - sapply(as.list(substitute({...})[-1]), deparse)
res - list(...)
names(res) - aname
res
}
I am a bit confused about how to get the format of the labels in xyplot
to show as dates rather than the numeric value.
e.g.,
x - as.Date(c(20/01/2001,20/02/2003,21/06/2004),%d/%m/%Y)
y - c(1,2,3)
plot(y~x) #produces formatted labels
xyplot(y~x) #doesn't give formatted x labels
On 9 February 2005 at 14:37, Denham Robert wrote:
| I am a bit confused about how to get the format of the labels in xyplot
| to show as dates rather than the numeric value.
|
| e.g.,
|
| x - as.Date(c(20/01/2001,20/02/2003,21/06/2004),%d/%m/%Y)
| y - c(1,2,3)
|
| plot(y~x) #produces
x - rep(seq(1,10),10)
y - rep(seq(1,10),each = 10)
z - matrix(100 + x*0.5 + y*0.5,ncol=10,byrow=T)
contour(z)
or back to where you were originally
x - 1:10
y - 1:10
# z - 100:110
z - 1:100
contour(matrix(z,ncol=10))
read the help on contour it states that
z: a matrix containing the values to
Dirk Eddelbuettel a écrit :
On 9 February 2005 at 14:37, Denham Robert wrote:
| I am a bit confused about how to get the format of the labels in xyplot
| to show as dates rather than the numeric value.
|
| e.g.,
|
| x - as.Date(c(20/01/2001,20/02/2003,21/06/2004),%d/%m/%Y)
| y - c(1,2,3)
|
thanks a lot for the information, reshape did the job
datars -reshape(data, timevar=TERRCODE, idvar=BID, direction=wide)
greetings
helli
BID TERRCODEANMCODE
200310413120660 22 0
200310413120660 273 0
200310413120660 280 0
200310413120660 467 0
200310413120660 468
80 matches
Mail list logo
