Re: [R] slightly extended lm class

2007-01-01 Thread Prof Brian Ripley
Have you looked at https://stat.ethz.ch/pipermail/r-devel/2006-December/044111.html ? That seems the appropriate mailing list if you want to continue this discussion. On Mon, 1 Jan 2007, ivo welch wrote: > obviously not. any other suggestion? of course, it would be even > better if someone wh

Re: [R] zero random effect sizes with binomial lmer

2007-01-01 Thread Ken Beath
On Sun, 31 Dec 2006, at 05:50 PM, Daniel Ezra Johnson <[EMAIL PROTECTED]> wrote: > > Gregor, > > Thanks for your replies. > > 1) Yes, I have tweaked the data to show as clearly as I can that > this is a > bug, that a tiny change in initial conditions causes the collapse of a > reasonable 'para

Re: [R] slightly extended lm class

2007-01-01 Thread ivo welch
obviously not. any other suggestion? of course, it would be even better if someone who knew what he is doing were to add this functionality on demand (e.g., optional parameters) to the standard lm classes. just a thought... On 1/1/07, Michael Kubovy <[EMAIL PROTECTED]> wrote: > On Jan 1, 2007,

Re: [R] zero random effect sizes with binomial lmer

2007-01-01 Thread Andrew Robinson
Good observation, Daniel. I had not picked up the different default fitting method. Cheers Andrew On Mon, Jan 01, 2007 at 08:59:55PM -0500, Daniel Ezra Johnson wrote: > > From: Andrew Robinson > > Date: Mon 01 Jan 2007 - 19:19:29 GMT > > > > I tried an earlier version of R, on a different pla

Re: [R] slightly extended lm class

2007-01-01 Thread Michael Kubovy
On Jan 1, 2007, at 9:00 PM, ivo welch wrote: > I have written a short lme.R function, which adds normalized > coefficients and White heteroskedasticity-adjusted statistics to the > standard output. Otherwise, it behaves like lm. Is it a good idea to call it lme, since there's a widely used lme()

[R] slightly extended lm class

2007-01-01 Thread ivo welch
Dear R readers: I have written a short lme.R function, which adds normalized coefficients and White heteroskedasticity-adjusted statistics to the standard output. Otherwise, it behaves like lm. This is of course trivial for experts, but for me and other amateur users perhaps helpful. y= rnorm

Re: [R] zero random effect sizes with binomial lmer

2007-01-01 Thread Daniel Ezra Johnson
> From: Andrew Robinson > Date: Mon 01 Jan 2007 - 19:19:29 GMT > > I tried an earlier version of R, on a different platform, and got quite > different results. Sadly, the *earlier* results are the ones that make > sense. Andrew, I tried installing R 2.3.1, it seemed to work as you said, but I

Re: [R] Does SQL group by have a heavy duty equivalent in R

2007-01-01 Thread Prof Brian Ripley
On Mon, 1 Jan 2007, BBands wrote: > On 12/30/06, Farrel Buchinsky <[EMAIL PROTECTED]> wrote: >> Instead I used RODBC >> >> sqlSave(channel,RawSeq) >> to push the table into a Microsoft Access database >> Then a sql query, courtesy of the Microsoft Access Query Wizard a la design >> mode. > > If SQ

Re: [R] Does SQL group by have a heavy duty equivalent in R

2007-01-01 Thread Farrel Buchinsky
Absolutely wonderful. Yes indeed, instead of having two variables being listed as group by in a sql statement one can use a combination of cbind, unique, duplicated and indexing to get the same data without having to go through second programs. [[alternative HTML version deleted]] ___

Re: [R] rimage package broken with fedora upgrade

2007-01-01 Thread John Kornak
Thanks to Jonathan and Oleg for their suggestions/comments. Their responses are added at the end of this email. I tried Oleg's suggestion of giving me write permissions on /usr/lib/R/library. This didn't seem to work (I tried various combinations of changing permissions and ownership but nothing

Re: [R] matrix size

2007-01-01 Thread roger koenker
On Jan 1, 2007, at 4:43 PM, Armelini, Guillermo wrote: > Hello everyone > Could anybody tell me how to set the following matrix? > > n2<-matrix(nrow=10185,ncol=10185,seq(0,0,length=103734225)) You can use: library(SparseM) as.matrix.coo(0,10185,10185) but then you need to find something intere

Re: [R] matrix size

2007-01-01 Thread Martin Maechler
> "GArmelini" == Armelini, Guillermo <[EMAIL PROTECTED]> > on Mon, 1 Jan 2007 23:43:31 +0100 writes: GArmelini> Hello everyone Could anybody tell me how to set GArmelini> the following matrix? GArmelini> n2<-matrix(nrow=10185,ncol=10185,seq(0,0,length=103734225)) GA

[R] matrix size

2007-01-01 Thread Armelini, Guillermo
Hello everyone Could anybody tell me how to set the following matrix? n2<-matrix(nrow=10185,ncol=10185,seq(0,0,length=103734225)) R answer was Error: cannot allocate vector of size 810423 Kb Are there any solution? I tried to increase the memory size but it didn't work G This message has b

Re: [R] advice on semi-serious attempt to extend summary

2007-01-01 Thread Gabor Grothendieck
Try: methods(print) stats:::print.summary.lm On 1/1/07, ivo welch <[EMAIL PROTECTED]> wrote: > Dear R wizards: > > I am trying (finally) to build a function that might be useful to > others. In particular, I want to create a summary.lme (extended lm) > method that [a] adds normalized coefficient

[R] advice on semi-serious attempt to extend summary

2007-01-01 Thread ivo welch
Dear R wizards: I am trying (finally) to build a function that might be useful to others. In particular, I want to create a summary.lme (extended lm) method that [a] adds normalized coefficients and [b] white heteroskedasticity adjusted se's and T's. I believe I already know how to do the progra

Re: [R] Any container in R?

2007-01-01 Thread Thomas Lumley
On Mon, 1 Jan 2007, Feng Qiu wrote: > Hi Duncan: > Thank you very much! I checked out unique(), it does exactly what I > want. But I'm still curious about if R provides "STL"(standard template > library). No. Some things the STL does aren't needed in R, others are implemented differentl

Re: [R] Help with filled.contour()

2007-01-01 Thread Michael Kubovy
On Jan 1, 2007, at 2:13 PM, Dieter Menne wrote: > Michael Kubovy virginia.edu> writes: > >> >> The following plot is a first approximation to what I need: > .. see below > >> *** >> It has four drawbacks. I would like: >> (1) to remove the white grid; >> (2) to rem

Re: [R] zero random effect sizes with binomial lmer

2007-01-01 Thread Andrew Robinson
According to pp 197-198 of MASS 4, the Hauck-Donner phenomenon refers to cases when the Wald approximations and the likelihood ratio tests have different p values because of the former underestimating the the change in log-likelihood on setting \beta_i = 0. This seems quite different to me than th

Re: [R] Types of objects?

2007-01-01 Thread Gabor Grothendieck
Try class, typeof and mode. On 1/1/07, Raymond Balise <[EMAIL PROTECTED]> wrote: > What I am looking for is a function that goes with section 2.8 of the > "Introduction to R" PDF that ships with R. It says that the R objects are > vectors, matrices, factors, lists, data frames and functions. Is

Re: [R] Help with filled.contour()

2007-01-01 Thread Dieter Menne
Michael Kubovy virginia.edu> writes: > > The following plot is a first approximation to what I need: .. see below > *** > It has four drawbacks. I would like: > (1) to remove the white grid; > (2) to remove the white contours; > (3) its appearance to be smooth;

Re: [R] Types of objects?

2007-01-01 Thread Raymond Balise
What I am looking for is a function that goes with section 2.8 of the "Introduction to R" PDF that ships with R. It says that the R objects are vectors, matrices, factors, lists, data frames and functions. Is there a function that returns this level of information or is that taxonomy just somethi

Re: [R] Problems with R

2007-01-01 Thread Feng Qiu
You got to supply the absolute path of your files and the seperating symbol in your csv file. Please read the help file by typing help(read.csv), you will have no problem. - Original Message - From: "Obinna Duru" <[EMAIL PROTECTED]> To: Sent: Monday, January 01, 2007 4:09 AM Subject:

Re: [R] Any container in R?

2007-01-01 Thread Feng Qiu
Hi Duncan: Thank you very much! I checked out unique(), it does exactly what I want. But I'm still curious about if R provides "STL"(standard template library). Best, Feng - Original Message - From: "Duncan Murdoch" <[EMAIL PROTECTED]> To: "Feng Qiu" <[EMAIL PROTECTED]> Cc:

Re: [R] Does SQL group by have a heavy duty equivalent in R

2007-01-01 Thread BBands
On 12/30/06, Farrel Buchinsky <[EMAIL PROTECTED]> wrote: > Instead I used RODBC > > sqlSave(channel,RawSeq) > to push the table into a Microsoft Access database > Then a sql query, courtesy of the Microsoft Access Query Wizard a la design > mode. If SQL does prove to be part of your approach you m

Re: [R] zero random effect sizes with binomial lmer

2007-01-01 Thread Dieter Menne
Daniel Ezra Johnson babel.ling.upenn.edu> writes: > > I am fitting models to the responses to a questionnaire that has > seven yes/no questions (Item). For each combination of Subject and > Item, the variable Response is coded as 0 or 1. > > I want to include random effects for both Subject and

Re: [R] Types of objects?

2007-01-01 Thread Dieter Menne
Raymond Balise gmail.com> writes: > > Is there a list of object types used by R? Is there a function that I can > use to tell the type of an object? typeof() Dieter __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

Re: [R] Any container in R?

2007-01-01 Thread Duncan Murdoch
On 1/1/2007 11:22 AM, Feng Qiu wrote: > Hi Duncan: > Thanks for your hints. > I'm trying to collect distinct elements in one column in a matrix. > If there is a map, I can easily build up such a collection. While if using > list, I have to check by myself if this element already

[R] Types of objects?

2007-01-01 Thread Raymond Balise
Is there a list of object types used by R? Is there a function that I can use to tell the type of an object? Thanks a bunch, Ray [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinf

[R] kernel density estimation with many variables (50 variables)

2007-01-01 Thread Serge Zaugg
Hello, Does any one know an R-function that performs kernel density estimation when there are many variables ? My data has up to 50 variables and several thousand data points. Does any one know an alternative method to perform nonparametric density estimation on high dimensional data ?

[R] zero random effect sizes with binomial lmer [sorry, ignore previous]

2007-01-01 Thread Daniel Ezra Johnson
I am fitting models to the responses to a questionnaire that has seven yes/no questions (Item). For each combination of Subject and Item, the variable Response is coded as 0 or 1. I want to include random effects for both Subject and Item. While I understand that the datasets are fairly small, and

[R] arima.sim with a periodic model

2007-01-01 Thread Philip
Hi all. I have a periodiv arma model and I want to simulate it. In S-plus, the following works for me: phi <- 0.9 theta <- 0 p <- 1 # period model <- list(ar=phi, ma=theta, period=p) Yt <- arima.sim(model, n=250) How do I do something like this "period=12" in R? I read help(arima.sim) but it d

[R] autofill for a tkentry field

2007-01-01 Thread Adrian Dragulescu
Hello all, I'm experimenting with tlk/tk and I find it very useful. I would like to have the text of a tkentry be automatically filled when the user starts typing, according to a predetermined list. I've found a solution on the tcl wiki pages, but I don't quite know how to translate it in R...

Re: [R] Any container in R?

2007-01-01 Thread Feng Qiu
Hi Duncan: Thanks for your hints. I'm trying to collect distinct elements in one column in a matrix. If there is a map, I can easily build up such a collection. While if using list, I have to check by myself if this element already exists in the collection every time I examine

Re: [R] [R-sig-Geo] plot methods in sp

2007-01-01 Thread Edzer J. Pebesma
Patrick, I can see this is confusing. Besides methods(plot), try after library(sp) a showMethods(plot) to get an overview of the S4-style plot methods; methods(plot) only shows the S3-style plot methods. Then, class?SpatialPolygons gives methods available for this class, among which plot, and

Re: [R] Any container in R?

2007-01-01 Thread Duncan Murdoch
On 1/1/2007 10:17 AM, Feng Qiu wrote: > R has list and array to contain elements. But does R have more powerful > container, such as "map" as in C++ STL? or is there such a package? In what way are maps more powerful than lists? You can use names to index lists. The other container in R is th

[R] Any container in R?

2007-01-01 Thread Feng Qiu
R has list and array to contain elements. But does R have more powerful container, such as "map" as in C++ STL? or is there such a package? Thanks and Happy 2007! Best, Feng __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listi

[R] Help with filled.contour()

2007-01-01 Thread Michael Kubovy
The following plot is a first approximation to what I need: *** mu1 <- 0 mu2 <- 5 s <- 1 x <- seq(-2.5, 7.5, length = 41) y <- seq(-2.5, 2.5, length = 41) f <- function(x,y){ term1 <- 1/(2*pi*sqrt(s*s)) term2 <- -1/2 term3 <- (x - mu1)^2/s term4 <- (y - mu1)^2/s ter

Re: [R] Problems with R

2007-01-01 Thread Duncan Murdoch
On 1/1/2007 4:09 AM, Obinna Duru wrote: > Hello, > > Please I have problems with R. I downloaded R from UC Berkeley CRAN > mirror, but each time I try to read a data file using the command > 'moi.data <- read.csv(choose.files())' or 'my.data <- choose.files() > #' , the program becomes unstable

Re: [R] New to R

2007-01-01 Thread Michael Kubovy
On Dec 31, 2006, at 10:15 PM, Obinna Duru wrote: > I have installed acepack but efforts to get started has been > unsuccessful. I can't seem to be able to load my data files because I > am yet to figure the syntax to use. Is there a work directory in R > where I can put my files and call them anyt

Re: [R] zero random effect sizes with binomial lmer

2007-01-01 Thread Gregor Gorjanc
Daniel Ezra Johnson babel.ling.upenn.edu> writes: > 1) Yes, I have tweaked the data to show as clearly as I can that this is a > bug, that a tiny change in initial conditions causes the collapse of a > reasonable 'parameter' estimate. I would not call this a bug, since this is related to data a

[R] Problems with R

2007-01-01 Thread Obinna Duru
Hello, Please I have problems with R. I downloaded R from UC Berkeley CRAN mirror, but each time I try to read a data file using the command 'moi.data <- read.csv(choose.files())' or 'my.data <- choose.files() #' , the program becomes unstable and exits itself. However, I just noticed from my