Re: [R] Using expression in Hmisc Key()

Michael Kubovy virginia.edu> writes: > > Dear r-helpers, > > How do I tell xYplot() and Key() that I want the labels in italic? > > Key(x = 0.667, y = 0.833, other = list(title = expression(italic(v)), > cex.title = 1, > labels = c(expression(italic(b)), expression(italic(c)), >

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

On Thu, 6 Dec 2007, Emmanuel Charpentier wrote: Prof Brian Ripley a écrit : Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit version of R somewhat moot. See chapter 8 of http://cran.r-project.org/doc/manuals/R-admin.html Thank you for this reminder|tip ! I didn't re-re

Re: [R] logistic regression using "glm",which "y" is set to be "1"

Dear all: By comparing glmresult$y and model.response(model.frame(glmresult)), I have found out which one is set to be "TRUE" and which "FALSE".But it seems that to fit a logistic regression , logit (or logistic) transformation has to be done before regression. Does anybody know how t

Re: [R] hclust in heatmap.2

Check the Rowv, Colv options to heatmap.2 data(mtcars) x <- as.matrix(mtcars) heatmap.2(x, Rowv=FALSE, dendrogram="column") -Ashoka Scientist - Pacific Northwest National Lab On Dec 5, 2007 4:20 PM, affy snp <[EMAIL PROTECTED]> wrote: > Dear list, > > I am using heatmap.2(x) to draw a hea

[R] Segmented regression

Hello all, I have 3 time series (tt) that I've fitted segmented regression models to, with 3 breakpoints that are common to all, using code below (requires segmented package). However I wish to specifiy a zero coefficient, a priori, for the last segment of the KW series (green) only. Is this possi

Re: [R] Working with "ts" objects

anscombe is built into R already so you don't need to read it in. An intercept is the default in lm so you don't have to specify it. opar <- par(mfrow = c(2,2)) plot(y1 ~ x1, anscombe) reg <- lm(y1 ~ x1, anscombe) reg abline(reg) ...etc... par(opar) Note that plot(anscombe[1:2]) and lm(anscombe[2

[R] snow package on multi core unix box

Is the rmpi package (or rpvm) needed to exploit multiple cores on a single unix box using the snow package. The documentation of the package does not provide info about setting up a single machine with multiple cores. Also, if how effective is it to run a bayesian simulation on parallel (or distrib

Re: [R] logistic regression using "glm",which "y" is set to be "1"

Dear Bin: you type ?glm in R console and you will find the Detail section of help file for glm i pasted it for you too Details A typical predictor has the form response ~ terms where response is the (numeric) response vector and terms is a series of terms which specifies a linear predictor for r

Re: [R] logistic regression using "glm",which "y" is set to be "1"

Dear Marc Schwartz: When I ask R2.6.0 for windows, the information it gives does not contain much about family=binomial . You said that there is a detail section of "?glm". I want to read it thoroughly. Could you tell me where and how I can find the detail section of "?glm". Thank you very m

Re: [R] logistic regression using "glm",which "y" is set to be "1"

On Wed, 2007-12-05 at 18:06 -0800, Bin Yue wrote: > Dear friends : > using the "glm" function and setting family=binomial, I got a list of > coefficients. > The coefficients reflect the effects of predicted variables on the > probability of the response to be "1". > My response variable consi

[R] Using expression in Hmisc Key()

Dear r-helpers, How do I tell xYplot() and Key() that I want the labels in italic? I tried this: pdf(file = "khw.pdf", width = 8, height = 8) with(subset(khw, length < 1.6), xYplot(lo ~ length|sub, groups = v, aspect="xy", pch = c(1, 2, 16), col = c(2, 4, 9), xlab=expression(frac(abs(bolditalic(

[R] logistic regression using "glm",which "y" is set to be "1"

Dear friends : using the "glm" function and setting family=binomial, I got a list of coefficients. The coefficients reflect the effects of predicted variables on the probability of the response to be "1". My response variable consists of "A" and "D" . I don't know which level of the response

Re: [R] HTML help search in R 2.6.0 v 2.6.1

Peter Dalgaard biostat.ku.dk> writes: > > > > Has something changed from 2.6.0 to 2.6.1 that may require JVM > 1.4.1? > > If > > so, I can use that information to request an upgrade of my JVM. > > > > > > > Hmm, could be. They got rebuilt on my system and committted at some > point in

Re: [R] HTML help search in R 2.6.0 v 2.6.1

Michael Bibo wrote: > I am running R on a corporate Windows XP SP2 machine on which I do not > have > administrator privileges or access to most settings in Control Panel. > R is > installed from my limited user account. The version of the JVM I have > > installed is perhaps best described as a

[R] prediction R-squared

Hi, I need to compute a prediction R-squared for a linear regression. I have figured out how to compute Allen's PRESS statistic (using the PRESS function in the MPV library), but also want to compute the R-squared that goes along with this statistic. I have read that this is computed like an ad

[R] correlation coefficient from qq plot

Hi, I am trying to figure out how to get the correlation coefficient for a QQ plot (residual plot). So to be more precise, I am creating the plot like this: qq.plot(rstudent(regrname), main = rformula, col=1) But want to also access (or compute) the correlation coefficient for that plot.

Re: [R] File based configuration

See the following (the email seems to have wrapped 2 lines onto one at one point but it should be obvious): http://tolstoy.newcastle.edu.au/R/e2/help/07/06/18853.html On Dec 5, 2007 5:05 PM, Thomas Allen <[EMAIL PROTECTED]> wrote: > I'm wanting to run R scripts non-interactively as part of a > "t

[R] hclust in heatmap.2

Dear list, I am using heatmap.2(x) to draw a heatmap. Ideally, I want to the matrix x clustered only by columns and keep the original order of rows unchanged. Is there a way to do that in heatmap.2()? Thanks a lot! Any suggestions will be appreciated! Best, Allen [[alternative HTM

[R] HTML help search in R 2.6.0 v 2.6.1

I am running R on a corporate Windows XP SP2 machine on which I do not have administrator privileges or access to most settings in Control Panel. R is installed from my limited user account. The version of the JVM I have installed is perhaps best described as antique: > system(paste("java -v

Re: [R] predict error for survreg with natural splines

Charles C. Berry wrote: > On Wed, 5 Dec 2007, Gad Abraham wrote: > >> Hi, >> >> The following error looks like a bug to me but perhaps someone can shed >> light on it: >> >> > library(splines) >> > library(survival) >> > s <- survreg(Surv(futime, fustat) ~ ns(age, knots=c(50, 60)), >> data=ovarian

Re: [R] Dealing with NA's in a data matrix

Henrique Dallazuanna gmail.com> writes: > x[is.na(x)] <- 0 > > On 05/12/2007, Amit Patel yahoo.co.uk> wrote: > > Hi I have a matrix with NA value that I would like to convert these to a value of 0. > > any suggestions also library(gdata) x <- matrix(rnorm(16), nrow=4, ncol=4) x[1, 1] <- NA NA

Re: [R] nearest correlation to polychoric

Dear Jens, I've submitted a new version (0.7-4) of the polycor package to CRAN. The hetcor() function now uses your nearcor() in sfsmisc to make the returned correlation matrix positive-definite if it is not already. I know that quite some time has elapsed since you raised this issue, and I apolo

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

Prof Brian Ripley a écrit : > Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit > version of R somewhat moot. See chapter 8 of > > http://cran.r-project.org/doc/manuals/R-admin.html Thank you for this reminder|tip ! I didn't re-read this document since ... oh, my ... very e

Re: [R] Junk or not Junk ???

Thank you As per advice from several R users I have set r-project.org, stat.math.ethz.ch,fhcrc.org, stat.ethz.ch, math.ethz.ch, hypatia.math.ethz.ch all to be "safe domains" But still some R emails go to Junk and require to be found manually I have explored the issue with Univ Wash computing

Re: [R] 2/3d interpolation from a regular grid to another regular grid

On 2007-December-05 , at 16:47 , Scionforbai wrote: >> I just read the description in ?Krig in the package fields which >> says: >> " Fits a surface to irregularly spaced data. " > > Yes, that is the most general case. Regular data location is a subset > of irregular. Anyway, kriging, just one g

Re: [R] Java parser for R data file?

David Coppit wrote: > Hi everyone, > > Has anyone written a parser in Java for either the ASCII or binary format > produced by save()? You might want to consider using the hdf5 package to save the array in HDF5 format. There are HDF5 libraries for Java as well

[R] coxme frailty model standard errors?

Hello, I am running R 2.6.1 on windows xp I am trying to fit a cox proportional hazard model with a shared Gaussian frailty term using coxme My model is specified as: nofit1<-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat) With x1-x3 being dummy variables, and isl being the communi

Re: [R] Learning to do randomized block design analysis

Dear R-Helpers (1) After a night's sleep, I realized why the other helpers think differently from me. I agree with others that it may be better to use multi-stratum model but I was a bit surprised since they seem to think 'block' variable should *not* be a fixed effect. A. Others seemed to think,

[R] File based configuration

I'm wanting to run R scripts non-interactively as part of a "technology independent" framework. I want control over the behaviour of these processes by specifying various global variables in a "configuration file" that would be passed as a command line argument. I'm wondering if you know of any R

Re: [R] Displaying numerics to full double precision

On 12/5/07, Jeff Delmerico <[EMAIL PROTECTED]> wrote: > Does anyone know if there is a way to change the display or storage settings > so that the values will be displayed to their full precision? Or does > rnorm only produce values to single precision? > > Any assistance would be greatly appreci

Re: [R] kalman filter random walk

You may want to look at package dlm. Giovanni > Date: Wed, 05 Dec 2007 12:05:00 -0600 > From: Alexander Moreno <[EMAIL PROTECTED]> > Sender: [EMAIL PROTECTED] > Precedence: list > > Hi, > > I'm trying to use the kalman filter to estimate the variable drift of a > random walk, given that I have

Re: [R] Interpretation of 'Intercept' in a 2-way factorial lm

You estimate a model with the Factors A or B either present (1) or not present (0) and with an intercept. Thus you would predict: For both A and B not present: Intercept For A only present: Intercept+coef(A) For B only preseent: Intercept+coef(B) For both present: Intercept+coef(A)+coef(B). Again

[R] coxme frailty model standard errors?

Hello, I am running R 2.6.1 on windows xp I am trying to fit a cox proportional hazard model with a shared Gaussian frailty term using coxme My model is specified as: nofit1<-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat) With x1-x3 being dummy variables, and isl being the commun

Re: [R] Displaying numerics to full double precision

Jeff Delmerico wrote: > > Thanks Ben, that fixed the display within R. However, even after changing > the display settings, the matrix elements still appear to be exported in > single precision. The matrix object is being passed into my C routines as > an SEXP Numeric type, and somewhere alon

Re: [R] newbie lapply question

I am not sure to understand your problem, but it seems to me that you can use directly the function on the range of the dates: > x=as.Date(c('2007-01-01','2007-01-02')) > fff=function(x){y=x+1;return(y)} > fff(x) [1] "2007-01-02" "2007-01-03" > class(fff(x)) [1] "Date" Perhaps your function

Re: [R] kalman filter random walk

sspir or ldm. -Roy M. On Dec 5, 2007, at 10:05 AM, Alexander Moreno wrote: > Hi, > > I'm trying to use the kalman filter to estimate the variable drift > of a > random walk, given that I have a vector of time series data. > Anyone have > any thoughts on how to do this in R? > > Thanks, > A

Re: [R] Interpretation of 'Intercept' in a 2-way factorial lm

Gustaf Granath wrote: > Hi all, > > I hope this question is not too trivial. I can't find an explanation > anywhere (Stats and R books, R-archives) so now I have to turn to the R-list. > > Question: > > If you have a factorial design with two factors (say A and B with two > levels each). What does

Re: [R] Is R portable?

On 05-Dec-07 18:57:58, Barry Rowlingson wrote: > Peter Dalgaard wrote: > >> The toolchain availability tends to get in the way. Linux-based >> gadgets could prove easier. I do wonder from time to time whether >> there really is a market for R on cellphones... > > As soon as someone writes library

Re: [R] Dealing with NA's in a data matrix

x[is.na(x)] <- 0 On 05/12/2007, Amit Patel <[EMAIL PROTECTED]> wrote: > Hi I have a matrix with NA value that I would like to convert these to a > value of 0. > any suggestions > > Kind Regards > Amit Patel > > > > > ___ > > ttp://uk.

[R] Working with "ts" objects

I am relatively new to R and object oriented programming. I have relied on SAS for most of my data analysis. I teach an introductory undergraduate forecasting course using the Diebold text and I am considering using R in addition to SAS and Eviews in the course. I work primarily with univariate or

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

On Wed, 05-Dec-2007 at 06:11PM +0100, Peter Dalgaard wrote: [] |> One oddity about Ubuntu is that there are no CRAN builds for 64bit. |> Presumably, the Debian packages work, or you can get the build |> script and make your own build. This is not really within my |> domain, though. I've alwa

[R] kalman filter random walk

Hi, I'm trying to use the kalman filter to estimate the variable drift of a random walk, given that I have a vector of time series data. Anyone have any thoughts on how to do this in R? Thanks, Alex [[alternative HTML version deleted]] __ R-h

[R] Dealing with NA's in a data matrix

Hi I have a matrix with NA value that I would like to convert these to a value of 0. any suggestions Kind Regards Amit Patel ___ ttp://uk.promotions.yahoo.com/forgood/ [[alternative HTML version deleted]] ___

[R] Displaying numerics to full double precision

I'm working on a shared library of C functions for use with R, and I want to create a matrix in R and pass it to the C routines. I know R computes and supposedly stores numerics in double precision, but when I create a matrix of random numerics using rnorm(), the values are displayed in single pr

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

>> One oddity about Ubuntu is that there are no CRAN builds for 64bit. > > Volunteers would be welcomed with open arms. > I don't think I know enough to volunteer (but I am open to suggestion off list), but I had no difficulty installing R from source on a 64 bit Ubuntu machine (7.10, with, I

Re: [R] Displaying numerics to full double precision

Jeff Delmerico wrote: > > I'm working on a shared library of C functions for use with R, and I want > to create a matrix in R and pass it to the C routines. I know R computes > and supposedly stores numerics in double precision, but when I create a > matrix of random numerics using rnorm(), th

[R] Interpretation of 'Intercept' in a 2-way factorial lm

Hi all, I hope this question is not too trivial. I can't find an explanation anywhere (Stats and R books, R-archives) so now I have to turn to the R-list. Question: If you have a factorial design with two factors (say A and B with two levels each). What does the intercept coefficient with treatm

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

On Wed, Dec 05, 2007 at 06:11:40PM +0100, Peter Dalgaard wrote: > One oddity about Ubuntu is that there are no CRAN builds for 64bit. Volunteers would be welcomed with open arms. Dirk, -- Three out of two people have difficulties with fractions. __ R

Re: [R] Is R portable?

Peter Dalgaard wrote: > The toolchain availability tends to get in the way. Linux-based gadgets > could prove easier. I do wonder from time to time whether there really > is a market for R on cellphones... As soon as someone writes library(ringtone) there might be :) And I think you'd have t

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

Prof Brian Ripley wrote: > Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit > version of R somewhat moot. See chapter 8 of > > http://cran.r-project.org/doc/manuals/R-admin.html > Only somewhat. The Opteron actually has 1GB too (Hey, it was bought in 2004! And the main point

Re: [R] R function for percentrank

On Wed, 2007-12-05 at 18:42 +0100, Martin Maechler wrote: > I'm coming late to this, but this *does* need a correction > just for the archives ! > > > "MS" == Marc Schwartz <[EMAIL PROTECTED]> > > on Sat, 01 Dec 2007 13:33:21 -0600 writes: > > MS> On Sat, 2007-12-01 at 18:40 +000

Re: [R] Is R portable?

michael watson (IAH-C) wrote: > I opened this hoping someone had installed R on windows mobile or simbian > :( > I think that's Symbian with a 'y'. The toolchain availability tends to get in the way. Linux-based gadgets could prove easier. I do wonder from time to time whether there really

Re: [R] confidence intervals for y predicted in non linear regression

On Wed, 5 Dec 2007, Peter Dalgaard wrote: > Prof Brian Ripley wrote: >> You mean the package nls2 at >> >> http://w3.jouy.inra.fr/unites/miaj/public/AB/nls2/install.html >> >> and not the unfortunately named nls2 that has just appeared on CRAN? >> >> The first is not really a 'package' in the R se

Re: [R] how to interpolate a plot with a logistic curve

On Wednesday 05 December 2007, Simone Gabbriellini wrote: > hello, > > I have this simple question. This is my dataset > > size > 1 57 > 2 97 > 3 105 > 4 123 > 5 136 > 6 153 > 7 173 > 8 180 > 9 193 > 10202 > 11213 > 12219 > 13224 > 14224

[R] how to interpolate a plot with a logistic curve

hello, I have this simple question. This is my dataset size 1 57 2 97 3 105 4 123 5 136 6 153 7 173 8 180 9 193 10 202 11 213 12 219 13 224 14 224 15 248 16 367 17 496 18 568 19 618 20

Re: [R] confidence intervals for y predicted in non linear regression

Prof Brian Ripley wrote: > You mean the package nls2 at > > http://w3.jouy.inra.fr/unites/miaj/public/AB/nls2/install.html > > and not the unfortunately named nls2 that has just appeared on CRAN? > > The first is not really a 'package' in the R sense. > Actually, it is one (sort of), but it is brok

Re: [R] R function for percentrank

I'm coming late to this, but this *does* need a correction just for the archives ! > "MS" == Marc Schwartz <[EMAIL PROTECTED]> > on Sat, 01 Dec 2007 13:33:21 -0600 writes: MS> On Sat, 2007-12-01 at 18:40 +, David Winsemius wrote: >> David Winsemius <[EMAIL PROTECTED]> wrot

Re: [R] newbie lapply question

On Wed, 5 Dec 2007, Prof Brian Ripley wrote: >[...] Thanks I'll read it more carefully. > Perhaps if you told us what you are trying to achieve we might be able to > help you achieve it. > I have a function which takes a date as an argument. I've tested it, and I'd like to run it over a ran

Re: [R] newbie lapply question

On Wed, 5 Dec 2007, Ranjan Bagchi wrote: > Hi -- > > I just noticed the following (R 2.6.1 on OSX) > >> lapply(c(as.Date('2007-01-01')), I) > [[1]] > [1] 13514 > > This is a bit surprising.. Why does lapply unclass the object? Sorry for > such a basic question, I don't seem able to produce the ri

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit version of R somewhat moot. See chapter 8 of http://cran.r-project.org/doc/manuals/R-admin.html I would install a i386 version of R on x86_64 Linux unless I had 2Gb or more of RAM. I don't know how easily that works on Ub

Re: [R] Java parser for R data file?

On Wed, 5 Dec 2007, David Coppit wrote: > Hi everyone, > > Has anyone written a parser in Java for either the ASCII or binary format > produced by save()? I need to parse a single large 2D array that is > structured like this: > > list( > "32609_1" = c(-9549.39231289146, -9574.07159324482, ... ),

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

8rino-Luca Pantani wrote: > Dear R-users. > I eventually bought myself a new computer with the following > characteristics: > > Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2) > Mother board ASR SK-AM2 2 > Ram Corsair Value 1 GB DDR2 800 Mhz > Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB > > I'

[R] newbie lapply question

Hi -- I just noticed the following (R 2.6.1 on OSX) > lapply(c(as.Date('2007-01-01')), I) [[1]] [1] 13514 This is a bit surprising.. Why does lapply unclass the object? Sorry for such a basic question, I don't seem able to produce the right google keywords. Ranjan ___

[R] Which Linux OS on Athlon amd64, to comfortably run R?

Hi Ottorino, I have been using R on 64-bit Ubuntu for about a year without problems, both Intel and AMD CPUs. Installing and using several packages (e1071, svmpath, survival) also works. However, I had to install R from source: $ gunzip -c R-2.6.1.tar.gz | tar xvf - $ cd R-2.6.1 $ ./configure --

Re: [R] Calculating large determinants

Hi Murray, A likely reason for the observed information matrix not to be positive definite is the inaccuracies in the numerical estimation of the scores. You might want to try more accurate methods (e.g. Richardson extrapolation) for approximating numerical derivatives, such as available in the p

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

> but unfortunately it did not work. What did not work? Provide some information more... By the way, isn't 'gfortran' the new GNU fortran compiler which replaced 'g77'? Or not on Ubuntu? > I immediately realize that an amd64 is not the right processor to make > life easy. ??? Example of bad extr

[R] Java parser for R data file?

Hi everyone, Has anyone written a parser in Java for either the ASCII or binary format produced by save()? I need to parse a single large 2D array that is structured like this: list( "32609_1" = c(-9549.39231289146, -9574.07159324482, ... ), "32610_2" = c(-6369.12526971635, -6403.99620977124,

Re: [R] Plotting error bars in xy-direction

Hans W Borchers wrote: > > Dear R-help, > > I am looking for a function that will plot error bars in x- or y-direction > (or > both), the same as the Gnuplot function 'plot' can achieve with: > > plot "file.dat" with xyerrorbars,... > > Rsite-searching led me to the functions 'errbar' a

Re: [R] 2/3d interpolation from a regular grid to another regular grid

> I just read the description in ?Krig in the package fields which says: > " Fits a surface to irregularly spaced data. " Yes, that is the most general case. Regular data location is a subset of irregular. Anyway, kriging, just one g, after the name of Danie Krige, the south african statistician w

Re: [R] Which Linux OS on Athlon amd64, to comfortably run R?

8rino-Luca Pantani wrote: >Dear R-users. >I eventually bought myself a new computer with the following >characteristics: > >Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2) >Mother board ASR SK-AM2 2 >Ram Corsair Value 1 GB DDR2 800 Mhz >Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB > >I'm a new

Re: [R] Multiple stacked barplots on the same graph?

And qplot(x=Categorie,y=Total,data=mydata,geom="bar",fill=Part) + coord_flip() makes it a bit easier to read the labels. Hadley On Dec 5, 2007 8:33 AM, Domenico Vistocco <[EMAIL PROTECTED]> wrote: > This command works: > > qplot(x=Categorie,y=Total,data=mydata,geom="bar",fill=Part) > > for your

[R] Which Linux OS on Athlon amd64, to comfortably run R?

Dear R-users. I eventually bought myself a new computer with the following characteristics: Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2) Mother board ASR SK-AM2 2 Ram Corsair Value 1 GB DDR2 800 Mhz Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB I'm a newcomer to the Linux world. I started u

[R] Asymmetrically dependent variables to 2D-map?

Hello, I'm searching for a method which maps variables of this kind of table, see below, to 2-dimensional space, like in multidimensional scaling. However, this table is asymmetric: for example, variable T1 affects T2 more than T2 affects T1(0.41 vs. 0.21). > DEPTABLE T1T2 T3

[R] how to estimate ARCH?

Dear R gurus, Could somebody enlighten me on how to estimate the parameters of ARCH process? I tried both garchFit and garch functions, but no one works for me. R version 2.6.1, Windows library(fGarch) set.seed(5) # simulate ARCH(3) rr=garchSim(model = list(alpha = c(0.1,0.2,0.6)),n=400) rrm=rr

Re: [R] Multiple stacked barplots on the same graph?

This command works: qplot(x=Categorie,y=Total,data=mydata,geom="bar",fill=Part) for your data. domenico vistocco Stéphane CRUVEILLER wrote: > Hi, > > the same error message is displayed with geom="bar" as parameter. > here is the output of dput: > > > dput(mydata) > structure(list(Categorie = s

Re: [R] Calculating large determinants

Hi Murray, A likely reason for the observed information matrix not to be positive definite is the inaccuracies in the numerical estimation of the scores. You might want to try more accurate methods (e.g. Richardson extrapolation) for approximating numerical derivatives, such as available in the p

[R] Significance of clarkevans in spatstat

Dear R-Users, I was wondering if there is a way to test the significance of the clarkevans statistic in spatstat package? I did not find any related function or the related values to calculate it by hand. does someone has any ideas? thank you, Jens -- +

Re: [R] about color palettes, colorRamp etc

[if you get this twice: it seems to have not made it through, yesterday] > "Earl" == Earl F Glynn <[EMAIL PROTECTED]> > on Mon, 3 Dec 2007 13:26:11 -0600 writes: Earl> "affy snp" <[EMAIL PROTECTED]> wrote in message Earl> news:[EMAIL PROTECTED] >> For example, it sh

[R] Mixed Integer Non Linear Programming.

Hi, I am attempting to solve a Mixed Integer Non Linear Program using R. I am aware of the GLPK package that solved mixed integer linear programs, and the Rdonlp2 package that solve linear programs. Is there any package that that solves the MINLP case? Any info or pointers would be greatly apprec

[R] Export to LaTeX using xtable() - Control the digits to the right of the separator [solved]

Hello everyone, The thread title speaks for itself. Here's the code that worked for me: > numSummary(finance[,"Employees"], statistics=c("mean", "sd", "quantiles")) mean sd 0% 25% 50% 75% 100% n NA 11492.92 29373.14 1777 3040 4267 6553 179774 53 5 > str(numSummary(finance[,"

Re: [R] weighted Cox proportional hazards regression

> I'm getting unexpected results from the coxph function when using > weights from counter-matching. For example, the following code > produces a parameter estimate of -1.59 where I expect 0.63: I agree with Thomas' answer wrt using offset instead of weight. One way to understand this is to l

Re: [R] Help installing taskPR

I suggest you ask the maintainer (as does the posting guide). Red Hat headers for mpi are in a different layout from the sources, and recent versions of mpi from earlier ones. But the path to mpi.h (which I guess is /usr/local/lam-7.1.4/include in your setup) is missing. You might try adding

Re: [R] rgdal for R.2.4.0?

On Wed, 5 Dec 2007, bernardo lagos alvarez wrote: > Hi, > > Know anyone where to find the package rgdal for R.2.4.0? On CRAN: the current version has Package: rgdal Title: Bindings for the Geospatial Data Abstraction Library Version: 0.5-20 Date: 2007-11-07 Depends: R (>= 2.3.0), methods, sp Or

Re: [R] confidence intervals for y predicted in non linear regression

You mean the package nls2 at http://w3.jouy.inra.fr/unites/miaj/public/AB/nls2/install.html and not the unfortunately named nls2 that has just appeared on CRAN? The first is not really a 'package' in the R sense. On Wed, 5 Dec 2007, Ndoye Souleymane wrote: Hi, Salut, You should use the pa

[R] rgdal for R.2.4.0?

Hi, Know anyone where to find the package rgdal for R.2.4.0? I would appreciate advice very much. Best wishes, Bernardo Lagos Department of Statistics University of Concepción __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] confidence intervals for y predicted in non linear regression

I don't think this is referring to the nls2 package on CRAN but rather something else. On Dec 5, 2007 7:37 AM, Ndoye Souleymane <[EMAIL PROTECTED]> wrote: > > Hi, Salut, > > You should use the package nsl2 (only for Linux distribution) > Vous pouvez utiliser le package nls2 (Linux seulement) > > R

Re: [R] alternatives to latex() or xtable() ?

On 12/5/07, Romain Francois <[EMAIL PROTECTED]> wrote: > > > > Hello, > > My guess is that you are actually talking about the numSummary function in > Rcmdr, not in abind. In that case, you can look how the structure of the > output is like: > > > str( numSummary( iris[,1:4] ) ) > List of 4 >

Re: [R] confidence intervals for y predicted in non linearregression

Hi Thanks for your suggestion, I'm trying to install this package in Ubuntu (7.10) but unsuccessfully. Also tried in MacOSX, and no success too. _ De: Ndoye Souleymane [mailto:[EMAIL PROTECTED] Enviado el: miércoles, 05 de diciembre de 2007 13:38 Para: [EMAIL PROTECTED]; Florencio

Re: [R] Quadratic programming

G'day Serge, On Wed, 5 Dec 2007 11:25:41 +0100 "de Gosson de Varennes Serge (4100)" <[EMAIL PROTECTED]> wrote: > I am using the quadprog package and its solve.QP routine to solve and > quadratic programming problem with inconsistent constraints, which > obviously doesn't work since the constraint

Re: [R] confidence intervals for y predicted in non linear regression

Hi, Salut, You should use the package nsl2 (only for Linux distribution) Vous pouvez utiliser le package nls2 (Linux seulement) Regards, Souleymane> Date: Tue, 4 Dec 2007 16:07:57 +0100> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]> CC: [EMAIL PROTECTED]> Subject: Re: [R] confidence inter

[R] alternatives to latex() or xtable() ?

Hello everyone, I have several problems with exporting to LaTeX the output of numSummary() from the abind package. > numSummary(finance[,"Cash_flow"], statistics=c("mean", "sd", "quantiles")) mean sd 0% 25% 50% 75%100% n NA 188070.9 414771.9 -426804 26743 53866 1509

Re: [R] seasonal time serie with missing values

Hello All, I don't know if help, but I find a package written in Java and .NET that have a implementation of tramo/seats. The link http://www.nbb.be/app/dqrd/index.htm If I have more information I will post it. Thanks for the reply s, they really help! João Santos CCordeiro wrote: > > Tha

[R] Help installing taskPR

Hello, I'd like to take a look at how the taskPR package performs on our computational cluster, and this morning I've spend some time trying to do a test install of this package on one of our linux machines. Unfortunately I've been unsuccessful in completing the install, and I wondered if someo

Re: [R] Information criteria for kmeans

This is not primarily an R question: if you tell us how you want to define it, we may be able to help you compute it. I presume you are talking about Schwarz (1978), which is not billed as an 'information criterion'. AFAIK, all Gideon Schwarz did was to define a criterion for linear regression

Re: [R] Is R portable?

I opened this hoping someone had installed R on windows mobile or simbian :( From: [EMAIL PROTECTED] on behalf of Duncan Murdoch Sent: Tue 04/12/2007 10:54 PM To: Tom Backer Johnsen Cc: [EMAIL PROTECTED] Subject: Re: [R] Is R portable? On 04/12/2007 5:25 PM

Re: [R] 2/3d interpolation from a regular grid to another regular grid

On 2007-December-04 , at 21:38 , Scionforbai wrote: >> - krigging in package fields, which also requires irregular spaced >> data > > That kriging requires irregularly spaced data sounds new to me ;) It > cannot be, you misread something (I feel free to say that even if I > never used that packa

[R] Information criteria for kmeans

Hello, how is, for example, the Schwarz criterion is defined for kmeans? It should be something like: k <- 2 vars <- 4 nobs <- 100 dat <- rbind(matrix(rnorm(nobs, sd = 0.3), ncol = vars), matrix(rnorm(nobs, mean = 1, sd = 0.3), ncol = vars)) colnames(dat) <- paste("var",1:4) (cl <-

[R] lme output

Dear all, I noticed the following in the call of lme using msVerbose. fm1 <- lme(distance ~ age, data = Orthodont, control = lmeControl(msVerbose=T)) 9 318.073: -0.567886 0.152479 1.98021 10 318.073: -0.567191 0.152472 1.98009 11 318.073: -0.567208 0.152473

Re: [R] converting factors to dummy variables

Try this also: table(cbind.data.frame(Price=my.dataset$Price, Colour=paste(my.dataset$Colour, my.dataset$Store, sep=":"))) On 05/12/2007, Tim Calkins <[EMAIL PROTECTED]> wrote: > Hi all - > > I'm trying to find a way to create dummy variables from factors in a > regression. I have been using big

[R] Quadratic programming

Hi, I'm quite new at R and I haven't found the answer to my question anywhere on the net, so either it is trivial or not documented. So, bare with be. I am using the quadprog package and its solve.QP routine to solve and quadratic programming problem with inconsistent constraints, which obviousl

  1   2   >