[R] an array of matrices

2007-05-21 Thread Scott Hyde
I'd like to have a three dimensional array of matrices. I thought I could construct a five dimensional array to have the three dimensional array of matrices. However, not all of the matrices in the array have the same dimensions, which seems to mean I can't use a five dimensional array.

[R] [OT/R+MikTeX] WinAnsi.enc fonts not found by ghostscript by default [solved]

2007-05-21 Thread Peter Ruckdeschel
Concerns: OS Win32 / R 2.5.0 Recently (IIRC not with R 2.3.0), I had some problems including postscript graphics produced in R into LaTeX with MikTeX [version 2.5.2574]: yap, respectively GhostScript/GSView (versions 8.54/4.7 resp.) threw errors. In fact, I realized this was just the same

Re: [R] an array of matrices

2007-05-21 Thread jim holtman
Create a matrix of list: A = matrix(1:4,2,2) B = matrix(1:25,5,5) C = matrix(1,3,3) D = matrix(1,4,4) x - list(A,B,C,D) dim(x) - c(2,2) x[[1,1]] [,1] [,2] [1,]13 [2,]24 x[[2,2]] [,1] [,2] [,3] [,4] [1,]1111 [2,]1111 [3,]11

Re: [R] an array of matrices

2007-05-21 Thread Petr Klasterecky
Scott Hyde napsal(a): I'd like to have a three dimensional array of matrices. I thought I could construct a five dimensional array to have the three dimensional array of matrices. However, not all of the matrices in the array have the same dimensions, which seems to mean I can't use a

Re: [R] Running an R script without running R

2007-05-21 Thread Uwe Ligges
Alberto Vieira Ferreira Monteiro wrote: Is there any way to run an R script without running R? As an example, suppose I have a tcl/tk interface that asks for a number (in a GUI) and displays its factorial. Is there a way to invoke this script without invoking R? I'm using R 2.4.1 in

Re: [R] Execute expression ( R -e) without close

2007-05-21 Thread Marcel.
Thank you for the fast answer ! -- View this message in context: http://www.nabble.com/Execute-expression-%28-R--e%29-without-close-tf3777652.html#a10715676 Sent from the R help mailing list archive at Nabble.com. __ R-help@stat.math.ethz.ch

[R] question about forward library

2007-05-21 Thread Manuele Pesenti
Dear R User, I'm a new member of this list, I'm using R language from a short time and I would like to ask you something I cannot solve... I'm using the fwdlm function from the forward library, I'm specially interested by the results of Cook's distances: I'd want to know which inputs are the

[R] Making links in per-session dir

2007-05-21 Thread Ronaldo Reis Junior
Hi, When I use help.start() R spends some time to make links in per-session dir. I have a hole R system installed. It is possible to make these links permanent to avoid spend this time in help.start()? Thanks Ronaldo -- There can be no twisted thought without a twisted molecule.

Re: [R] Running an R script without running R

2007-05-21 Thread Thomas Adams
Alberto, Below I have two scripts: (1) cpc2fgroup (2) R.cpc.6_10day.outlook.batch Bash shell script (1) calls the R script (2). (1) also sets some environment variables for (2). I make no claims that this is good coding and I might write things a little differently now (rather than a few

[R] more simplified output from glht object

2007-05-21 Thread Ronaldo Reis Junior
Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo -- Human kind cannot bear very much reality. -- T. S. Eliot, Four Quartets: Burnt Norton -- Prof.

Re: [R] using lm() with variable formula

2007-05-21 Thread Vladimir Eremeev
I was solving similar problem some time ago. Here is my script. I had a data frame, containing a response and several other variables, which were assumed predictors. I was trying to choose the best linear approximation. This approach now seems to me useless, please, don't blame me for that.

Re: [R] Sample correlation coefficient question NOT R question

2007-05-21 Thread Bruce Willy
Hmm I'm not too sure but the correlation coefficient must not be taken for its empirical counterpart and the law of large numbers tells you roughly you can approximate a mean by its empirical counterpart when the variables are identically distributed and independant. If they are not

Re: [R] Running an R script without running R

2007-05-21 Thread Alberto Monteiro
Thomas Adams wrote: Below I have two scripts: (1) cpc2fgroup (2) R.cpc.6_10day.outlook.batch Bash shell script (1) calls the R script (2). It works with R CMD BATH r-script I know how to do this (I do it inside make...) I guess I did not make myself clear. I want to write some GUI

[R] Chosing a subset of a non-sorted vector

2007-05-21 Thread Christoph Scherber
Dear all, I have a tricky problem here: I have a dataframe with biodiversity data in which suplots are a repeated sequence from 1 to 4 (1234,1234,...) Now, I want to randomly pick two subplots each from each diversity level (DL). The problem is that it works up to that point - but if I try

[R] ... and 5 cents more.

2007-05-21 Thread Vladimir Eremeev
Vladimir Eremeev wrote: I was solving similar problem some time ago. Here is my script. I had a data frame, containing a response and several other variables, which were assumed predictors. I was trying to choose the best linear approximation. This approach now seems to me useless,

[R] accessing unexported functions

2007-05-21 Thread Erich Neuwirth
I need to access some internal function of a package temporarily. These functions normally that are not visible outside of the package. I know the MyPackage:::myfunction mechanism, but that does not solve my problem. Instead, I need to add the internal namespace of the package to the search path,

Re: [R] Making links in per-session dir

2007-05-21 Thread Prof Brian Ripley
The short answer is no. When you say you have the 'hole R system installed', what did you mean? This is not slow if you have all of R installed, but it might be if you have all of CRAN installed. But how do you know that is where the time is being spent? You get a message that it is making

[R] Questions about bwplot

2007-05-21 Thread Klaus Nordhausen
Dear R-experts, I have some questions about boxplots with lattice. My data is similar as in the example below, I have two factors (Goodness of Fit and Algorithms) and data values but in each panels the scales are quite different, therefore the normal boxplots produced by set.seed(1) GOF -

[R] title line

2007-05-21 Thread elyakhlifi mustapha
hello, I know how to do to write down the columns in a data frame but I also wanna write down the lines . Do you know how could I do this? thanks _ [[alternative HTML version deleted]]

Re: [R] FW: R

2007-05-21 Thread Ricky Rankin
AJ, Can you test and let us know if this is OK Ricky __ Principal Analyst Information Services Queen's University Belfast   tel: 02890 974824 fax: 02890 976586 email: [EMAIL PROTECTED] -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: 21

Re: [R] FW: R

2007-05-21 Thread Peter Dalgaard
Ricky Rankin wrote: AJ, Can you test and let us know if this is OK Erm, I think you might want to keep me and a few thousand readers of r-help out of the loop at this stage... -pd Ricky __ Principal Analyst Information Services Queen's University Belfast

Re: [R] can I get same results using lme and gls?

2007-05-21 Thread Douglas Bates
On 5/20/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: I was wondering how to get the same results with gls and lme. In my lme, the design matrix for the random effects is (should be) a identity matrix and therefore G should add up with R to produce the R matrix that gls would report

Re: [R] title line

2007-05-21 Thread Michael Kubovy
On May 21, 2007, at 8:21 AM, elyakhlifi mustapha wrote: I know how to do to write down the columns in a data frame but I also wanna write down the lines . Do you know how could I do this? I'm not sure what you mean by write down the lines. If you mean that you want to know what the row

[R] absent

2007-05-21 Thread bcutayar
Je serai absent(e) à partir du 21/05/2007 de retour le 28/05/2007. En mon absence, vous pouvez contacter M.Bidart (5503). Si vous n'etes pas destinataires de ce message, merci d'aver...{{dropped}} __ R-help@stat.math.ethz.ch mailing list

Re: [R] Questions about bwplot

2007-05-21 Thread Richard M. Heiberger
## I simplified this panel function by taking the calculations outside ## the panel function to allow use of the original variable names. my.boxpanel3 - function(x, y, subscripts, ..., medians, best) { cols - grey(1-best[,panel.number()]/3) panel.bwplot(x, y, fill=cols, ...)

Re: [R] Questions about bwplot

2007-05-21 Thread Charilaos Skiadas
On May 21, 2007, at 8:39 AM, Klaus Nordhausen wrote: Dear R-experts, I have some questions about boxplots with lattice. My data is similar as in the example below, I have two factors (Goodness of Fit and Algorithms) and data values but in each panels the scales are quite different,

Re: [R] more simplified output from glht object

2007-05-21 Thread Ronaldo Reis Junior
Em Segunda 21 Maio 2007 08:14, Ronaldo Reis Junior escreveu: Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo I try to use multcompLetters, but dont work Look

Re: [R] Making links in per-session dir

2007-05-21 Thread Ronaldo Reis Junior
Em Segunda 21 Maio 2007 08:31, Prof Brian Ripley escreveu: The short answer is no. When you say you have the 'hole R system installed', what did you mean? This is not slow if you have all of R installed, but it might be if you have all of CRAN installed. But how do you know that is where

[R] help with this indexing

2007-05-21 Thread javier garcia-pintado
Hi all, Let's say I have a long data frame and a short one, both with three colums: $east, $north, $value And I need to fill in the short$value, extracting the corresponding value from long$value, for coinciding $east and $north in both tables. I know the possibility: for (i in

[R] Boostrap p-value in regression [indirectly related to R]

2007-05-21 Thread Viechtbauer Wolfgang \(STAT\)
Hello All, Despite my preference for reporting confidence intervals, I need to obtain a p-value for a hypothesis test in the context of regression using bootstrapping. I have read John Fox's chapter on bootstrapping regression models and have consulted Efron Tibshirani's An Introduction to the

[R] Equivalent of 'Retain' statement in R

2007-05-21 Thread Shubha Vishwanath Karanth
Hi R, What is the equivalent keyword for 'retain' statement of SAS in R? Thank you, Shubha [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] How to create an R RPM from source

2007-05-21 Thread Waichler, Scott R
I use some computers that run older versions of Redhat Linux such as EL3. RPMs for the current version of R are no longer provided via CRAN for these older operating system versions. How can I compile my own RPM for the current version of R, which I could then use to install R quickly and easily

Re: [R] How to compare linear models with intercept and those withoutintercept using minimizing adjs R^2 strategy

2007-05-21 Thread Lucke, Joseph F
One issue is whether you want your estimators to be based on central moments (covariances) or on non-central moments. Removing the intercept changes the statistics from central to non-central moments. The adjusted R2, by which I think you mean Fisher's adjusted R2, is based on central moments

Re: [R] accessing unexported functions

2007-05-21 Thread Seth Falcon
Erich Neuwirth [EMAIL PROTECTED] writes: I need to access some internal function of a package temporarily. These functions normally that are not visible outside of the package. I know the MyPackage:::myfunction mechanism, but that does not solve my problem. Instead, I need to add the

[R] How to create an R RPM from source

2007-05-21 Thread R P Herrold
On Mon, 21 May 2007, Waichler, Scott R wrote: I use some computers that run older versions of Redhat Linux such as EL3. RPMs for the current version of R are no longer provided via CRAN for these older operating system versions. How can I compile my own RPM for the current version of R,

Re: [R] help with this indexing

2007-05-21 Thread Adaikalavan Ramasamy
merge() javier garcia-pintado wrote: Hi all, Let's say I have a long data frame and a short one, both with three colums: $east, $north, $value And I need to fill in the short$value, extracting the corresponding value from long$value, for coinciding $east and $north in both tables. I know

Re: [R] How to compare linear models with intercept and those withoutintercept using minimizing adjs R^2 strategy

2007-05-21 Thread 李俊杰
2007/5/21, Lucke, Joseph F [EMAIL PROTECTED]: One issue is whether you want your estimators to be based on central moments (covariances) or on non-central moments. Removing the intercept changes the statistics from central to non-central moments. The adjusted R2, by which I think you mean

[R] coefficients regression

2007-05-21 Thread romi
Hi, I would like to calculate a polynomial regression with R, but I don't get the same coefficients as when using SPSS. Is there a way to transform the coefficients? I use: regression - lm (biomass ~ poly (temperature, 2)) Thank you, Romana Limberger

[R] Source code of add1

2007-05-21 Thread Antonio José Sáez Castillo
We are developing a new regression model. We now need to develop an analogous to add1 in stat package. My problem is about update.formula, but I think I could solve this problem if I can get source code of add1. So, can I get the source code of add1? Prof.

[R] sequentially process a list

2007-05-21 Thread Robert
Hi dear R users, I'm a R beginner and I have a basic question about sequential treatments of lists. I have a time based (i.e. events are consecutive) list of values of a biological property. Like : time value 15 210 3 7 4 10 5 19 6 21 7 20 8

[R] comparing fit of cubic spline

2007-05-21 Thread Thomas Reed
I want to compare the fit of a quadratic model to continuous data, with that of a cubic spline fit. Is there a way of computing AIC from for e.g. a GAM with a smoothing spine, and comparing this to AIC from a quadratic model? Cheers ** Tom Reed PhD

[R] generate multivariate non-normal distribution

2007-05-21 Thread Jiao Yang
I was wondering whether there's a code in R to generate a multivariate non-normal distribution for a given covariance matrix other than rmvt. Thanks a lot! Regards, Jiao __ R-help@stat.math.ethz.ch mailing list

Re: [R] How to create an R RPM from source

2007-05-21 Thread Waichler, Scott R
I use some computers that run older versions of Redhat Linux such as EL3. RPMs for the current version of R are no longer provided via CRAN for these older operating system versions. How can I compile my own RPM for the current version of R, which I could then use to install R

Re: [R] coefficients regression

2007-05-21 Thread Viechtbauer Wolfgang \(STAT\)
Try: regression - lm (biomass ~ poly (temperature, degree=2, raw=TRUE)) See the help page for poly what raw=TRUE does. Best, -- Wolfgang Viechtbauer  Department of Methodology and Statistics  University of Maastricht, The Netherlands  http://www.wvbauer.com/ Original Message From:

Re: [R] How to compare linear models with intercept and those withoutintercept using minimizing adjs R^2 strategy

2007-05-21 Thread 李俊杰
I have a question about what you've wrote in your pdf file. Why must we view my problem in the viewpoint of hypothesis testing? Is testing the original philosophy of maximizing Fisher's A-statistic to choose a optimum model? Thanks. 2007/5/21, Lucke, Joseph F [EMAIL PROTECTED]: I taken the

Re: [R] Equivalent of 'Retain' statement in R

2007-05-21 Thread Peter Dalgaard
Shubha Vishwanath Karanth wrote: What is the equivalent keyword for 'retain' statement of SAS in R? There is none. This is strongly linked to the sequential operation of the DATA step i SAS and R just does things differently. The lag() function sometimes helps, but you may need an explicit

Re: [R] comparing fit of cubic spline

2007-05-21 Thread Ben Bolker
Thomas Reed tom.reed at ed.ac.uk writes: I want to compare the fit of a quadratic model to continuous data, with that of a cubic spline fit. Is there a way of computing AIC from for e.g. a GAM with a smoothing spine, and comparing this to AIC from a quadratic model? Why don't you try

[R] Stratified Cox proportional Hazard Model

2007-05-21 Thread Dora Ilieva
Hello everyone, I am a new user of R. Does anybody know how hazard ratios are extracted for each factor level in a stratified Cox proportional hazard regression model? I have a cancer data set where the variable ?differentiation? is a factor with three levels: poor, intermediate and good. I would

Re: [R] density

2007-05-21 Thread Emilio Tortosa Ausina
There's a nice package ('ks') which even allows you to specify a matrix of bandwidths (not only one bandwidth for each coordinate direction). Hope this helps, Emili Missatge citat per Bruce Willy [EMAIL PROTECTED]: Hello, I have a n*2 matrix, called plan, which contains n observations

Re: [R] Source code of add1

2007-05-21 Thread Duncan Murdoch
On 21/05/2007 4:57 AM, Antonio José Sáez Castillo wrote: We are developing a new regression model. We now need to develop an analogous to add1 in stat package. My problem is about update.formula, but I think I could solve this problem if I can get source code of add1. So, can I get the

Re: [R] Source code of add1

2007-05-21 Thread Uwe Ligges
Antonio José Sáez Castillo wrote: We are developing a new regression model. We now need to develop an analogous to add1 in stat package. My problem is about update.formula, but I think I could solve this problem if I can get source code of add1. So, can I get the source code of add1?

[R] quartz() on MAC OSX

2007-05-21 Thread Rolf Turner
I am (desperately) trying to get used to using a Mac here at my new location. (Why *anyone* would ever use anything other than Linux, except under duress as I am, totally escapes me, but that's another story.) Fortunately much of the Mac OSX is actually Unix, so a civilized person can manage to

Re: [R] sequentially process a list

2007-05-21 Thread jim holtman
This should do it for you. It makes use of 'rle' to find the 'runs' of thresholds and then changes the values to the sequence numbers x - time value + 15 + 210 + 3 7 + 4 10 + 5 19 + 6 21 + 7 20 + 8 18 + 9 10 + 10 7 + 11 8 + 12

Re: [R] sequentially process a list

2007-05-21 Thread Mike Lawrence
Try nested ifelse() statements to label the group ex. time$group=ifelse(time$time5,NA, ifelse(time$time9,1, ifelse(time$time13,NA, ifelse(time$time17,2,NA) ) ) ) Then use aggregate to find the max value. ex.

[R] How to create an R RPM from source

2007-05-21 Thread R P Herrold
On Mon, 21 May 2007, Waichler, Scott R wrote: I should have said that I do have root privileges. root privleges or not, 'better' practice mandates against building as root to avoid inadvertent pollution of one's permanent environment; non-root building is the only sensible way to go, short

Re: [R] Equivalent of 'Retain' statement in R

2007-05-21 Thread Daniel Nordlund
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Shubha Vishwanath Karanth Sent: Monday, May 21, 2007 7:40 AM To: r-help Subject: [R] Equivalent of 'Retain' statement in R Hi R, What is the equivalent keyword for 'retain' statement of SAS in

Re: [R] R2 always increases as variables are added?

2007-05-21 Thread Paul Lynch
Junjie, First, a disclaimer: I am not a statistician, and have only taken one statistics class, but I just took it this Spring, so the concepts of linear regression are relatively fresh in my head and hopefully I will not be too inaccurate. According to my statistics textbook, when

Re: [R] question about forward library

2007-05-21 Thread Manuele Pesenti
I answer about my own question... On Monday 21 May 2007 12:40, Manuele Pesenti wrote: Dear R User, I'm a new member of this list, I'm using R language from a short time and I would like to ask you something I cannot solve... I'm using the fwdlm function from the forward library, I'm

Re: [R] How to create an R RPM from source

2007-05-21 Thread Martyn Plummer
On Mon, 2007-05-21 at 07:46 -0700, Waichler, Scott R wrote: I use some computers that run older versions of Redhat Linux such as EL3. RPMs for the current version of R are no longer provided via CRAN for these older operating system versions. How can I compile my own RPM for the current

Re: [R] How to compare linear models with intercept and those withoutintercept using minimizing adjs R^2 strategy

2007-05-21 Thread 李俊杰
So when I am using the adjusted R2 and as a penalized optimality criterion, and I have to compare models with intercept and those without intercept to decide the final model selected, does my crierion in my first email make sense? Because we know that in leaps(leaps), if we want to select a model

Re: [R] {10,20,30}={25,30,15}

2007-05-21 Thread Greg Snow
To decide on which way to do a one tailed test, you should really ask a few questions (before looking at the data) either of yourself, your client, or other expert in the field. I would start by asking the 3 questions: 1. What will I/you do if A is less than B? 2. What will I/you do if A is

[R] size limit in R?

2007-05-21 Thread j-simpson
Hi, Please see the email exchanges below. I am having trouble generating output that is large enough for our needs, specifically when using the GaussRF function. However, when I wrote Dr. Schlather (the author of the GaussRF function), he indicated that there is also a limit imposed by R

Re: [R] sequentially process a list

2007-05-21 Thread Mike Lawrence
Sorry, obviously misread your help request. Jim's solution is what you're looking for. On 21-May-07, at 12:56 PM, Mike Lawrence wrote: Try nested ifelse() statements to label the group ex. time$group=ifelse(time$time5,NA, ifelse(time$time9,1, ifelse(time$time13,NA,

Re: [R] R2 always increases as variables are added?

2007-05-21 Thread Alberto Monteiro
Paul Lynch wrote: I don't think it makes sense to compare models with and without an intercept term. (Also, I don't know what the point of using a model without an intercept term would be, but that is probably just my ignorance.) Suppose that you are 100% sure that the intercept term is

[R] plot problems

2007-05-21 Thread Markus voigt
Hello I'm a newbie to the R program and I have problems to plot two functions in one plot window. I know that I can use the plot command for the first function and the lines command for the second function but the problem is that the first function has only values beginnig from 1 to 7 on the

[R] Comparing multiple distributions

2007-05-21 Thread jiho
Hello eveybody, I am studying the vertical distribution of plankton and want to study its variations relatively to several factors (time of day, species, water column structure etc.). So my data is special in that, at each sampling site (each observation), I don't have *one* number, I have

[R] 'varFunc' class with additive variance? (was: can I get same results using lme and gls?)

2007-05-21 Thread Spencer Graves
Hi, Doug and others: What might be the best tools for modeling an additive variance structure for residuals, something like the following: var(resid) = s0.2*(1+var.pred) + daysFromTraining*var(process migration per day), where var.pred = relative variance of

[R] Installing packages from command line on Linux RHEL4

2007-05-21 Thread Kermit Short
Greetings. I am a System Administrator, and thus have very little knowledge of R itself. I have been asked to install a list of some 200 packages (from CRAM) to R. Rather than installing each package manually, I was hoping I could script this. I've written a BASH script that hopefully will

Re: [R] plot problems

2007-05-21 Thread Roland Rau
Markus voigt wrote: because it should start at 2. Is there a parameter where you can define the start on the x-axis? Maybe this will help you? datax1 - 1:8 datax2 - 2:9 datay1 - runif(length(datax1)) datay2 - runif(length(datax2)) plot(x=datax1, y=datay1, type=l, col=blue, xlab=X,

Re: [R] size limit in R?

2007-05-21 Thread Roland Rau
[EMAIL PROTECTED] wrote: Hi, Please see the email exchanges below. I am having trouble generating output that is large enough for our needs, specifically when using the GaussRF function. However, when I wrote Dr. Schlather (the author of the GaussRF function), he indicated that there

Re: [R] Installing packages from command line on Linux RHEL4

2007-05-21 Thread Kermit Short
Dirk- Many thanks for your reply. As I mentioned, I know very little about programming in 'R' and what I've got is a BASH script. If needs be, I'll look up how to read in a text file through R and add that into your script in lieu of the (argv) stuff, but you wouldn't happen to know how

Re: [R] Installing packages from command line on Linux RHEL4

2007-05-21 Thread Dirk Eddelbuettel
Hi Kernit, On 21 May 2007 at 11:37, Kermit Short wrote: | Greetings. | |I am a System Administrator, and thus have very little knowledge of R | itself. I have been asked to install a list of some 200 packages (from | CRAM) to R. Rather than installing each package manually, I was hoping I

Re: [R] more simplified output from glht object

2007-05-21 Thread Richard M. Heiberger
I don't yet understand what you are looking for. Can you post a reproducible example, preferably an example from ?glht, of what you are getting and of what you want. Rich __ R-help@stat.math.ethz.ch mailing list

Re: [R] quartz() on MAC OSX

2007-05-21 Thread Kyle.
I don't know if it's civilized enough for you, but running R from the x11 command line will allow you to open an x11() device without a problem. Kyle H. Ambert Graduate Student, Dept. Behavioral Neuroscience Oregon Health Science University [EMAIL PROTECTED] On May 20, 2007, at 11:14 PM,

Re: [R] plot cyrillic characters on MacOS X

2007-05-21 Thread Vlad Skvortsov
Prof Brian Ripley wrote: On Sat, 12 May 2007, Vlad Skvortsov wrote: Hi! I'm running R version 2.4.1 (2006-12-18) on MacOS X 10.4; my $LANG is ru_RU.KOI8-R. But we are told that all locales on MacOS X are actually UTF-8: it will not matter as R knows the charset it is using. I'm not sure

[R] quantreg and sparseM will not load

2007-05-21 Thread S . Hodgkinson
I have recently started using R and want to use the quantreg (and sparseM) packages. I downloaded the .tar files for each, and placed the subsequent folders into the library folder in the frameworks/R.framework/ resources/library folder with all the other packages. When I try to load either

Re: [R] can I get same results using lme and gls?

2007-05-21 Thread toby909
Douglas Bates wrote: On 5/20/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: I was wondering how to get the same results with gls and lme. In my lme, the design matrix for the random effects is (should be) a identity matrix and therefore G should add up with R to produce the R matrix that gls

Re: [R] quartz() on MAC OSX

2007-05-21 Thread hadley wickham
Hi Rolf, Two possible solutions: * DISPLAY=0.0.0.0:0 R - and then X11() should work without having to use xterm * install.packages(CarbonEL); library(CarbonEL); quartz() Hadley On 5/21/07, Rolf Turner [EMAIL PROTECTED] wrote: I am (desperately) trying to get used to using a Mac here at

[R] Graphically show population density for a specific geographic area using R

2007-05-21 Thread Liu, Ningwei
Hi, we all know that R is powerful in dealing with graphics. I am now trying to use R to show some geographic attributes for a specific geographic area on a graph. For example, if I would like to generate a graph showing the population densities (in terms of the intensity of colors) for

Re: [R] plot problems

2007-05-21 Thread Markus voigt
Hi Roland, thanks for your help. I can see that your exapmle do that what I want to have. Here my source code because your solution doesn't work for my data: currentarray etc. are vectors table(currentarray) currentarray 2 3 4 5 6 7 8 13 70 260 339 240 72 6

Re: [R] Graphically show population density for a specific geographicarea using R

2007-05-21 Thread Greg Snow
1. Yes you need a shape file with that information. 2. One place to look is: http://www.census.gov/geo/www/cob/bdy_files.html There are shapefiles (the .shp ones) for counties and zip codes (as well as many other things). These work pretty directly with the tools in the maptools and sp packages.

Re: [R] quartz() on MAC OSX

2007-05-21 Thread jiho
On 2007-May-21 , at 08:14 , Rolf Turner wrote: I am (desperately) trying to get used to using a Mac here at my new location. (Why *anyone* would ever use anything other than Linux, except under duress as I am, totally escapes me, but that's another story.) Oh that's harsh, Mac OS X is

Re: [R] lattice contourplot error

2007-05-21 Thread Deepayan Sarkar
On 5/20/07, Mihai Bisca [EMAIL PROTECTED] wrote: Hello, A weird thing is happening to me. I have a 128x128 elements matrix m with values between 0 and 1. If I do: R library(lattice) R contourplot (m, cuts=9) I get the following error message: Error in validDetails.text(x) :

[R] how to abort script with message

2007-05-21 Thread Blew, Ted
it is desired to abort an R script with a message, returning to the R prompt, pending 'if' results, as follows: first part of script . . if (condition) { action } else { 'error' abort } . . remainder of script --- note: 'stop' aborts the current

[R] PLS in R and SAS

2007-05-21 Thread Kyle Rogers
Dear all: I am comparing the PLS outputs of R and SAS for the following data set: Y x1 x2 x3 3 6 2 2 3 1 5 5 4 7 4 1 5 6 5 6 2 4 3 2 8 5 0 9 where Y is the dependent variable and x1, x2, x3 are the independent variables. I found several PLS algorithms in R

Re: [R] quartz() on MAC OSX

2007-05-21 Thread hadley wickham
On 5/21/07, jiho [EMAIL PROTECTED] wrote: On 2007-May-21 , at 08:14 , Rolf Turner wrote: I am (desperately) trying to get used to using a Mac here at my new location. (Why *anyone* would ever use anything other than Linux, except under duress as I am, totally escapes me, but that's

Re: [R] Graphically show population density for a specific geographicarea using R

2007-05-21 Thread Roger Bivand
On Mon, 21 May 2007, Greg Snow wrote: 1. Yes you need a shape file with that information. 2. One place to look is: http://www.census.gov/geo/www/cob/bdy_files.html There are shapefiles (the .shp ones) for counties and zip codes (as well as many other things). These work pretty directly

Re: [R] how to abort script with message

2007-05-21 Thread Henrik Bengtsson
Hi, I think the behavior that you outline is due to the fact that you cut'n'paste the script to the R prompt, is that correct? If so, use source() instead to run your script, then stop() will do what you want. /Henrik On 5/21/07, Blew, Ted [EMAIL PROTECTED] wrote: it is desired to abort an R

Re: [R] PLS in R and SAS

2007-05-21 Thread Lukasz Komsta
Dnia 2007-05-21 21:59, Użytkownik Kyle Rogers napisał: Any clarification on this matter will be greatly appreciated. The result obtained with R is calculated with centering the x1-x3 matrix around the column means. The SAS result is calculating without any preprocessing of data.

[R] A subscript out of bonds and write.table problem on manipulating a large size dataset

2007-05-21 Thread Yong Wang
Dear all: Described below is a large data set problem (data size 2G after unzipping, table delimited). I know R is not the appropriate tool for such task, anyway I did it on a server and get some straightforward problems. 1. The first is count.fields can count all the rows, however, when I

Re: [R] how to abort script with message

2007-05-21 Thread Blew, Ted
Henrik, yes, probably right, since the way i load the script is via the 'File/Open Script' command in the windows R interface. i tried 'File/Source R code...' and that works fine. thanks very much for your prompt response! ted. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

[R] Transfer entropy not implemented in R?

2007-05-21 Thread Atte Tenkanen
Hi, I couldn't find the 'transfer entropy' implemented in any R packages. I can't follow the math needed here (www.macalester.edu/~kaplan/knoxville/PRL00461.pdf). If it is not implemented, could somebody do it? Or write a garage level code for testing purposes? Atte Tenkanen

Re: [R] plot problems

2007-05-21 Thread Ben Bolker
Markus voigt markus-voigt at gmx.net writes: plot(table(currentarray)/1000,typ=b, col=c(red),ylim=c(0,0.8),xlim=c(1,8)) Not at all obvious but here's what you have to do to use the proper x values: t2 = table(currentarray2)/1000 x2 = as.numeric(names(t2)) t3 = table(probearray)/1000 x3 =

Re: [R] lattice contourplot error

2007-05-21 Thread Mihai Bisca
On 5/21/07, Deepayan Sarkar [EMAIL PROTECTED] wrote: I get the following error message: Error in validDetails.text(x) : Invalid 'rot' value This depends on the data (and I haven't been able to reproduce it), but I'm pretty sure the culprit is the calculation of the rotation for the contour

[R] Special characters on axis lables

2007-05-21 Thread Arne Brutschy
Hello, is there a way to specify special characters as a axis label? I need to have a tick marking the infinity in my plot... Is there a way? Or possibly a way to include latex-style formulas? Thanks in advance, Arne __ R-help@stat.math.ethz.ch

Re: [R] quartz() on MAC OSX

2007-05-21 Thread Rolf Turner
On 22/5/07 6:48 AM, hadley wickham [EMAIL PROTECTED] wrote: Hi Rolf, Two possible solutions: * DISPLAY=0.0.0.0:0 R - and then X11() should work without having to use xterm * install.packages(CarbonEL); library(CarbonEL); quartz() It is clear that life is determined to frustrate

Re: [R] lattice contourplot error

2007-05-21 Thread Deepayan Sarkar
On 5/21/07, Mihai Bisca [EMAIL PROTECTED] wrote: On 5/21/07, Deepayan Sarkar [EMAIL PROTECTED] wrote: I get the following error message: Error in validDetails.text(x) : Invalid 'rot' value This depends on the data (and I haven't been able to reproduce it), but I'm pretty sure the

[R] Selecting complementary colours

2007-05-21 Thread John Fox
Dear r-helpers, I wonder whether, given the #rrggbb representation of a colour, there is a simple way to select the complementary colour, also expressed as a #rrggbb string. Any suggestions would be appreciated. John John Fox, Professor Department of Sociology

[R] Need 64-bit integers on 32-bit platform

2007-05-21 Thread Earl F. Glynn
I'm using the latest R on Windows XP: R.version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 5.0 year 2007 month 04 day23 svn rev41293 language

Re: [R] Selecting complementary colours

2007-05-21 Thread Deepayan Sarkar
On 5/21/07, John Fox [EMAIL PROTECTED] wrote: Dear r-helpers, I wonder whether, given the #rrggbb representation of a colour, there is a simple way to select the complementary colour, also expressed as a #rrggbb string. Any suggestions would be appreciated. You want rgb2col. The following

  1   2   >