Re: [R] install error uder Cygwin

2006-09-23 Thread Prof Brian Ripley
On Fri, 22 Sep 2006, Duncan Murdoch wrote: On 9/22/2006 10:29 PM, X.H Chen wrote: Dear R users, I have a question about R installation under Cygwin. When running ./configure, I can't pass the checking phase due to an error: --with-readline=yes (default) and headers/libs are not available...

Re: [R] logistic + neg binomial + ...

2006-09-23 Thread Prof Brian Ripley
On Fri, 22 Sep 2006, [EMAIL PROTECTED] wrote: Hi Folks, I've just come across a kind of problem which leads me to wonder how to approach it in R. Basically, each a set of items is subjected to a series of impacts until it eventually fails. The force of each impact would depend on

Re: [R] Update to Dillo browser question

2006-09-23 Thread Rainer M Krug
I start the html help with help.start(browser=dillo) and all subsequent requests with ? open in dillo Rainer Dylan Beaudette wrote: hi, how exactly do you get R to open the page in dillo ... ? thanks! On Thursday 21 September 2006 23:54, Rainer M Krug wrote: Hi I asked about

Re: [R] nls

2006-09-23 Thread Spencer Graves
I used debug to walk through your example line by line, I found that the error message was misleading. By making as.vector(semivariance) and as.vector(h) columns of a data.frame, I got it to work. My revised code appears below. Thanks for providing a self-contained and

[R] recursive question

2006-09-23 Thread X.H Chen
Hi all, How can store the recursive results from each step without using global operator -? Thanks ahead. Xiaohui Chen Dept. of Statistics UBC, Canada _ Don’t waste time standing in line—try shopping online. Visit Sympatico /

[R] frailty in coxph

2006-09-23 Thread Mohammad Ehsanul Karim
Dear all, I have been doing some frailty calculations and been facing some difficulties. I can extract coefficients, value of theta and the following things library(survival) fit-coxph(Surv(time,status)~covariate+frailty(group), data=simulated.data) fit$coef fit$history[[1]]$theta

Re: [R] recursive question

2006-09-23 Thread Peter Dalgaard
X.H Chen [EMAIL PROTECTED] writes: Hi all, How can store the recursive results from each step without using global operator -? Thanks ahead. By returning them to the caller, recursively... -- O__ Peter Dalgaard Øster Farimagsgade 5, Entr.B c/ /'_ --- Dept. of

Re: [R] Lattice strip labels for two factors

2006-09-23 Thread Rafael Duarte
Thank you for your suggestion. This could be a solution that I didn't think of. But I forgot to say that I didn't want to change the original data frame (I have other code that depends on the original df and on the original factor levels). I was looking more for an implementation directly in

Re: [R] Fatal error: unable to restore saved data in .RData --- no package called 'nlme'

2006-09-23 Thread Prof Brian Ripley
On Fri, 22 Sep 2006, michael papenfus wrote: I am using Windows XP and R 2.3.1. During my lastest session I updated my packages and now when I try to start R 2.3.1 I get the following error message: Fatal error: unable to restore saved data in .RData in an error window and Error in

Re: [R] Y-axis on pbinom

2006-09-23 Thread Maciej Bliziński
On Sat, 2006-09-23 at 01:37 -0400, Ethan Johnsons wrote: How do you change the Y-axis from 0 ~ 0.6? plot(..., ylim = c(0, 0.6)) -- Maciej Bliziński http://automatthias.wordpress.com __ R-help@stat.math.ethz.ch mailing list

Re: [R] proj4R library will not install

2006-09-23 Thread Roger Bivand
On Fri, 22 Sep 2006, Philip Bermingham wrote: proj4R was never a CRAN package. Its functionality is included in the CRAN package rgdal, please use that instead. I'm hoping someone can help me. I have downloaded the proj4R.zip and under my version of R (2.3.1) I install the package from

[R] Fitdistr() versus nls()

2006-09-23 Thread Luca Telloli
Hello R-Users, I'm new to R so I apologize in advance for any big mistake I might be doing. I'm trying to fit a set of samples with some probabilistic curve, and I have an important question to ask; in particular I have some data, from which I calculate manually the CDF, and then I

Re: [R] nls

2006-09-23 Thread Martin Maechler
SpG == Spencer Graves [EMAIL PROTECTED] on Sat, 23 Sep 2006 00:52:30 -0700 writes: SpG I used debug to walk through your example line by line, I found SpG that the error message was misleading. By making SpG as.vector(semivariance) and as.vector(h) columns of a data.frame,

Re: [R] logistic + neg binomial + ...

2006-09-23 Thread Ted Harding
On 23-Sep-06 Prof Brian Ripley wrote: On Fri, 22 Sep 2006, [EMAIL PROTECTED] wrote: Hi Folks, I've just come across a kind of problem which leads me to wonder how to approach it in R. Basically, each a set of items is subjected to a series of impacts until it eventually fails. The force

[R] Create a vector of indices from a matrix of start and end points

2006-09-23 Thread George W. Gilchrist
I have a very large dataframe and wish to extract a subset of rows. I have a two column matrix listing the starting and ending indices of one subset on each row. My idea is to create a vector of indices that could be applied to the dataframe and I have a solution using a for loop (below).

Re: [R] Create a vector of indices from a matrix of start and end points

2006-09-23 Thread jim holtman
try this: x - rbind(c(2,5), c(7,9), c(15,20)) x [,1] [,2] [1,]25 [2,]79 [3,] 15 20 unlist(mapply(seq, x[,1], x[,2])) [1] 2 3 4 5 7 8 9 15 16 17 18 19 20 On 9/23/06, George W. Gilchrist [EMAIL PROTECTED] wrote: I have a very large dataframe and wish to

Re: [R] Create a vector of indices from a matrix of start and end points

2006-09-23 Thread David Barron
unlist(apply(tmp1,1,function(x) seq(x[1],x[2]))) On 23/09/06, George W. Gilchrist [EMAIL PROTECTED] wrote: I have a very large dataframe and wish to extract a subset of rows. I have a two column matrix listing the starting and ending indices of one subset on each row. My idea is to create a

[R] problem with summary of data

2006-09-23 Thread Nel N.
Greetings everybody, I'm new with R and I'm trying to do some statistical analysis with R using my data. It is a DNA microarray data with 69 samples(row) and 7070 (column) variables. I read the file using the command 'file.table' and it can be done with no error.However, when I tried running a

Re: [R] Lattice strip labels for two factors

2006-09-23 Thread Gabor Grothendieck
1. You can write a custom strip function: my.strip - function(which.given, ..., factor.levels) { levs - if (which.given == 1) factor.levels else c(faro, porto, lisbon, setubal) strip.default(which.given, ..., factor.levels = levs) }

Re: [R] [BioC] two questions associated with heatmap

2006-09-23 Thread Weiwei Shi
oh, i mean how it defines the group? i just want to confirm if it is based on the whole matrix? On 9/23/06, Sean Davis [EMAIL PROTECTED] wrote: Sean Davis wrote: Weiwei Shi wrote: hi, there: i have 2 questions associated with heatmap in heatmap.2{gplot}, there is a bar called raw

Re: [R] [BioC] two questions associated with heatmap

2006-09-23 Thread Sean Davis
Weiwei Shi wrote: oh, i mean how it defines the group? i just want to confirm if it is based on the whole matrix? The process is typically done within each probe, not across the entire matrix, but it depends a bit on the parameters given, at least for heatmap and heatmap.2. Reading the help

Re: [R] logistic + neg binomial + ...

2006-09-23 Thread Ted Harding
On 22-Sep-06 Ted Harding wrote: I've just come across a kind of problem which leads me to wonder how to approach it in R. Basically, each a set of items is subjected to a series of impacts until it eventually fails. The force of each impact would depend on covariates X,Y say; [...] ...

Re: [R] Y-axis on pbinom

2006-09-23 Thread Ethan Johnsons
thx so much..Maciej . ej On 9/23/06, Maciej Bliziński [EMAIL PROTECTED] wrote: On Sat, 2006-09-23 at 01:37 -0400, Ethan Johnsons wrote: How do you change the Y-axis from 0 ~ 0.6? plot(..., ylim = c(0, 0.6)) -- Maciej Bliziński http://automatthias.wordpress.com

Re: [R] Creating Movies with R

2006-09-23 Thread Henrik Bengtsson
On 9/22/06, Gabor Grothendieck [EMAIL PROTECTED] wrote: See the flag= argument on formatC: n - 10 formatC(1:n, dig = nchar(n)-1, flag = 0) # Here is another way: n - 10 sprintf(paste(%0, nchar(n), .0f, sep = ), 1:n) sprintf(%0*.0f, nchar(n), 1:n) or even sprintf(%0*d, nchar(n), 1:n)

[R] pbinom with rnorm

2006-09-23 Thread Ethan Johnsons
Is there a way to plug in the random samples, rnorm(100), in this plot? I am using 100, but want it to be randome samples. plot(seq(from=0,to=100,by=1),1-pbinom (seq(from=0,to=100,by=1),size=100,prob=0.05),pch=15) thx much __ R-help@stat.math.ethz.ch

Re: [R] problem with summary of data

2006-09-23 Thread David Barron
You haven't shown the command you used to input the data, but it looks as though you haven't told R not to convert the input into factors. summary on factors shows the frequencies of the levels for each variable, which is what your output shows. See ?read.table for how to prevent conversion into

Re: [R] problem with summary of data

2006-09-23 Thread Peter Dalgaard
David Barron [EMAIL PROTECTED] writes: You haven't shown the command you used to input the data, but it looks as though you haven't told R not to convert the input into factors. Looks even more like R wasn't told that the first line of data was the variable names. summary on factors shows

Re: [R] Lattice strip labels for two factors

2006-09-23 Thread Deepayan Sarkar
On 9/23/06, Gabor Grothendieck [EMAIL PROTECTED] wrote: 1. You can write a custom strip function: my.strip - function(which.given, ..., factor.levels) { levs - if (which.given == 1) factor.levels else c(faro, porto, lisbon, setubal)

Re: [R] Adding percentage to Pie Charts

2006-09-23 Thread Jim Lemon
Hi all, Anupam Tyagi mentioned an interesting idea a few days ago. A modification in a pie chart that draws overlapping areas with a common start point at the top of the circle, can make is more informative than a dot-chart. Something like: * Start drawing at the top of the circle, as zero

[R] Data frames questions

2006-09-23 Thread Jonathan Greenberg
Hi there, couple of questions on data frames: 1) Is there a way to build an empty data frame, containing nothing but the data frame variable names? 2) Is there a way to reorder the variables in a data frame, e.g. When I go to write out a data frame using write.table or write.matrix, I want the

[R] Question about merge()

2006-09-23 Thread Jonathan Greenberg
If I want to do a join based on *two* matching fields in two data frames, can merge() handle this? It appears to only handle a single matching column -- do I need to make a metacolumn or is there some way to do this? E.g.: Dataframe 1 contains columns A,B,C and Dataframe 2 contains A,B,D I

[R] plotting grouped data object

2006-09-23 Thread Afshartous, David
All, I'd like to plot the main relationship of a grouped data object for all levels of a factor in a single panel. The sample code below creates a separate panel for each level of the factor. I realize that this could be done in other ways, but I'd like to do it via plotting the grouped data

Re: [R] Data frames questions

2006-09-23 Thread X.H Chen
1) Is there a way to build an empty data frame, containing nothing but the data frame variable names? Yes, you can do it one way as followings: df1-as.data.frame(matrix(nrow=2,ncol=2),row.names=c(R1,R2)) 2) Is there a way to reorder the variables in a data frame, e.g. When I go to write out a

Re: [R] plotting grouped data object

2006-09-23 Thread X.H Chen
I don't get your meaning in what is in datagrp, anyway, try: X11() par(new=T) before calling: plot(data.grp, outer = ~ y) Xiaohui Chen Dept. of Statistics UBC, Canada From: Afshartous, David [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] plotting grouped data object Date: Sat,

[R] really irregular time series

2006-09-23 Thread phil henshaw
I built some reasonably successful tools in a graphical database for reconstructing the developmental turning points for feedback loop driven natural processes.I'm trying to move it to R and am having difficulty with the very basics, i.e. a) defining time series using time as a natural number

Re: [R] Question about merge()

2006-09-23 Thread X.H Chen
merge() can handle this if you specify the by parameters to the same length vector for x and y. The parameters corresponds the columns you want to filter out under certain conditions and add the values to the merging data.frame. For example(modified from No.1 example from ?merge): authors -

[R] strange warning message

2006-09-23 Thread Tong Wang
Hi Everyone, I recently start to get this warning message, while loading files in to R. Could someone tell me what does it mean ? I am using R 2.3.0 with Emacs on WinXP. use of NULL environment is deprecated Thanks a lot for any help best

[R] variance-covariance structure of random effects in lme

2006-09-23 Thread Mi, Deming
Dear R users, I have a question about the patterned variance-covariance structure for the random effects in linear mixed effect model. I am reading section 4.2.2 of Mixed-Effects Models in S and S-Plus by Jose Pinheiro and Douglas Bates. There is an example of defining a compound symmetry

Re: [R] really irregular time series

2006-09-23 Thread Gabor Grothendieck
The zoo package supports irregular time series with arbitrary time classes. It include na.approx and na.spline for interpolation. library(zoo) vignette(zoo) vignette(zoo-quickref) On 9/23/06, phil henshaw [EMAIL PROTECTED] wrote: I built some reasonably successful tools in a graphical

Re: [R] Question about merge()

2006-09-23 Thread Wensui Liu
MERGE DATA FRAMES BY 2 OR MORE VARIABLES ### # MERGE 2 DATA FRAMES BASED ON# # 2 OR MORE VARIABLES # ### data1-data.frame(x.id1 = 1:10, x.id2 = (1:10) * 2, x = rnorm(length(1:10)));

Re: [R] Data frames questions

2006-09-23 Thread Wensui Liu
3) How to I append to the bottom of a dataframe? APPEND ROWS ITERATIVELY TO A DATA FRAME # # APPEND ROWS ITERATIVELY TO A DATA FRAME # # # LOOP FROM 1 TO 10 for(i in 1:10) { # SET THE RANDOM

Re: [R] strange warning message

2006-09-23 Thread Duncan Murdoch
On 9/23/2006 7:15 PM, Tong Wang wrote: Hi Everyone, I recently start to get this warning message, while loading files in to R. Could someone tell me what does it mean ? I am using R 2.3.0 with Emacs on WinXP. use of NULL environment is deprecated The files were saved in an

Re: [R] Adding percentage to Pie Charts

2006-09-23 Thread Gabor Grothendieck
It might also be nice to be able to align the fans at the left or right, not just the center. On 9/23/06, Jim Lemon [EMAIL PROTECTED] wrote: Hi all, Anupam Tyagi mentioned an interesting idea a few days ago. A modification in a pie chart that draws overlapping areas with a common start

[R] s language version and complie languages

2006-09-23 Thread Taka Matzmoto
Hi R users Which version of S does the current R version use? Venables and Ripley's book (2000) says R is version 3 of S language. Is there any change since 2000 ? I searched R-help archive for using C or C++ with R. Although there are some postings, I am looking for up-to-date answers.

Re: [R] s language version and complie languages

2006-09-23 Thread xchen
Taka Matzmoto wrote: Hi R users Which version of S does the current R version use? Venables and Ripley's book (2000) says R is version 3 of S language. Is there any change since 2000 ? I searched R-help archive for using C or C++ with R. Although there are some postings, I am looking