Re: [R] how to draw a spherical quadrant

2010-08-10 Thread [Ricardo Rodriguez] Your XEN ICT Team
Thanks for pointing me on that direction! Extremely feature-rich this, as usual, package! Duncan Murdoch wrote: The rgl package can do general 3D plotting, but I don't know of anyone who has put together what you're looking for, so you'll need to compute the individual line segments or polyg

[R] question about bayesian model selection for quantile regression

2010-08-10 Thread Xin__ Li
Hi All: Recently I am researching my dissertation about the quantile model selection by bayesian approach. I have the dependent variable(return) and 16 independent variables and I need to select the best variable for each quantile of return. And the method I used is the bayesian approach, which is

Re: [R] Multiple imputation, especially in rms/Hmisc packages

2010-08-10 Thread Mark Seeto
Hi Frank (and others), Thank you for your comments in reply to my questions. I had not encountered contrast tests before. I've looked in a few texts, but the only place I could find anything about contrast tests was your Regression Modeling Strategies book. You wrote that the "leave some var

[R] Bigmemory: Error Running Example

2010-08-10 Thread harsh yadav
Hi, I am trying to run the bigmemory example provided on the http://www.bigmemory.org/ The example runs on the "airline data" and generates summary of the csv files:- library(bigmemory) library(biganalytics) x <- read.big.matrix("2005.csv", type="integer", header=TRUE, backingfile="airline.bin",

[R] Variance inflation factor

2010-08-10 Thread Grant Gillis
Hello all and thanks in advance for any advice. I would like to calculate the variance inflation factor for a linear model (lm) with 4 explanatory variables. I would then like to use this to calculate QAIC. I have used the function vif() in the car package and I get values for each variable howe

Re: [R] svydesign syntax and deviance!

2010-08-10 Thread Thomas Lumley
On Sat, 31 Jul 2010, aline uwimana wrote: Thank you Thomas,  if i understand what you say below i need to use glm in order to have the deviance and consider it for the model with svyglm and choose the best model. I use  glm and svyglm for two models but don't have deviance for the svyglm. As i

Re: [R] Intersecting list vs rows in matrix

2010-08-10 Thread GL
Very cool. Thanks! -Original Message- From: "Henrique Dallazuanna [via R]" To: Lipori, Gigi Sent: 08/10/2010 05:18:25 PM Subject: Re: Intersecting list vs rows in matrix Try this: colSums(apply(List_2, 1, is.element, List_1)) > 0 On Tue, Aug 10, 2010 at 5:42 PM, GL wrote: >

[R] TRUE FALSE issue

2010-08-10 Thread Alexander Eggel
How can I extract the samples (S1-S5) containing a TRUE value in their row? Solution should apply to a much bigger data frame. > a Samples A B C D . . . 1 S1 FALSE FALSE FALSE FALSE 2 S2 FALSE FALSE NA TRUE 3 S3 FALSE FALSE FALSE FALSE 4 S4 FALSE TRUE FALSE FALSE 5 S5 FALSE FALSE FALSE FALSE . . .

Re: [R] TRUE FALSE issue

2010-08-10 Thread Erik Iverson
On 08/10/2010 09:04 PM, Alexander Eggel wrote: How can I extract the samples (S1-S5) containing a TRUE value in their row? Solution should apply to a much bigger data frame. a Samples A B C D . . . 1 S1 FALSE FALSE FALSE FALSE 2 S2 FALSE FALSE NA TRUE 3 S3 FALSE FALSE FALSE FALSE 4 S4 FALSE TR

Re: [R] TRUE FALSE issue

2010-08-10 Thread David Winsemius
On Aug 10, 2010, at 10:16 PM, Erik Iverson wrote: On 08/10/2010 09:04 PM, Alexander Eggel wrote: How can I extract the samples (S1-S5) containing a TRUE value in their row? Solution should apply to a much bigger data frame. a Samples A B C D . . . 1 S1 FALSE FALSE FALSE FALSE 2 S2 FALSE F

Re: [R] TRUE FALSE issue

2010-08-10 Thread Bill.Venables
> with(a, Samples[apply(a[,-1], 1, any)]) [1] S2 S4 Levels: S1 S2 S3 S4 S5 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Alexander Eggel Sent: Wednesday, 11 August 2010 12:05 PM To: r-help@r-project.org Subject: [R] TRUE FALSE is

Re: [R] TRUE FALSE issue

2010-08-10 Thread Bill.Venables
For huge cases this might be a whiff faster > with(a, Samples[rowSums(a[, -1], na.rm = TRUE) > 0]) [1] S2 S4 Levels: S1 S2 S3 S4 S5 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of bill.venab...@csiro.au Sent: Wednesday, 11 August

[R] Where to download files for ANOVA test?

2010-08-10 Thread Stephen Liu
Hi folks, I'm learning ANOVA needing sample data files for the test. > data() displays many data files. But I have no idea which of them is suitable for my application. I also found following thread; R and Analysis of Variance http://www.personality-project.org/R/r.anova.html The sample d

Re: [R] Where to download files for ANOVA test?

2010-08-10 Thread David Winsemius
On Aug 10, 2010, at 10:28 PM, Stephen Liu wrote: Hi folks, I'm learning ANOVA needing sample data files for the test. data() displays many data files. But I have no idea which of them is suitable for my application. I also found following thread; R and Analysis of Variance http://www.

Re: [R] Where to download files for ANOVA test?

2010-08-10 Thread Stephen Liu
> - Original Message > From: David Winsemius > To: Stephen Liu >Cc: r-help@r-project.org >Sent: Wed, August 11, 2010 11:42:55 AM >Subject: Re: [R] Where to download files for ANOVA test? - snip - >> >> Please advise where can I find a bigger size sample files >> for my test? > http:

Re: [R] Where to download files for ANOVA test?

2010-08-10 Thread David Winsemius
On Aug 11, 2010, at 12:36 AM, Stephen Liu wrote: - Original Message From: David Winsemius To: Stephen Liu Cc: r-help@r-project.org Sent: Wed, August 11, 2010 11:42:55 AM Subject: Re: [R] Where to download files for ANOVA test? - snip - Please advise where can I find a bigger

[R] sweep and zoo objects

2010-08-10 Thread steven mosher
rc<-list(c( 123,321,234,543,654,768,986,987,246,284),c("Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec")) # the matrix has rownames that are used as identifiers and columns # of time. 1 years worth of data. Thats the native format test<-matrix(seq(1,120, by=1), nrow=10,dim

Re: [R] package for measurement error models

2010-08-10 Thread Prof Brian Ripley
On Tue, 10 Aug 2010, Carrie Li wrote: Thank you both! I found that the model with both x and y have measurement errors should be pretty common in practice. But it seems to me that there is no a simple solution for it...(I mean, a ready-to-use package or program handing this model fitting proble

[R] Error: the leading minor of order 6 is not positive definite

2010-08-10 Thread stompper33
Hey guys, I'm trying to run a Canonical Correlation Analysis (CCA) between two data sets. But for my case I am dealing with data sets(which are stored as two matrices X and Y) in which the number of experimental units is greater than the number of variables, so I want to use a sample from my "mot

Re: [R] [BUGS] [R-BUGS] error while plotting

2010-08-10 Thread Anamika Chaudhuri
Trevor, Thanks for your reply. That doesnot help Any other suggestions? Anamika On Tue, Aug 10, 2010 at 11:58 AM, Trevor Davies wrote: > I think you need to load is R2WinBUGS again. > > require(R2WinBUGS) > > Trevor > > > Hi All: > > > > I am getting an error while trying to plot in R2Winbugs >

[R] problem with Bitmap

2010-08-10 Thread kayj
Hi All, I am trying to run the following script but I am getting and error message one=read.table("sample.txt",sep="\t") bitmap(file="sample.JPG",type="jpeg",width=5,height=5,res=300,pointsize=10) “Error in (st + 1):(en - 1) : argument of length 0” plot(one$V1,one$V2) I tried to google the error

[R] Std. error of correlation coefficients

2010-08-10 Thread fusun gonul
How do I get the std. errors of correlation coefficients? When I use the cor function it only prints the correlation. Thanks, Fusun. -- View this message in context: http://r.789695.n4.nabble.com/Std-error-of-correlation-coefficients-tp2320609p2320609.html Sent from the R help mailing list arch

[R] axis labels defaulting to scientific notation

2010-08-10 Thread Watkins, Janice
The labels on the x-axis are defaulting to scientific notation no matter how small cex.axis is. How can I override scientific notation to get the labels to print out as specified? Here is the code (UNIT here is 0.0105): plot(xm,yv,log="xy",ylim=c(0.1,20)/UNIT,xlim=c(0.004,20)*UNIT,xaxt="n",t ype=

<    1   2