Re: [R] Problem installing packages in R 2.4.1

2007-03-21 Thread Uwe Ligges
Paul Gottemoller wrote: > I am attempting to install a anchoring vingettes package for R and I seem to > have a problem with unzipping the files in R. I use the command that can be > found here . > > Here is the dialog I get after entering the command: > >>

Re: [R] problem with RCurl install on Unix

2007-03-21 Thread Prof Brian Ripley
What is 'Unix' here? It seems you mean ppc64 Linux (which is not Unix). On Wed, 21 Mar 2007, Mark W Kimpel wrote: I am having trouble getting an install of RCurl to work properly on a Unix server. The steps I have taken are: 1. installed cUrl from source without difficulty 2. installed RCurl

Re: [R] dynamic linear models in R

2007-03-21 Thread Gavin Simpson
On Wed, 2007-03-21 at 23:58 -0400, Johann Hibschman wrote: > Hi all, > > I've just started working my way through Mike West and Jeff Harrison's > _Bayesian Forecasting and Dynamic Models_, and I was wondering if > there were any publically-available packages to handle dynamic linear > models, as t

[R] Arrellano Bond and Blundell and Bond estimators in R

2007-03-21 Thread justin bem
Hi all, Is there a package to estimate Arrellano Bond estimator and blundelll and bond system estimator for dynamic panel in R ? Best regards. Justin BEM Elève Ingénieur Statisticien Economiste BP 294 Yaoundé. Tél (00237)9597295. __ R-help@stat.math.e

Re: [R] R difftime function: How can we fix the difftime unit?

2007-03-21 Thread Marc Schwartz
On Thu, 2007-03-22 at 11:33 +0530, d. sarthi maheshwari wrote: > Hi, > > I am trying to take difference of two time objects. I want to fix the > result's unit to minutes. How can I do that? > > Here is an example: > > > difftime(x, y) > Time difference of 2.030720 hours > > difftime(x, z) > Time

Re: [R] R difftime function: How can we fix the difftime unit?

2007-03-21 Thread Frede Aakmann Tøgersen
That's easy: difftime(x, y, units = "min") Best regards Frede Fra: [EMAIL PROTECTED] på vegne af d. sarthi maheshwari Sendt: to 22-03-2007 07:03 Til: r-help@stat.math.ethz.ch Emne: [R] R difftime function: How can we fix the difftime unit? Hi, I am tryin

[R] R difftime function: How can we fix the difftime unit?

2007-03-21 Thread d. sarthi maheshwari
Hi, I am trying to take difference of two time objects. I want to fix the result's unit to minutes. How can I do that? Here is an example: > difftime(x, y) Time difference of 2.030720 hours > difftime(x, z) Time difference of 30.34672 mins where x = '2007-03-05 08:32:58' y = '2007-03-

Re: [R] export table

2007-03-21 Thread Stephen Tucker
You can also use write.table(y,"D:/Desktop/export_table.txt", col.names=NA) with this function call, the row names get printed but the column names are offset so that they are aligned in the right column. --- Chuck Cleland <[EMAIL PROTECTED]> wrote: > Sergio Della Franca wrote: > > Dear R-He

Re: [R] export table

2007-03-21 Thread Stephen Tucker
You can also use write.table(y,"D:/Desktop/export_table.txt", col.names=NA) with this function call, the row names get printed but the column names are offset so that they are aligned in the right column. --- Chuck Cleland <[EMAIL PROTECTED]> wrote: > Sergio Della Franca wrote: > > Dear R-He

Re: [R] Ticks on barplots

2007-03-21 Thread Stephen Tucker
Hi Mike, you can try using axTicks as in this example (you can also use pretty instead). Below, instead of a barplot() I have used plot() but with type="h" and lend=3 (see ?par for details on lend) which I have found to be useful at times when making barchart-like plots with time-series. In any ca

[R] DBI + RSQLite or SQLiteDF?

2007-03-21 Thread Bing Ho
Hello, I've finally reached the wall - my 2gb RAM machine simply can't handle the datasets that I am working with. I've tried a 64-bit compile of R on a 8gb RAM machine but that's not always available to use. Now there are several proposed ways around this, but it seems the most general solution

[R] dynamic linear models in R

2007-03-21 Thread Johann Hibschman
Hi all, I've just started working my way through Mike West and Jeff Harrison's _Bayesian Forecasting and Dynamic Models_, and I was wondering if there were any publically-available packages to handle dynamic linear models, as they describe. I found the "dynlm" package, but either I don't yet unde

Re: [R] problem with RCurl install on Unix

2007-03-21 Thread Steven McKinney
I get the same problem, and haven't figured it out yet. Is it a 32bit/64bit clash? (Similarly, I don't have RMySQL up and running cleanly.) > library(RCurl) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library '/share/apps/R/R-2.4.1/library/RCurl/

Re: [R] Labelling a second y-axis

2007-03-21 Thread Marc Schwartz
On Wed, 2007-03-21 at 18:14 -0700, Jesse Sinclair wrote: > Hi, > I am using the following code as an example function to create > multiple y-axes on one plot. I have it working fine however, I can't > seem to add a label on the second (right) axis. I have tried adding > ylab="y2" in the axis

[R] Labelling a second y-axis

2007-03-21 Thread Jesse Sinclair
Hi, I am using the following code as an example function to create multiple y-axes on one plot. I have it working fine however, I can't seem to add a label on the second (right) axis. I have tried adding ylab="y2" in the axis call but, that didn't work; any ideas? Thanks, Jesse Code: funct

[R] problem with RCurl install on Unix

2007-03-21 Thread Mark W Kimpel
I am having trouble getting an install of RCurl to work properly on a Unix server. The steps I have taken are: 1. installed cUrl from source without difficulty 2. installed RCurl from source using the command ~/R_HOME/R-devel/bin/R CMD INSTALL -l ~/R_HOME/R-devel/library ~/RCurl_0.8-0.ta

[R] Clarification

2007-03-21 Thread Sender
Hello, Could someone help me understand why one of these "structures" would be preferred over the other? x <- setClass( "patient", representation( data = "numeric", fname = "character", lname="character", disease = "character" )) jd <- new("patient", data = rn

[R] Snow Package and SPRNG: Will it solve my problem?

2007-03-21 Thread robert . robinson
Hello and thanks in advance for your time. I currently have a simulation running on my cluster with the help of snow that relies on global variables being changes regularly to random values. It uses these values, lets call them x1 x2 and x3, in custom functions for logliklyhood and score

Re: [R] Stepwise Logistic Regression

2007-03-21 Thread Cody_Hamilton
Rich Ulrich has compiled some posts (I believe from the S list) on some of the dangers of stepwise regression: http://www.pitt.edu/~wpilib/statfaq/regrfaq.html Regards, -Cody "Sergio Della

[R] Problem installing packages in R 2.4.1

2007-03-21 Thread Paul Gottemoller
I am attempting to install a anchoring vingettes package for R and I seem to have a problem with unzipping the files in R. I use the command that can be found here . Here is the dialog I get after entering the command: > install.packages("anchors", dependen

Re: [R] Integer Optimization

2007-03-21 Thread apjaworski
Dimitri, I am assuming you mean integer linear optimization problem. Try lpSolve. I think you will find it much easier to use. Cheers, Andy __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail: [EMAIL PROTECTED] Tel: (651)

Re: [R] Ticks on barplots

2007-03-21 Thread Marc Schwartz
On Wed, 2007-03-21 at 14:40 -0400, Mike Prager wrote: > Marc Schwartz <[EMAIL PROTECTED]> wrote: > > > On Tue, 2007-03-20 at 18:04 -0400, Michael H. Prager wrote: > > > I am generating stacked barplots of age-composition of fish populations > > > (Y) over time (X). As there are many years, not e

Re: [R] RMySQL *was* working...

2007-03-21 Thread Seth Falcon
Pete Cap <[EMAIL PROTECTED]> writes: > List, > > Last week with the help of Uwe and some other folks I was able to get > RMySQL 0.5-7 compiled against R 2.4.1 and MySQL 5.0.27. It was > working fine--I was able to send select queries to the db, put the > results in a data frame, and so forth. > >

Re: [R] Implementing trees in R

2007-03-21 Thread Bálint Czúcz
Dear Kevin, for an example how this idea is implemented in a working package, have a look at the help page of BinaryTree class in package party. The @tree node of such an object is a recursive list of this kind. HTH, Bálint On 3/16/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > Let me reph

Re: [R] Stepwise Logistic Regression

2007-03-21 Thread Frank E Harrell Jr
Sergio Della Franca wrote: > Dear R-Helpers, > > I'd like to perform a Logistic Regression whit a Stepwise method. > > > Can you tell me how can i proceed to develop this procedure? > > > Thank you in advance. > > > Sergio Della Franca. If the number of events is not incredibly large, you c

Re: [R] Stepwise Logistic Regression

2007-03-21 Thread David Farrar
Stella, An obvious answer is ?step However, I'm having a bit of a problem with it, lately. I got it wording with backwards selection, then it didn't work when I changed the "direction" from "backward" to "both" (backwards and forwards). I would like to know whatever you find that work

Re: [R] Ticks on barplots

2007-03-21 Thread Mike Prager
Marc Schwartz <[EMAIL PROTECTED]> wrote: > On Tue, 2007-03-20 at 18:04 -0400, Michael H. Prager wrote: > > I am generating stacked barplots of age-composition of fish populations > > (Y) over time (X). As there are many years, not every bars is labeled. > > When looking at the plot, it becomes

[R] Stepwise Logistic Regression

2007-03-21 Thread Sergio Della Franca
Dear R-Helpers, I'd like to perform a Logistic Regression whit a Stepwise method. Can you tell me how can i proceed to develop this procedure? Thank you in advance. Sergio Della Franca. [[alternative HTML version deleted]] __ R-help@stat.

Re: [R] Gaussian Adaptive Quadrature

2007-03-21 Thread Earl F. Glynn
"Ravi Varadhan" <[EMAIL PROTECTED]> wrote in message news:[EMAIL PROTECTED] > The terminology "Gaussian quadrature" is not restricted to Gauss-Hermite > quadrature (where exp(-x^2) is the weight function), but applies more > broadly to Gauss-Legendre, Gauss-Laguerre, etc., where the abscissa are >

Re: [R] export table

2007-03-21 Thread Chuck Cleland
Sergio Della Franca wrote: > Dear R-Helpers, > > I have a problem. > > I want to export from R to .txt my data set(y): > > YEARS PRODUCTS > 1990 10 > 1995 15 > 1997 26 > 1998 29 > 2000 34 > > > I used this code: > > write.table(y,"D:/Desktop/export_table.txt"). > > This p

[R] export table

2007-03-21 Thread Sergio Della Franca
Dear R-Helpers, I have a problem. I want to export from R to .txt my data set(y): YEARS PRODUCTS 1990 10 1995 15 1997 26 1998 29 2000 34 I used this code: write.table(y,"D:/Desktop/export_table.txt"). This procedure run correctly, but i acquired this result: YEARS PRODU

Re: [R] question on suppressing error messages with Rmath library

2007-03-21 Thread Ranjan Maitra
Hi Luke, Thanks! Sorry, my error which I did not realize until after sending out the program. I think I will just extricate the pnt code and compile that separately and that should be fine. Thanks very much again, to you and everybody else who replied. Best wishes, Ranjan On Wed, 21 Mar 2007

Re: [R] RMySQL *was* working...

2007-03-21 Thread Pete Cap
List, Please disregard. I merely forgot the quotes: > dbDriver("MySQL") Thanks, Pete Pete Cap <[EMAIL PROTECTED]> wrote: List, Last week with the help of Uwe and some other folks I was able to get RMySQL 0.5-7 compiled against R 2.4.1 and MySQL 5.0.27. It was working fine--I was able to se

Re: [R] any way to append a table in SQL server

2007-03-21 Thread Jim Porzak
see... > library(RODBC) > ?sqlSave and rest of RODBC docs. You could also pass a SQL INSERT statement through sqlQuery() On 3/20/07, Wensui Liu <[EMAIL PROTECTED]> wrote: > Dear Lister, > Is there an interface in R with SQL server that allows me to append > records to table in the DB? Might I do

[R] RMySQL *was* working...

2007-03-21 Thread Pete Cap
List, Last week with the help of Uwe and some other folks I was able to get RMySQL 0.5-7 compiled against R 2.4.1 and MySQL 5.0.27. It was working fine--I was able to send select queries to the db, put the results in a data frame, and so forth. Today, dbDriver() threw an error: > dbDriver(My

Re: [R] how to get "lsmeans"?

2007-03-21 Thread Frank E Harrell Jr
Chuck Cleland wrote: > Liaw, Andy wrote: >> I verified the result from the following with output from JMP 6 on the >> same data (don't have SAS: don't need it): >> >> set.seed(631) >> n <- 100 >> dat <- data.frame(y=rnorm(n), A=factor(sample(1:2, n, replace=TRUE)), >> B=factor(sam

Re: [R] how to get "lsmeans"?

2007-03-21 Thread Prof Brian Ripley
On Wed, 21 Mar 2007, Chuck Cleland wrote: > Liaw, Andy wrote: >> I verified the result from the following with output from JMP 6 on the >> same data (don't have SAS: don't need it): >> >> set.seed(631) >> n <- 100 >> dat <- data.frame(y=rnorm(n), A=factor(sample(1:2, n, replace=TRUE)), >>

Re: [R] for loop in assigning column names

2007-03-21 Thread Dimitris Rizopoulos
try df.new <- cbind(df, ObJeCt[1:10]) names(df.new) <- c(names(df), paste("St", 1:10, sep = "")) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel:

Re: [R] question on suppressing error messages with Rmath library

2007-03-21 Thread Ranjan Maitra
Thanks, Jim and Ronggui! This would work, of course! But I apologize for not being clear in the first post. I am calling Rmath from a C program. I could not figure out how to set this without going back and recompiling: the message in in pnt.c and pulls it from the header file nmath.h. Thanks,

Re: [R] Gaussian Adaptive Quadrature

2007-03-21 Thread Ravi Varadhan
Hi, The function "integrate" is based on QUADPACK Fortran package by Piessens. It does indeed perform adaptive Gauss-Kronrod quadrature, where Kronrod's method allows the re-use of abscissa from the previous iteration, thus enabling the estimation of the quadrature error and its control. In contr

Re: [R] for loop in assigning column names

2007-03-21 Thread ONKELINX, Thierry
Have you tried something like df <- cbind(df, ObJeCT[, 1:10]) Cheers, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature and Forest Cel biometrie, methodologie en kwalit

[R] for loop in assigning column names

2007-03-21 Thread Luis Ridao Cruz
R-help, I have a data frame (df) and I want to add some columns whose names should correspond to the "i" index in the loop below. for(i in 1:10) { df$eval(paste("St", as.character(i), sep = "" )) <- ObJeCt[i] } An error message comes out : "Error: attempt to apply non-function" How can I

Re: [R] question on suppressing error messages with Rmath library

2007-03-21 Thread Luke Tierney
You might get less noise in the replies if you were explicit about using Rmath stand-alone and asked on r-devel. As far as I can see you would need to compile a version of the stand-alone library that defines the macros for handling of warning messages differently -- the current one just calls pri

Re: [R] package:AlgDesign and .Random.seed

2007-03-21 Thread Uwe Ligges
So let me read the help page of ?RNG for all of us: "Initially, there is no seed; a new one is created from the current time when one is required. Hence, different sessions will give different simulation results, by default." Hence this is definitely a bug in AlgDesign. Further on also remarka

Re: [R] Select the last two rows by id group

2007-03-21 Thread Muenchen, Robert A (Bob)
Marc, thanks for so many great variations! I especially like: tail(sort(table(DF$County))) I often have frequency tables that are of interest only towards the end. Cheers, Bob = Bob Muenchen (pronounced Min'-chen), Manager Statistica

[R] ar predict without noise

2007-03-21 Thread Michael John Shaw
Hi. I am just getting started with R. I would like to start by using one of the autoregressive models (beginning at the very basic, perhaps with just ar) to forecast a multivariate time series one step ahead without noise and, while digging around for how to do so, maybe someone can just tell me

[R] He wilma

2007-03-21 Thread Sonia KMohammed
Check this company out C E O AMERICA I N C SYm-CE OA Currently : 6 Cents, CHEAP!!! Add this to your radar AN ALL AMERICAN COMPANY Get IN Before the rush TOMORROW led the Pirates to a 68-55 record, a trip to the round of 16 in the 2000 games and eight of their next nine. ''This does a lot for

Re: [R] package:AlgDesign and .Random.seed [Broadcast]

2007-03-21 Thread Liaw, Andy
From: Michael Kubovy > > On Mar 21, 2007, at 4:16 AM, Uwe Ligges wrote: > > > Michael Kubovy wrote: > >> Dear r-helpers, > >> Could you please help me solve the following problem: When I run > >> require(AlgDesign) > >> trt <- LETTERS[1:5] > >> blk <- 10 > >> trtblk <- 3 > >> BIB <- optBlock(~.,

Re: [R] question on suppressing error messages with Rmath library

2007-03-21 Thread James W. MacDonald
Hi Ranjan, Ranjan Maitra wrote: > Dear list, > > I have been using the Rmath library for quite a while: in the current > instance, I am calling dnt (non-central t density function) > repeatedly for several million. When the argument is small, I get the > warning message: > > full precision was no

Re: [R] newbie upgrade from 2.4.0 to 2.4.1 question

2007-03-21 Thread Prof Brian Ripley
On Wed, 21 Mar 2007, Mark Fisher wrote: > I am new to R. I installed version 2.4.0 some time ago and I find that > some packages I want to use require 2.4.1. I am using Windows XP. Do I > need to uninstall R first before running the setup file for 2.4.1 or > does the setup file "do the right thing

Re: [R] question on suppressing error messages with Rmath library

2007-03-21 Thread ronggui
op <- options(warn=-1) [main codes here] options(op) On 3/21/07, Ranjan Maitra <[EMAIL PROTECTED]> wrote: > Dear list, > > I have been using the Rmath library for quite a while: in the current > instance, I am calling dnt (non-central t density function) repeatedly for > several million. When

[R] question on suppressing error messages with Rmath library

2007-03-21 Thread Ranjan Maitra
Dear list, I have been using the Rmath library for quite a while: in the current instance, I am calling dnt (non-central t density function) repeatedly for several million. When the argument is small, I get the warning message: full precision was not achieved in 'pnt' which is nothing unexpect

Re: [R] newbie upgrade from 2.4.0 to 2.4.1 question

2007-03-21 Thread Alberto Monteiro
Mark Fisher wrote: > > I am new to R. I installed version 2.4.0 some time ago and I find > that some packages I want to use require 2.4.1. I am using Windows > XP. Do I need to uninstall R first before running the setup file for > 2.4.1 or does the setup file "do the right thing"? > You don't

Re: [R] package:AlgDesign and .Random.seed

2007-03-21 Thread Michael Kubovy
On Mar 21, 2007, at 4:16 AM, Uwe Ligges wrote: > Michael Kubovy wrote: >> Dear r-helpers, >> Could you please help me solve the following problem: When I run >> require(AlgDesign) >> trt <- LETTERS[1:5] >> blk <- 10 >> trtblk <- 3 >> BIB <- optBlock(~., withinData = trt, blocksizes = rep(trtblk, b

Re: [R] how to get "lsmeans"?

2007-03-21 Thread Chuck Cleland
Liaw, Andy wrote: > I verified the result from the following with output from JMP 6 on the > same data (don't have SAS: don't need it): > > set.seed(631) > n <- 100 > dat <- data.frame(y=rnorm(n), A=factor(sample(1:2, n, replace=TRUE)), > B=factor(sample(1:2, n, replace=TRUE)), >

[R] substitution matrix for amino acid

2007-03-21 Thread Nora Muda
Hi, I want to compute the pairwise distances of amino acid sequences and then construct the phylogenetic tree. I noticed that ape package support only for DNA sequences and compute the pairwise distances by using dist.dna function. Can anyone suggest what package that support for protein sequenc

Re: [R] Gaussian Adaptive Quadrature

2007-03-21 Thread Douglas Bates
On 3/21/07, Doran, Harold <[EMAIL PROTECTED]> wrote: > The function integrate() uses AGQ. There are other functions for > gaussian quadrature in the statmod() package that I really like. I think that integrate does adaptive quadrature but not adaptive Gaussian quadrature (which probably should hav

[R] newbie upgrade from 2.4.0 to 2.4.1 question

2007-03-21 Thread Mark Fisher
Hi, I am new to R. I installed version 2.4.0 some time ago and I find that some packages I want to use require 2.4.1. I am using Windows XP. Do I need to uninstall R first before running the setup file for 2.4.1 or does the setup file "do the right thing"? --Mark.

Re: [R] how to get "lsmeans"?

2007-03-21 Thread Liaw, Andy
I verified the result from the following with output from JMP 6 on the same data (don't have SAS: don't need it): set.seed(631) n <- 100 dat <- data.frame(y=rnorm(n), A=factor(sample(1:2, n, replace=TRUE)), B=factor(sample(1:2, n, replace=TRUE)), C=factor(sample

Re: [R] Over-writing functions from other packages? What is a good strategy??

2007-03-21 Thread Gabor Grothendieck
Sometimes but its also easy to forget about simple solutions. On 3/21/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote: > Gabor Grothendieck wrote: > > > > Could you call yours Edges? > > > Is it a good idea to have two different functions, whose name > differs by case sensitivity? I believe this is

Re: [R] Detailed legend in mathplot ...

2007-03-21 Thread Petr Pikal
Hi you can put quite sophisticated legend using legend command e.g. legend(locator(1), legend=c("first", "second"), pch=c(21,22), lty = c(2,1)) Regards Petr On 21 Mar 2007 at 8:38, Petar Milin wrote: From: Petar Milin <[EMAIL PROTECTED]> To: r-help Organ

Re: [R] Detailed legend in mathplot ...

2007-03-21 Thread Uwe Ligges
Petar Milin wrote: > You are absolutely right, and I apologize. Here is the data.frame: > no day concentrated distributed > 1 4 7 6.8 > 2 5 3 6 > 3 6 2.3 5.65 > 4 7 1.6 5.30 > 9 8 0.9 4.95 > > And then, if DAT = data.frame, mathplot() should do the job. So what is the problem with creating the

Re: [R] Gaussian Adaptive Quadrature

2007-03-21 Thread Doran, Harold
The function integrate() uses AGQ. There are other functions for gaussian quadrature in the statmod() package that I really like. > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Caio > Lucidius Naberezny Azevedo > Sent: Wednesday, March 21, 2007 5

[R] Integer Optimization

2007-03-21 Thread Mahieux Dimitri
Hello everybody, I am looking for a function (or package) which can solve a integer optimization problem. I have found the glpk package but I don't known very well how to use it. Someone has another solution ? thx __ R-help@stat.math.ethz.ch mailing

Re: [R] Detailed legend in mathplot ...

2007-03-21 Thread Petar Milin
You are absolutely right, and I apologize. Here is the data.frame: no day concentrated distributed 1 4 7 6.8 2 5 3 6 3 6 2.3 5.65 4 7 1.6 5.30 9 8 0.9 4.95 And then, if DAT = data.frame, mathplot() should do the job. Thank you in advance, Petar On Wed, 2007-03-21 at 09:37 +0100, Uwe Ligges wro

Re: [R] Over-writing functions from other packages? What is a good strategy??

2007-03-21 Thread Alberto Monteiro
Gabor Grothendieck wrote: > > Could you call yours Edges? > Is it a good idea to have two different functions, whose name differs by case sensitivity? I believe this is a shortcut to Chaos :-) IMHO, case sensitivity should _only_ be used for aesthetical purposes, else it may result in errors that

[R] Gaussian Adaptive Quadrature

2007-03-21 Thread Caio Lucidius Naberezny Azevedo
Hi all, Does anybody know any function that performs gaussian adapative quadrature integration of univariate functions? Thanks in advance, Regards, Caio __ [[alternative HTML version deleted]] ___

Re: [R] Ticks on barplots

2007-03-21 Thread Jim Lemon
Michael H. Prager wrote: > Dear Gurus, > > Using R 2.4.1 on Windows XP > > I am generating stacked barplots of age-composition of fish populations > (Y) over time (X). As there are many years, not every bars is labeled. > When looking at the plot, it becomes difficult to associate labels with

Re: [R] multcomp

2007-03-21 Thread Chuck Cleland
Nair, Murlidharan T wrote: > Thanks for the hint. > Here is what I had done earlier according to the old package > > mult.comp<-simint(a~b, data=z, conf.level=0.99, type=c("Tukey") ) > isoforms<-as.vector(rownames(mult.comp$estimate)) > estimate<-as.vector(mult.comp$estimate) > lower<-as.vec

Re: [R] Strange integer result on Debian/amd64

2007-03-21 Thread Dave Ewart
On Wednesday, 21.03.2007 at 09:29 +, Prof Brian Ripley wrote: > >OK, I understand the problem now. I was wondering if it something along > >those lines, but initially dismissed that as a possibility since I > >expected 'choose' to be an integer calculation. > > It is not, because it would be

Re: [R] Strange integer result on Debian/amd64

2007-03-21 Thread Prof Brian Ripley
On Wed, 21 Mar 2007, Dave Ewart wrote: > On Wednesday, 21.03.2007 at 14:07 +0800, Berwin A Turlach wrote: > >>> All help/suggestions/appreciated! >> >> The calculations are done in floating point arithmetic, not integer >> arithmetic. From the help page on `choose' one might already guess >> so mu

Re: [R] Strange integer result on Debian/amd64

2007-03-21 Thread Dave Ewart
On Wednesday, 21.03.2007 at 09:10 +, Dave Ewart wrote: > However, the underlying problem that gave rise to the difficulty was as > follows. A colleague wishes to create a matrix, where one of the > dimensions of the matrix is the result of the 'choose' function, i.e. > > mycols<-choose(11,

Re: [R] Strange integer result on Debian/amd64

2007-03-21 Thread Dave Ewart
On Wednesday, 21.03.2007 at 14:07 +0800, Berwin A Turlach wrote: > > All help/suggestions/appreciated! > > The calculations are done in floating point arithmetic, not integer > arithmetic. From the help page on `choose' one might already guess > so much, but reading "R-2.4.1/src/nmath/choose.c" d

Re: [R] Detailed legend in mathplot ...

2007-03-21 Thread Uwe Ligges
Petar Milin wrote: > Hello, > Recently, I have asked for a help with building graphs, and I got few > great advices. Now, my appetite is growing :) and I wander how to add > legend for two (or more) lines in following example: > > matplot(DAT[, c(3,4)], type="b", ylim=c(0,8), xaxt="n", yaxt="n",

Re: [R] package:AlgDesign and .Random.seed

2007-03-21 Thread Uwe Ligges
Michael Kubovy wrote: > Dear r-helpers, > > Could you please help me solve the following problem: When I run > > require(AlgDesign) > trt <- LETTERS[1:5] > blk <- 10 > trtblk <- 3 > BIB <- optBlock(~., withinData = trt, blocksizes = rep(trtblk, blk)) > > In response to the last command, R comp

[R] how to get "lsmeans"?

2007-03-21 Thread Xingwang Ye
Dear all, I search the mail list about this topic and learn that no simple way is available to get "lsmeans" in R as in SAS. Dr.John Fox and Dr.Frank E Harrell have given very useful information about "lsmeans" topic. Dr. Frank E Harrell suggests not to think about lsmeans

Re: [R] error installing packages

2007-03-21 Thread Prof Brian Ripley
On Tue, 20 Mar 2007, Bernardo Rangel Tura wrote: > On Wed, 2007-03-07 at 09:33 -0800, Bricklemyer, Ross S wrote: >> I was finally able to get R to 'configure', 'make', and 'install' on >> Mandriva 2007. Itried to install gnomeGUI and I received an error. See >> below. At what step do I make R