[R] disaggregating table

2006-10-22 Thread Jeff Miller
Hi all, This should be easy, but I can't seem to figure it out. I have a table like this named newtable a1 a2 a3 a4 Cnts Score 1 100 4 3.28 1 011 2 2.63 I want the following: a1 a2 a3 a4 Cnts Score 1 100 4

Re: [R] Copula fitting

2006-10-22 Thread Uwe Ligges
Puskás László wrote: > Hi, > > Is anybody using Copula package for fitting copulas to own > datas? I always got this message: > > Error in qnorm(p, mean, sd, lower.tail, log.p) : > Non-numeric argument to mathematical function Please specify version of package and R as well as the

[R] setting options in a list

2006-10-22 Thread Michael Toews
Hi, I have a question regarding setting options that use a list. There are none that I'm aware of in R-base, however, it can make sense for some custom situations. For example: options(myoptions=list(lwd=1,density=NULL,angle=45,col="green")) Getting values is simple: getOption("myoptions")$col

[R] Hi,about the png device

2006-10-22 Thread lidaof
Hi, i have a problem with png device that is my error log "> png(file="figure1.png");plot(significant,type="scatter");dev.off() Error in X11(paste("png::", filename, sep = ""), width, height, pointsize, : unable to start device PNG In addition: Warning message: unable to open connection to

[R] Copula fitting

2006-10-22 Thread Puskás László
Hi, Is anybody using Copula package for fitting copulas to own datas? I always got this message: Error in qnorm(p, mean, sd, lower.tail, log.p) : Non-numeric argument to mathematical function thanks, Psl ___ Ossza meg barát

[R] R tcl/tk

2006-10-22 Thread Charles Annis, P.E.
Friends: I am a long-time R user, learning tcl/tk, and am tying myself in knots over something that should be simple. I want to create a frame and put that frame inside the toplevel frame. This works (i.e. it places text in col 1, and the corresponding entry box in column 2) and later frame.2

Re: [R] Key colour problem with lattice plot

2006-10-22 Thread Geoff Russell
Thanks Deepayan, On 10/23/06, Deepayan Sarkar <[EMAIL PROTECTED]> wrote: > On 10/22/06, Geoff Russell <[EMAIL PROTECTED]> wrote: > > Many thanks Deepayan, > > > > .. > > I'm not sure what you mean by that, but my intention in pointing you > towards show.settings() was to make you realize that you

Re: [R] Key colour problem with lattice plot

2006-10-22 Thread Deepayan Sarkar
On 10/22/06, Geoff Russell <[EMAIL PROTECTED]> wrote: > Many thanks Deepayan, > > On 10/23/06, Deepayan Sarkar <[EMAIL PROTECTED]> wrote: > > On 10/22/06, Geoff Russell <[EMAIL PROTECTED]> wrote: > > > Hi, > > > > > > I seem to have a key colour problem with a lattice barchart. The colours > > > o

Re: [R] Multilevel model ("lme") question

2006-10-22 Thread Lukas Rode
Thanks Douglas for the detailed explanation; it was very helpful (I'll definitely give lme4 a try as well). For the record: I also found it helpful to play with model.matrix (e.g. p. 14ff and p.28 of the book by Jose Pinheiro and Douglas) to see how the formulas for fixed and random effects are rel

Re: [R] Key colour problem with lattice plot

2006-10-22 Thread Geoff Russell
Many thanks Deepayan, On 10/23/06, Deepayan Sarkar <[EMAIL PROTECTED]> wrote: > On 10/22/06, Geoff Russell <[EMAIL PROTECTED]> wrote: > > Hi, > > > > I seem to have a key colour problem with a lattice barchart. The colours on > > my > > key rectangle don't match the colours on the barchart. > > T

Re: [R] Key colour problem with lattice plot

2006-10-22 Thread Deepayan Sarkar
On 10/22/06, Geoff Russell <[EMAIL PROTECTED]> wrote: > Hi, > > I seem to have a key colour problem with a lattice barchart. The colours on my > key rectangle don't match the colours on the barchart. That's because you are using the wrong colors in the key (see the output of show.settings()). I w

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Marc Schwartz
On Sun, 2006-10-22 at 14:03 -0500, Marc Schwartz wrote: > On Sun, 2006-10-22 at 14:37 -0400, Wensui Liu wrote: > > Thank you so much, Marc and Peter, > > > > Your method works great if I want to convert N dummies into N-level > > factor. But what if I want to convert N dummies into (N+1)-level > >

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Wensui Liu
Thank you so much, Marc and Peter, Your method works great if I want to convert N dummies into N-level factor. But what if I want to convert N dummies into (N+1)-level factor? I tried both ways but none works. Again, thank you so much! On 10/22/06, Marc Schwartz <[EMAIL PROTECTED]> wrote: > On

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Wensui Liu
Great! It seems everyone is having fun with R in the weekend afternoon. Thank you so much, Marc and Peter. On 10/22/06, Marc Schwartz <[EMAIL PROTECTED]> wrote: > On Sun, 2006-10-22 at 14:03 -0500, Marc Schwartz wrote: > > On Sun, 2006-10-22 at 14:37 -0400, Wensui Liu wrote: > > > Thank you so

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Peter Dalgaard
"Wensui Liu" <[EMAIL PROTECTED]> writes: > Thank you so much, Marc and Peter, > > Your method works great if I want to convert N dummies into N-level > factor. But what if I want to convert N dummies into (N+1)-level > factor? I tried both ways but none works. Mine should work, except for the

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Marc Schwartz
On Sun, 2006-10-22 at 14:37 -0400, Wensui Liu wrote: > Thank you so much, Marc and Peter, > > Your method works great if I want to convert N dummies into N-level > factor. But what if I want to convert N dummies into (N+1)-level > factor? I tried both ways but none works. > > Again, thank you so

Re: [R] problem with mode of marginal distriubtion of rdiri chlet{gtools}

2006-10-22 Thread Rob Carnell
hong qin gmail.com> writes: > > Hi all, > > I have a problem using rdirichlet{gtools}. > For Dir( a1, a2, ..., a_n), its mode can be found at $( a_i -1)/ ( > \sum_{i}a_i - n)$; > The means are $a_i / (\sum_{i} a_i ) $; > > I tried to study the above properties using rdirichlet from gtools.

Re: [R] Book for Maximum Likelihood Methods in R

2006-10-22 Thread Ben Bolker
Milton Cezar Ribeiro yahoo.com.br> writes: > > Hi R-guys, > > Is there a good book for maximum likelihood methods in R? > I´m looking mainly to Ecological approaches to apply ML. > ** BLATANT PLUG warning **: I'm working on a book on statistical/ecological modeling in R, to be published by

Re: [R] Possible bugs in 'seek' and 'readBin'

2006-10-22 Thread Zepu Zhang
Sorry I can't provide reproducible code because it involves data files. I saw no warning messages. I'm using Mac 10.3, but the R is compiled by myself from command line. I'm reading in data from a XDR formatted file. The first 4 bytes is a 'long' integer, the remainder is double. --- problem wit

Re: [R] how to create a vector with different categories in a simple way?

2006-10-22 Thread Marc Schwartz
On Sun, 2006-10-22 at 17:38 +0200, Jenny persson wrote: > Hi R-users, > > I have a matrice called layout which contains 5 columns:id, name, > row, column and block. The column called "block" has totally 48 blocks > and looks like > > > 1 2 3 4 > 5 6 7 8 > 9 10 11 12 > 1

Re: [R] Recipe: Print a Color Sampler

2006-10-22 Thread Gabor Grothendieck
Plotting points with a pch = NA also plots nopoints so we could remove the type= and the switch statement. Also it needs an on.exit to reset par back to what it was. The cex seems to have no effect so I removed it. getColorName <- function(colorNumber) colors()[colorNumber] # pch = NA means no

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Marc Schwartz
On Sun, 2006-10-22 at 09:35 +0200, Peter Dalgaard wrote: > Marc Schwartz <[EMAIL PROTECTED]> writes: > > > On Sat, 2006-10-21 at 21:04 -0400, Wensui Liu wrote: > > > Dear Listers, > > > > > > I am wondering how to convert multiple dummy variables to 1 factor > > > variable. > > > > > > Thanks.

Re: [R] Recipe: Print a Color Sampler

2006-10-22 Thread Ana Nelson
Thanks, guys, for all the mods! :-) Here is another version to allow people to display the dots if they wish and also to set a different background color. I think this is useful just for identifying the colors, but also to see what a particular shape/color/background color combination will l

Re: [R] object(s) are masked from package - what does it mean?

2006-10-22 Thread tom soyer
Thanks Peter. That make sense. On 22 Oct 2006 09:26:27 +0200, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > > "tom soyer" <[EMAIL PROTECTED]> writes: > > > Thanks Marc and Gregor for the detailed explanation. You are right, the > > masking is potentially dangerous. Since R is object oriented, I am >

Re: [R] Suitability of R for Algorithm simulations

2006-10-22 Thread François Pinard
Hi, people. A correspondent puts me in front of a reply I sent to r-help, a few weeks ago, and quoted below. I should have been tired when I sent it. Please replace "Eiffel" by "Erlang" all over. Sorry for this error. Date: 2006-10-05 00:43:36 Message-ID: [EMAIL PROTECTED] [Ethan B. Fin

[R] how to create a vector with different categories in a simple way?

2006-10-22 Thread Jenny persson
Hi R-users, I have a matrice called layout which contains 5 columns:id, name, row, column and block. The column called "block" has totally 48 blocks and looks like 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 3

Re: [R] Multilevel model ("lme") question

2006-10-22 Thread Douglas Bates
On 10/21/06, Lukas Rode <[EMAIL PROTECTED]> wrote: > Dear list, > > I'm trying to fit a multilevel (mixed-effects) model using the lme function > (package nlme) in R 2.4.0. As a mixed-effects newbie I'm neither sure about > the modeling nor the correct R syntax. You may also want to consider using

Re: [R] Recipe: Print a Color Sampler

2006-10-22 Thread Antonio, Fabio Di Narzo
Really fine. Should this go on the wiki? Antonio. 2006/10/22, Gabor Grothendieck <[EMAIL PROTECTED]>: > Here it is again corrected: > > getColorName <- function(colorNumber) colors()[colorNumber] > printColorSampler <- function(n = 0) { >i <- seq(colors()) >k <- ceiling(sqrt(length(i))) >

[R] Book for Maximum Likelihood Methods in R

2006-10-22 Thread Milton Cezar Ribeiro
Hi R-guys, Is there a good book for maximum likelihood methods in R? I´m looking mainly to Ecological approaches to apply ML. Any help? Kind regards, miltinho BRAZIL - Você quer respostas para suas perguntas? Ou você sabe muito e quer com

Re: [R] Recipe: Print a Color Sampler

2006-10-22 Thread Gabor Grothendieck
Here it is again corrected: getColorName <- function(colorNumber) colors()[colorNumber] printColorSampler <- function(n = 0) { i <- seq(colors()) k <- ceiling(sqrt(length(i))) xy <- cbind(floor(i/k)*2, i %% k) plot(xy, type = "n", axes = FALSE, xlab = "", ylab = "") text(xy[,1]+.5,

Re: [R] Recipe: Print a Color Sampler

2006-10-22 Thread Gabor Grothendieck
I have removed the dots, vectorized it and changed the argument to the number points to be identified (default 0): getColorName <- function(colorNumber) colors()[colorNumber] printColorSampler <- function(n = 0) { i <- seq(colors()) k <- ceiling(sqrt(length(i))) xy <- cbind(floor(i/k)*2,

Re: [R] Recipe: Print a Color Sampler

2006-10-22 Thread Stefano Calza
Thanks, might be very useful. I attached a slightly modified version that uses identify (whith an argument to choose if use it), to return the name of the colors on the panel. Left-click on the choosen colors, and right click to end. Stefano On Sun, Oct 22, 2006 at 12:51:19PM +0100, Ana Nelso

[R] Recipe: Print a Color Sampler

2006-10-22 Thread Ana Nelson
I wrote this to help me choose a suitable plot color. You can set the value of pch to whichever type of point you are using in your plot, and then you can see what the various colors will look like. When this plot it stretched to A4/Letter landscape it prints quite nicely, even though it wi

[R] Key colour problem with lattice plot

2006-10-22 Thread Geoff Russell
Hi, I seem to have a key colour problem with a lattice barchart. The colours on my key rectangle don't match the colours on the barchart. Here is my data frame: "LandFill" "Ruminants" "United States (USA)" .21428 5528.16 "France" 200.527083 1299.87 "Australia" 185.878368 2448.17 "Russian Fed

Re: [R] Question:shardsplot (package:klaR)

2006-10-22 Thread Uwe Ligges
I am forwarding this question to Nils Raabe, the author of shardsplot(). Nils? Uwe Ligges ebi wrote: > Dear all, > > > I have a question on the shardsplot package:klaR(see the below Example). > > Plese tell me the meanings of " logstand <- t((t(logcount) / sdlogcount) * > c(1,2,6,5,5,3))",

Re: [R] how to convert multiple dummy variables to 1 factor variable?

2006-10-22 Thread Peter Dalgaard
Marc Schwartz <[EMAIL PROTECTED]> writes: > On Sat, 2006-10-21 at 21:04 -0400, Wensui Liu wrote: > > Dear Listers, > > > > I am wondering how to convert multiple dummy variables to 1 factor variable. > > > > Thanks. > > > > wensui > > I was thinking of a function that is essentially the revers

Re: [R] object(s) are masked from package - what does it mean?

2006-10-22 Thread Peter Dalgaard
"tom soyer" <[EMAIL PROTECTED]> writes: > Thanks Marc and Gregor for the detailed explanation. You are right, the > masking is potentially dangerous. Since R is object oriented, I am surprised > that this is an issue. Does this mean that encapsulation does not exist in > R? R does have namespaces

Re: [R] Problem with installing Hmisc and Design: gfrotran: error

2006-10-22 Thread Prof Brian Ripley
On Sat, 21 Oct 2006, Stephan Lindner wrote: > cannot exec 'f951': No such file or directory > Reply-To: > > Hi all, > > I'm using R on an Suse Linux system. Since Hmisc and Design need both > fortran, I installed gfortran (through rpm). However, I still get an > error message, namely: > > gfrotran

Re: [R] Possible bugs in 'seek' and 'readBin'

2006-10-22 Thread Prof Brian Ripley
Please see the message at the bottom of this posting, and supply the information requested, including your OS. (If the latter is Windows, did you notice the warnings in ?seek.) On Sat, 21 Oct 2006, Zepu Zhang wrote: > I found that > seek(..., origin = 'current', ...) > and > readBin(..., wha