Re: [R] browseURL question

2004-02-28 Thread Prof Brian Ripley
The problem is the interpretation by the shell used on Unix: on Windows no shell is used. It seems to me that the Unix version of browseURl should be quoting `url' in remoteCmd - if (isLocal) switch(basename(browser), gnome-moz-remote = , open = url, galeon = paste(-x,

RE: [R] LME, where is the package?

2004-02-28 Thread Ko-Kang Kevin Wang
Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ] On Behalf Of Yong Wang Sent: Saturday, February 28, 2004 8:48 PM To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Subject: [R] LME, where is the package? Hi, all My R 1.7.1 can not find lme, I just downloaded

Re: [R] logististic regression (GLM). How to get 95 pct. confidence limits?

2004-02-28 Thread Prof Brian Ripley
On Sat, 28 Feb 2004, Niels Steen Krogh wrote: Dear R-list. I'm doing af logistic analyses using gml. The model explaines variations in Adverse events infections (0 og 1) using age as explanatory variable. model2d-glm(formula=AEorSAEInfecBac~Age,family=binomial(logit),data=emrisk) I

Re: [R] logististic regression (GLM). How to get 95 pct. confidence limits?

2004-02-28 Thread Christian Schulz
Hi, in library(MASS) the function profile.glm should help you. christian Am Samstag, 28. Februar 2004 08:22 schrieb Niels Steen Krogh: Dear R-list. I'm doing af logistic analyses using gml. The model explaines variations in Adverse events infections (0 og 1) using age as explanatory

Re: [R] LME, where is the package?

2004-02-28 Thread Uwe Ligges
Ko-Kang Kevin Wang wrote: Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ] On Behalf Of Yong Wang Sent: Saturday, February 28, 2004 8:48 PM To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Subject: [R] LME, where is the package? Hi, all My R 1.7.1 can not find lme,

[R] Stepwise regression and partial correlation for wildlife census time series

2004-02-28 Thread Tomas Willebrand
Dear List; We are doing a time series analysis of wildlife census data. We use a stepwise regression of the annual per capita rate of increase against pervious years population size (log transformed) as suggested by Berryman Turchin (2001, Oikos 92:265-270). How can we obtain the partial

Re: [R] load data for mypkg-Ex.R

2004-02-28 Thread Prof Brian Ripley
On Fri, 27 Feb 2004, Roger Koenker wrote: On a related note is there a convention for cleaning up the detritus after running example(foo) I suppose sometimes users would like to have access to the objects that were created in the course of this, but perhaps more likely they

[R] when .Call can safely modify its arguments

2004-02-28 Thread Vadim Ogranovich
Hi, Writing R Extensions Guide clearly states that a C-function interfaced via .Call() should not modify any of its arguments. However I wonder if there are exceptions to this rule, i.e. when .Call can safely modify the arguments. For example when a function creates a list that is then populated

Re: [R] when .Call can safely modify its arguments

2004-02-28 Thread Prof Brian Ripley
On Sat, 28 Feb 2004, Vadim Ogranovich wrote: Hi, Writing R Extensions Guide clearly states that a C-function interfaced via .Call() should not modify any of its arguments. However I wonder if there are exceptions to this rule, i.e. when .Call can safely modify the arguments. For example

[R] Basic general statistical problem.

2004-02-28 Thread Fulvio Copex
Hello everyone, I'd like to have suggestions about a common basic statistical approach, hope to be useful also for other R beginners. When you first get some data (i.e length of river) you may want to look at its distribution. Then you probably want to find which law follows this distribution,

Re: [R] Basic general statistical problem.

2004-02-28 Thread Prof Brian Ripley
On Sat, 28 Feb 2004, Fulvio Copex wrote: Hello everyone, I'd like to have suggestions about a common basic statistical approach, hope to be useful also for other R beginners. Isn't this about statistics beginners? You seem to be assuming univariate data and a very specific distribution

[R] cluster-gruop-match with other attributes after na.omit

2004-02-28 Thread Christian Schulz
Hi, i want a cluster-analysis with clara, but getting an error because in cldat are NA's. Error in clara(cldat[, 1:3], 4) : Each of the random samples contains objects between which no distance can be computed. cldatx - subset(cldat,select=c(A,B,C)) cldaty - na.omit(cldatx) Now , clara

[R] Re: cluster-gruop-match with other attributes after na.omit

2004-02-28 Thread Christian Schulz
i get an idea: cldatx1 - subset(cldat,select=c(ID,A,B,C)) cldatx2 - subset(cldat,select=c(A,B,C)) The id didn't cause na.omit , so both have same obs. cl1 - clara(cldatx2,4) cluster - cl1$clustering match1 - cbind(cldatx1,cluster) resultdata - merge(cldat,macth1,by.x=ID,by.y=ID,all.x=T) Sure

RE: [R] matrix inverse in C

2004-02-28 Thread Liaw, Andy
In section 5.7 of the `Writing R Extensions' manual, it mentions the interface to numerical linear algebra routines (LINPACK, etc.). You will need to figure out which are the ones you need and how to call them. As these are expressed as Fortran routines, you can call them from C using the

[R] SVD/Eigenvector confusion

2004-02-28 Thread Philip Warner
My understanding of SVD is that, for A an mxn matrix, m n: A = UWV* where W is square root diagonal eigenvalues of A*A extended with zero valued rows, and U and V are the left right eigen vectors of A. But this does not seem to be strictly true and seems to require specific

RE: [R] importing S-Plus data files

2004-02-28 Thread Liaw, Andy
As I said, I believe at least in S-PLUS 6.1, the files in _Data are binary compatible across platforms. (As Prof. Ripley says, if that's the `data' you meant.) If the S-PLUS on Linux is so broken, and the Windows version can't read the data, then there's a good chance the data is corrupted

Re: [R] SVD/Eigenvector confusion

2004-02-28 Thread Douglas Bates
Philip Warner [EMAIL PROTECTED] writes: My understanding of SVD is that, for A an mxn matrix, m n: A = UWV* where W is square root diagonal eigenvalues of A*A extended with zero valued rows, and U and V are the left right eigen vectors of A. But this does not seem to be strictly

Re: [R] SVD/Eigenvector confusion

2004-02-28 Thread Philip Warner
At 01:01 AM 29/02/2004, Douglas Bates wrote: Is there a reason you are doing the SVD in such a complicated way? Why not use the svd function directly? I am using it to debug other code that is deigned to compute SVDs, so I actually want to understand the intermediate steps in constricting the

Re: [R] SVD/Eigenvector confusion

2004-02-28 Thread Prof Brian Ripley
On Sun, 29 Feb 2004, Philip Warner wrote: My understanding of SVD is that, for A an mxn matrix, m n: A = UWV* where W is square root diagonal eigenvalues of A*A extended with zero valued rows, and U and V are the left right eigen vectors of A. But this does not seem to be

Re: [R] SVD/Eigenvector confusion

2004-02-28 Thread Philip Warner
At 01:17 AM 29/02/2004, Prof Brian Ripley wrote: On Sun, 29 Feb 2004, Philip Warner wrote: My understanding of SVD is that, for A an mxn matrix, m n: A = UWV* where W is square root diagonal eigenvalues of A*A extended with zero valued rows, and U and V are the left right eigen

Re: [R] SVD/Eigenvector confusion

2004-02-28 Thread Prof Brian Ripley
On Sun, 29 Feb 2004, Philip Warner wrote: At 01:17 AM 29/02/2004, Prof Brian Ripley wrote: On Sun, 29 Feb 2004, Philip Warner wrote: My understanding of SVD is that, for A an mxn matrix, m n: A = UWV* where W is square root diagonal eigenvalues of A*A extended with zero

Re: [R] browseURL question

2004-02-28 Thread Gabor Grothendieck
One way to avoid this problem is to create a one line javascript file which redirects to your url. That way the shell never gets its hands on your url. After setting the url: url - http://search.yahoo.com/search?p=Bioconductorei=UTF-8fr=fp-tab-web-tn=20fl=0x=wrt; Run these three lines of

Re: [R] SVD/Eigenvector confusion

2004-02-28 Thread Spencer Graves
The documentation ?svd says U and V are orthogonal, i.e., that the transpose is the inverse. hope this helps. spencer graves Philip Warner wrote: At 01:01 AM 29/02/2004, Douglas Bates wrote: Is there a reason you are doing the SVD in such a complicated way? Why not use the svd function

[R] Stepwise regression and partial correlation for wildlife census time series

2004-02-28 Thread tomas . willebrand
Dear List; We are doing a time series analysis of wildlife census data. We use a stepwise regression of the annual per capita rate of increase against pervious years population size (log transformed) as suggested by Berryman Turchin (2001, Oikos 92:265-270). How can we obtain the partial

[R] questions about anova

2004-02-28 Thread Carlos Henrique Grohmann
Hello all, I have two questions about anova (one is probably VERY basic...) 1 - when one asks for a summary of a trend surface created with surf.ls, he/she gets: summary(g3r) Analysis of Variance Table Model: surf.ls(np = 3, x = gradiente$east, y = gradiente$north, z = gradiente$num1)

Re: [R] questions about anova

2004-02-28 Thread Roger Bivand
On Sat, 28 Feb 2004, Carlos Henrique Grohmann wrote: Hello all, I have two questions about anova (one is probably VERY basic...) About surf.ls() in the spatial package, and documented in Venables Ripley Modern Applied Statistics with S ... 1 - when one asks for a summary of a trend

[R] Stepwise regression and pacf

2004-02-28 Thread Erin Hodgess
Hello Tomas! There are functions for pacf and plot.acf. They are in library(ts) Hope this helps! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: [EMAIL PROTECTED] To: [EMAIL PROTECTED] Date: Sat, 28 Feb 2004

Re: [R] limit to complex plots?

2004-02-28 Thread Paul Murrell
Hi You could do this sort of thing easily using viewports in the grid package. And you can combine base plots with grid viewports using the gridBase package. Take a look at Volume 3(2) of R News for some examples and feel free to contact me directly if you have more questions; I'm always

Re: [R] writing polygons/segments to shapefiles (.shp) or other ArcGIS compatible file

2004-02-28 Thread Paul Murrell
Hi Patrick Giraudoux wrote: Hi, Thanks for all those information and the most valuable library you have developed. If somebody is to develop something to write shapefiles from polygon coordinates within R (most welcome), I don't think that the attribute file (dbf) will be an important issue.

[R] LCG with modulo 2^30

2004-02-28 Thread Enayetur Raheem
I can't run a function which generates random numbrers using linear congruential generator. My multiplier is a=5+8^6, increment is b=1 and modulo is m=2^30. the code I have written works for modulo upto m=2^28. For m= 2^29 , it says, can not allocate memory for the vector or something like

[R] Phase Plane

2004-02-28 Thread Erin Hodgess
Dear R People: Is there a function available to for phase plane plotting, please? Thanks, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: [EMAIL PROTECTED] PS R 1.8.1 on Windows XP 2000

Re: [R] LCG with modulo 2^30

2004-02-28 Thread Sean O'Riordain
Hi Raheem, Firstly - fair warning...I'm not an R expert at all! However it is my understanding that the expression for i in 1:m creates a full vector in memory of all the consecutive numbers 1 to m... (i presume these are 4-byte ints here otherwise it would have fallen over before 2^29), but

[R] Rcmd SHLIB

2004-02-28 Thread Gabriel Lawson
Ok, I think I may have a path or permissions problem (below). Anyone know which settings I should check? When I use Rcmd SHLIB filename I get: C:\Program Files\R\rw1081\binRcmd SHLIB info.diffusion.c process_begin: CreateProcess((null), dlltool -k --as as --dllname R.dll --def R. exp

Re: [R] LCG with modulo 2^30

2004-02-28 Thread Prof Brian Ripley
One quick addition: 1:2^31 fails because it is trying to create a sequence of integers, and the largest integer in R is 2^31-1. For the same sort of reason, the largest vector you can create is 2^31-1 elements, but on a 32-bit machine you will run out of address space at 2^29 doubles, and very

Re: [R] Rcmd SHLIB

2004-02-28 Thread Prof Brian Ripley
Most likely you don't have dlltool in your path. Further, I suspect you don't have the correct make in your path since AFAIK that does not use CreateProcess. Please check and check again that you have followed exactly all the instructions in the file readme.packages. See Q3.1 in the rw-FAQ.