Re: [R] polr problem solved

2004-10-09 Thread Prof Brian Ripley
On Fri, 8 Oct 2004, Peter Flom wrote: I'd like to thank John Fox and Chuck Cleland for their help in resovling this issue. It turned out to be something simple, but perhaps others have had similar problems In my original data frame, I had 4 categories of race/ethnicity. One of the

Re: [R] Survey of moving window statistical functions - still looking f or fast mad function

2004-10-09 Thread Prof Brian Ripley
On Fri, 8 Oct 2004, Tuszynski, Jaroslaw W. wrote: [...] Finally a question: I still need to get moving windows mad function faster my runmad function is not that much faster than apply/embed combo, and that I used before, and this is where my code spends most of its time. I need

Re: [R] Adding factor variable to a CoxPH

2004-10-09 Thread Göran Broström
On Sat, Oct 09, 2004 at 09:03:23AM +0800, Neil Leonard wrote: I think this is a pretty basic statistic question: If I have a variable which has 4 factors how can I add it to a coxph model? I have other variables added which have digits 1-4 instead of four factors. If I recode the variable

[R] which() and value replacement in a matrix

2004-10-09 Thread Patrick Giraudoux
Hi, I cannot go through the archives with which() as key-word... so common. Though I am sure to have seen something about this subject in the past could somebody put me on the track. I have a matrix (actually a data.frame) in which I would replace the non-null values by 1. I tried the

Re: [R] polr problem solved

2004-10-09 Thread Peter Flom
Prof Brian Ripley [EMAIL PROTECTED] 10/09/04 3:18 AM asked How did you use `glm or lm' for an order factor response? An empty factor level will certainly cause glm problems, depending which one it is. An empty level will always cause polr problems, as there is no MLE under those

[R] Shangai ranking of world Universities

2004-10-09 Thread Jean lobry
Dear all, this is somewhat off-topic, but given the audience of this mailling list this might interest some people here. I have imported under R the data used in 2004 by Liu et al. from Shanghai Jiao Tong University to build their Academic Ranking of World Universities, see:

RE: [R] which() and value replacement in a matrix

2004-10-09 Thread Liaw, Andy
Use the index vector directly, rather than breaking it up: x - matrix(sample(30), 10, 3) idx - which(x 25, arr.ind=TRUE) idx row col [1,] 6 1 [2,] 9 1 [3,] 4 2 [4,] 6 2 [5,] 4 3 x[idx] - 999 x [,1] [,2] [,3] [1,]7 14 16 [2,] 20 248 [3,] 17

RE: [R] problem Installing R 2.0.0 on SuSe 9.1

2004-10-09 Thread Christopher Albert
Hello: I first tried it by using ./configure with the source files, but found that there is a huge problem with the gcc-g77 package, then I tried using the RPM obtained from CRAN, getting the following, the first of which appears to be another manifestation of the g77 problem: linux:/tmp

[R] RE: zero-inflated count models (was polr problem solved)

2004-10-09 Thread John Fox
Dear Peter, -Original Message- . . . The analyses were part of a paper I am writing, illustrating that, when the DV is oddly distributed (the DV in question was a count, with many 0's, and a long right tail) that the 'usual' methods not only are wrong for statisically reasaons

Re: [R] which() and value replacement in a matrix

2004-10-09 Thread Patrick Giraudoux
Thanks for the clue. Actually the trouble comes when refering to a data.frame. If I use the matrix from the data.frame (matrix(mydataframe)), everything goes smoothly... So I wrote: indices-which(myforetbin 0,arr.ind=T) myforetbin-as.matrix(myforetbin) myforetbin[indices]-1

[R] conversion of a data.frame of numerics to a data.frame of factors

2004-10-09 Thread Patrick Giraudoux
Hi, I am trying to convert a data.frame of numerics (this could be a matrix as well in this case) into a data.frame of factors. I did it in a way that is less than direct... myforet2-t(myforet) for (i in 1:length(myforet2[1,])) { if (i == 1)myforetfact-list(as.factor(myforet2[,i]))

[R] inst directory

2004-10-09 Thread Bob Wheeler
R CMD check now balks at my inst directory. It contains a single folder doc, but apparently any folder causes the problem. If inst is empty, the project checks OK. This was not a problem before 1.9. I've checked the documentation, but don't see a change. What am I missing. -- Bob Wheeler ---

Re: [R] conversion of a data.frame of numerics to a data.frame of factors

2004-10-09 Thread Prof Brian Ripley
On Sat, 9 Oct 2004, Patrick Giraudoux wrote: mydf[] - lapply(mydf, as.factor) would appear to be what you want. For a matrix, I presume you want a factor matrix as the result. Something like my - matrix(1:12, 3, 4) dmy - dim(my) my - as.factor(my) dim(my) - dmy my which does not print as a

[R] [R-pkgs] Publishing R package descriptions in JSS

2004-10-09 Thread Jan de Leeuw
More and more R packages come with a corresponding article in the Journal of Statistical Software (www.jstatsoft.org). Achim Zeileis, our TeXnical Editor, has recently contributed jss style files for issues, bookreviews, software reviews, and code snippets. They can be downloaded from

[R] [R-pkgs] verfication package announcement

2004-10-09 Thread Matt Pocernich
The verification package has recently been posted to CRAN. This package was initially developed by people in the Verification Group at the National Center for Atmospheric Research to verfiy and study weather models and forecasts. It has been written in general terms to be applicable to other

[R] [R-pkgs] R interface for MINPACK least squares optimization library

2004-10-09 Thread Timur Elzhov
Hello guys. I've built and uploaded to CRAN an R interface to MINPACK Fortran library, which solves non-linear least squares problem by modification of the Levenberg-Marquardt algorithm. The package includes one R function, which passes all the necessary control parameters to the

Re: [R] inst directory

2004-10-09 Thread Prof Brian Ripley
On Sat, 9 Oct 2004, Bob Wheeler wrote: R CMD check now balks at my inst directory. It contains a single folder doc, but apparently any folder causes the problem. If inst is empty, the project checks OK. This was not a problem before 1.9. I've checked the documentation, but don't see a

Re: [R] conversion of a data.frame of numerics to a data.frame of factors

2004-10-09 Thread Patrick Giraudoux
This work perfectly with a data.frame. Actually my so-called data.frame was an output of t(), and thus not a data.frame (what I realise after trying Prof. Ripley's example). If applied to an output of t() (eg t(df)), the result is quite unxpected (to me): a strange list of 989 elements, each of

[R] Clustering of variable

2004-10-09 Thread clara gonalves
Hi! Does anybody know if there is any way to do clustering of variables? Thank you in advance, Clara. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Chernoff faces

2004-10-09 Thread Shigeru Mase
There is also an code which is more faithful (?) to the original one by S. Aoki although it is in a Japanese HP http://aoki2.si.gunma-u.ac.jp/R/face.html Quality? Judge yourself! Wolfram Fischer wrote: Kenneth == Kenneth Cabrera [EMAIL PROTECTED] writes: Kenneth Hello everybody: Does any one has a

[R] RE: zero-inflated count models (was polr problem solved)

2004-10-09 Thread Peter Flom
John Fox wrote From your description, it seems possible that there are too many zeros for a Poisson or negative-binomial model. Since the focus of your paper is the methodology, you might want to try a zero-inflated Poisson or negative-binomial model. Though I haven't tried them, I'm aware of two

[R] R-2.0.0 and tcltk package

2004-10-09 Thread Jean Eid
there does not seem to be a package tcltk on CRAN for 2.0.0. I have successfully installed the same package for 1.9.1. In essence I require the package for a GUI interface to setwd. All work fine with 1.9.1. This is on a Linux Debian unstable kernel 2.4.20 The build version of R 2.0.0 (issuing

[R] Is it safe? Cochran etc

2004-10-09 Thread Dan Bolser
I have the following contingency table dat - matrix(c(1,506,13714,878702),nr=2) And I want to test if their is an association between events A:{a,not(a)} and B:{b,not(b)} | b | not(b) | +-++ a | 1 | 13714 | +-++ not(a) | 506 |

Re: [R] Is it safe? Cochran etc

2004-10-09 Thread Dan Bolser
Why can't I just use Log odds? Does the standard error of the logs score depend on a similar chisq assumption? On Sat, 9 Oct 2004, Dan Bolser wrote: I have the following contingency table dat - matrix(c(1,506,13714,878702),nr=2) And I want to test if their is an association between events

[R] RSQLite query error

2004-10-09 Thread Mikkel Grum
Dear R-helpers, I ran the following little test on RSQLite and got the data below from the query. Unless I've made some mistake, the results of both the where and order by statements have problems: library(RSQLite) con-dbConnect(dbDriver(SQLite),dbname=test) data(USArrests)

Re: [R] Is it safe? Cochran etc

2004-10-09 Thread Frederico Zanqueta Poleto
Dan, I don't know what is the theory behind this hybrid option and what consists the Cochran conditions. However, I think even if you suppose the asymptotic distribution is not too accurate, because your sampled 1, there is a too strong association of A and B, as this can be noticed by

Re: [R] RSQLite query error

2004-10-09 Thread David James
Mikkel Grum wrote: Dear R-helpers, I ran the following little test on RSQLite and got the data below from the query. Unless I've made some mistake, the results of both the where and order by statements have problems: This is due to the fact that SQLite as of version 2.8 is typeless and

Re: [R] R-2.0.0 and tcltk package

2004-10-09 Thread Prof Brian Ripley
On Sat, 9 Oct 2004, Jean Eid wrote: there does not seem to be a package tcltk on CRAN for 2.0.0. It is part of the R tarball, so library(tcltk) should do something. It has never been on CRAN, to my knowledge. See the R FAQ 5.1.1 (as the posting guide suggests). -- Brian D. Ripley,

Re: [R] Is it safe? Cochran etc

2004-10-09 Thread Kjetil Brinchmann Halvorsen
Dan Bolser wrote: I have the following contingency table dat - matrix(c(1,506,13714,878702),nr=2) And I want to test if their is an association between events A:{a,not(a)} and B:{b,not(b)} | b | not(b) | +-++ a | 1 | 13714 | +-++ not(a) |

Re: [R] Clustering of variable

2004-10-09 Thread Kjetil Brinchmann Halvorsen
clara gonçalves wrote: Hi! Does anybody know if there is any way to do clustering of variables? Thank you in advance, varclust() in Hmisc Kjetil Clara. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

RE: [R] which() and value replacement in a matrix

2004-10-09 Thread Liaw, Andy
Use matrix indexing, instead of generating the index vector: x - as.data.frame(x) x[x 25] - -999 x V1 V2 V3 1 7 14 16 220 248 317 18 11 419 -999 -999 5234 15 6 -999 -9995 7219 12 82223 9 -999 131 106 10

[R] Mac: importing saved PDF figures into Illustrator CS

2004-10-09 Thread Tilo Blenk
Rene Bertin wrote: Hello, This is a usage question for others with experience of R under the Aqua Mac OS X interface. Basically, I don't succeed in importing PDF files (created with the 'Save As' menu to Quartz device windows) into Illustrator. Versions up to (and including?) 10 loose

[R] functions

2004-10-09 Thread Benjamin M. Osborne
Does anyone know of a list of functions that R already knows? I can't seem to find this anywhere in the help documentation. For example, I want to count the number of occurences of a certain value in a column of a data frame: What do I have to do to tell R to Count? Thanks, Ben Osborne --

Re: [R] functions

2004-10-09 Thread Duncan Murdoch
On Sat, 9 Oct 2004 18:12:20 -0400, Benjamin M. Osborne [EMAIL PROTECTED] wrote: Does anyone know of a list of functions that R already knows? I can't seem to find this anywhere in the help documentation. For example, I want to count the number of occurences of a certain value in a column of a

[R] Fatal error: unable to restore saved data in .Rdata

2004-10-09 Thread Susan Yeh
Dear Sir/Madam: I'm currently running R 1.9.1 on a Mac OS 10.3. I have been able to use it fine for a while now. But for some reason, it now crashes each time I open it with the error: Fatal error: unable to restore saved data in .Rdata I read on the online help forum that I can just delete the

[R] How to install a package that needs to see oher pkg dependencies:

2004-10-09 Thread Aldi Kraja
Hi, I am trying to install the genetics package on a server with Linux, Fedora. I installed it in a PC and worked fine. In a server since I am not used with R, I am not sure what do I need to change so genetics pkg can see some package dependencies: Any suggestion is appreciated, Aldi Note:

Re: [R] R-2.0.0 and tcltk package

2004-10-09 Thread Dirk Eddelbuettel
On Sat, Oct 09, 2004 at 12:36:51PM -0400, Jean Eid wrote: there does not seem to be a package tcltk on CRAN for 2.0.0. I have successfully installed the same package for 1.9.1. In essence I require the package for a GUI interface to setwd. All work fine with 1.9.1. This is on a Linux