Re: [R] Is there a way to print branch distances for hclust function?

On Sun, Dec 11, 2011 at 8:43 PM, kbrownk wrote: > The R function hclust is used to do cluster analysis, but based on R > help I see no way to print the actual fusion distances (that is, the > vertical distances for each connected branch pairs seen in the cluster > dendrogram). > > Any ideas? I'd l

[R] Colours for sunflowerplot

Dear fellow R users, I would like to draw a "sunflowerplot" because I have data (decade by month) that plots multiple times on the same x-y co-ordinates. Further I would like to colour each of the points/sunflower leaves on the plot according to the group they belong to (i.e. which type of even

Re: [R] question about spaces in r

Hi > > First of all, it's R, not r, and on this mailing list people care about > this kind of thing. > > Second, you will need to provide more information in order to get better > help. Please read the posting guide. > > There are a number of introductory level documents available via CRAN, > pl

[R] Is there a way to print branch distances for hclust function?

The R function hclust is used to do cluster analysis, but based on R help I see no way to print the actual fusion distances (that is, the vertical distances for each connected branch pairs seen in the cluster dendrogram). Any ideas? I'd like to use them test for significant differences from the me

[R] Color2D.matplot uniform color range

Hi all. I am having a problem using color2D.matplot I am trying to visualize several different matrices under two color ranges. One color range corresponds to values less than 1 and the second color range for values greater than 1. However, the minimum value of each matrix differs and is automatica

Re: [R] color2D.maplot fixed color range

Hello, Were you ever able to find a solution to your problem? I have been trying on and off to solve this problem for several months now. I would love to hear if you found a solution. -- View this message in context: http://r.789695.n4.nabble.com/color2D-maplot-fixed-color-range-tp3298679p41

Re: [R] how to assign a value?

On 2011-12-12 10:58, David Winsemius wrote: On Dec 11, 2011, at 9:07 PM, Jinsong Zhao wrote: On 2011-12-12 0:00, David Winsemius wrote: On Dec 11, 2011, at 10:27 AM, Jinsong Zhao wrote: Hi there, I hope to modify values in a vector or matrix in the following code: for (i in 1:9) { assign

Re: [R] Labelling R plots using Greek letters

A short answer is you can consider the tikzDevice package; then a long (really really long) answer is this: http://yihui.github.com/knitr/demo/graphics/ Regards, Yihui -- Yihui Xie Phone: 515-294-2465 Web: http://yihui.name Department of Statistics, Iowa State University 2215 Snedecor Hall, Ames,

Re: [R] how to assign a value?

On Dec 11, 2011, at 9:07 PM, Jinsong Zhao wrote: On 2011-12-12 0:00, David Winsemius wrote: On Dec 11, 2011, at 10:27 AM, Jinsong Zhao wrote: Hi there, I hope to modify values in a vector or matrix in the following code: for (i in 1:9) { assign(paste("a.", i, sep = ""), 1:i) get(paste("a.

Re: [R] Question about fitting seasonal ARIMA in R?

sorry for the re-post... On Sun, Dec 11, 2011 at 8:22 PM, Michael wrote: > Hi all, > > I just couldn't find a R function which can fit multiple seasonal > patters... i.e. in the following code: > > *arima(x = data, order = c(p, d, q), seasonal = list(order = c(P, D, Q), > period = S), ... > ***

Re: [R] how to assign a value?

On 2011-12-12 1:16, Patrick Burns wrote: You are basically in R Inferno Circle 8.1.40. http://www.burns-stat.com/pages/Tutor/R_inferno.pdf Thanks. R_inferno is a good material for me. In this document, there is several sections titled "string not the name". I try try to change get(paste("

[R] Question about fitting seasonal ARIMA in R?

Hi all, I just couldn't find a R function which can fit multiple seasonal patters... i.e. in the following code: *arima(x = data, order = c(p, d, q), seasonal = list(order = c(P, D, Q), period = S), ... *** * there can be only one "period", am I right? What if the data seem to have three differ

[R] Lagged values problem requiring short solution time

Hello, I am hoping someone can help tackle the problem below, for which I require a fast solution. It feels like there should be an elegant approach, but I am drawing blanks. Take a vector 'x' with random values > 0: x = runif(10,1,5) Assume some reasonably small positive value 'delta': delt

Re: [R] how to assign a value?

On 2011-12-12 0:00, David Winsemius wrote: On Dec 11, 2011, at 10:27 AM, Jinsong Zhao wrote: Hi there, I hope to modify values in a vector or matrix in the following code: for (i in 1:9) { assign(paste("a.", i, sep = ""), 1:i) get(paste("a.", i, sep = ""))[i] <- i+50 } Just one matrix? The

[R] limit ranges in hexbin

Hello everybody, I hope you can give me some help with limiting the ranges in x, y, and z for a hexbin plot. All I have found on the net is an unanswered message to this list from last year, so I hope my problem is not too stupid. I would like to plot some data using hexagonal binning. Currently

Re: [R] Labelling R plots using Greek letters

On Dec 11, 2011, at 6:38 PM, john james wrote: Dear R Users, Please I have the following query. I want to label one of the axes of my graph with the follwing latex expression- \beta^{\prime}x, i.e I have the transpose of beta. How do I go abt this. My second query is similar but it has

Re: [R] Kaplan-Meier survfit problem

On Dec 11, 2011, at 5:48 PM, Esteban Cervetto wrote: 2011/12/11 Esteban Cervetto I am working with uncensored data. I have duration of workers compensation. Then I have for each the number of days thet it doesn't work. This sample is not censored at right because I query only work acci

[R] Labelling R plots using Greek letters

Dear R Users, Please I have the following query. I want to label  one of the axes of my graph with the follwing latex expression- \beta^{\prime}x, i.e I have the transpose of beta. How do I go abt this. My second query is similar but it has to do with the conditioning symbol. The other axis

[R] using dcast to reshape a DF from long to wide with multiple measured variables per obs

I have data in the following format:   person<- c(1,1,1,1,2,2,2,2,2,3,3,3,3,3,3) v2<- c("2011-01-01", "2011-02-01", "2011-03-01", "2011-04-01", "2011-01-01", "2011-02-01", "2011-03-01", "2011-04-01", "2011-05-01", "2011-01-01", "2011-02-01", "2011-03-01", "2011-04-01", "2011-05-01", "2011-06-01")

[R] Kaplan-Meier survfit problem

2011/12/11 Esteban Cervetto > I am working with uncensored data. > > I have duration of workers compensation. Then I have for each the number > of days thet it doesn't work. > This sample is not censored at right because I query only work accidents > with date of return of work (saned) > > That i

[R] Boston/Cambridge -- Statistical Programming Language Technology Breakthroughs

If you are a statistician or researcher working in Boston/Cambridge, and you have a strong interest in breakthroughs in statistical programming language technology, contact me. Robert __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/lis

[R] nls start values

I'm using nls to fit periodic gene-expression data to sine waves. I need to set the upper and lower boundaries, because I do not want any negative phase and amplitude solutions. This means that I have to use the "port" algorithm. The problem is, that depending on what start value I choose for phase

Re: [R] how to assign a value?

I find that get() and assign() are awkward to use and that the syntax is easier if you put your objects into a list or environment. To me, it also makes it clearer what the code is doing and keeps the output of objects() shorter and easier to manage. E.g., nResults <- 9 results <- vector("list",

Re: [R] some CRAN mirrors not accessible

Le 11/12/2011 18:57, Prof Brian Ripley a écrit : First, see the status links in the first para of http://cran.r-project.org/mirrors.html . So it seems that the issue is not that the mirror is not accessible, but the part you are looking for is not current/available. yes indeed. For one "R fo

Re: [R] some CRAN mirrors not accessible

First, see the status links in the first para of http://cran.r-project.org/mirrors.html . So it seems that the issue is not that the mirror is not accessible, but the part you are looking for is not current/available. Second, the mirror list is http://cran.r-project.org/CRAN_mirrors.csv and t

Re: [R] multiple comparison of interaction of ANCOVA

Thank you for you use of HH. I think the right graph for this data is the much simpler ancova function library(HH) ancova(y ~ year * Trt, data=mydata) where we see that the three treatments have totally different slopes. The WoodEnergy example doesn't apply here. The WoodEnergy example illus

[R] some CRAN mirrors not accessible

Dear all, Since some weeks, look like the following CRAN mirrors are no longer accessible for package update: http://cran.univ-lyon1.fr http://mirror.ibcp.fr/pub/CRAN/ > update.packages(ask='graphics',checkBuilt=TRUE) Warning: unable to access index for repository http://cran.univ-lyon1.fr/b

Re: [R] efficiently finding the integrals of a sequence of functions

Thanks,Hans! I agree that this is a good way of solving this problem. Here is another way. Instead of defining a vector of uni-dimensional functions and trying to integrating each component (a uni-dimensional function), we can do something below my.integrand<-function(x,k) { return(f[x,k]) ## u

Re: [R] how to assign a value?

You are basically in R Inferno Circle 8.1.40. http://www.burns-stat.com/pages/Tutor/R_inferno.pdf On 11/12/2011 15:27, Jinsong Zhao wrote: Hi there, I hope to modify values in a vector or matrix in the following code: for (i in 1:9) { assign(paste("a.", i, sep = ""), 1:i) get(paste("a.", i, s

Re: [R] how to assign a value?

On Dec 11, 2011, at 10:27 AM, Jinsong Zhao wrote: Hi there, I hope to modify values in a vector or matrix in the following code: for (i in 1:9) { assign(paste("a.", i, sep = ""), 1:i) get(paste("a.", i, sep = ""))[i] <- i+50 } Just one matrix? Then you seem to have inappropriately borro

[R] how to assign a value?

Hi there, I hope to modify values in a vector or matrix in the following code: for (i in 1:9) { assign(paste("a.", i, sep = ""), 1:i) get(paste("a.", i, sep = ""))[i] <- i+50 } I get the following error message: Error in get(paste("a.", i, sep = ""))[i] <- i + 50 : target of assignment

Re: [R] multiple comparison of interaction of ANCOVA

Inline below. -- Bert On Sun, Dec 11, 2011 at 4:15 AM, Jinsong Zhao wrote: > Hi there, > > The following data is obtained from a long-term experiments. > >> mydata <- read.table(textConnection(" > +        y year Trt > +     9.37 1993   A > +     8.21 1995   A > +     8.11 1999   A > +     7.22 2

Re: [R] Overlaying density plot on forest plot

At 07:16 10/12/2011, Frank Peter wrote: Dear R User, Please, I am new to R. I want to overlay density plot for predictive interval pooled result in meta-analysis. http://addictedtor.free.fr/graphiques/graphcode.php?graph=114 It is hard to be sure from your rather brief question but does the

Re: [R] no non-missing arguments to max; returning -Inf

Yes - E-mail us reproducible code... Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English) -

Re: [R] Titelpage pdf-manual, packaging

On 11/12/2011 12:34, Johannes Radinger wrote: Am 10.12.2011 um 16:50 schrieb Uwe Ligges: On 09.12.2011 14:15, Johannes Radinger wrote: Hi, I am trying to get write my first own package. I followed the instructions in "Creating R Packages: A Tutorial" and "Writing R Extensions". So far ever

[R] no non-missing arguments to max; returning -Inf

Hi, I just wanted to object of class xts, zoo and I got the error no non-missing arguments to max; returning -Inf After all the values were printed. Tried searching on this forum, but no luck. Anyone has any idea ? Thanks. George -- View this message in context: http://r.789695.n4.nabble.com/no

[R] Bioconductor. MA plot for qPCR array

Dear all, Is there anyway too generate MA plot for 2 qPCR assays (an array of 2x 400). -- View this message in context: http://r.789695.n4.nabble.com/Bioconductor-MA-plot-for-qPCR-array-tp4182805p4182805.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Titelpage pdf-manual, packaging

Am 10.12.2011 um 16:50 schrieb Uwe Ligges: > > > On 09.12.2011 14:15, Johannes Radinger wrote: >> Hi, >> >> I am trying to get write my first own package. >> I followed the instructions in "Creating R Packages: A Tutorial" >> and "Writing R Extensions". So far everything works really >> fine,

[R] multiple comparison of interaction of ANCOVA

Hi there, The following data is obtained from a long-term experiments. > mydata <- read.table(textConnection(" +y year Trt + 9.37 1993 A + 8.21 1995 A + 8.11 1999 A + 7.22 2007 A + 7.81 2010 A +10.85 1993 B +12.83 1995 B +13.21 1999 B +

Re: [R] incomplete final line found warning

Xiaobo.Gu wrote > > Hi, > > I saved the following as a UTF-8 encoded file named amberutil.r > > as.factor.loop <- function(df, cols){ > > if (!is.null(df) && !is.null(cols) && length(cols) > 0) > { > for(col in cols) > { > df[[col]] <- as.fac

[R] Any packages or ways available providing un-/a-/non-symmetric centered distributions?

Dear R folks, I would like to use un-/a-/non-symmetric centered, i. e. the expectation vanishes/is zero, distributions. Searching for that is not that easy since all three words un-, a-, non-symmetric seem to be common. Do I need to create such a distribution myself by for example composing a d

Re: [R] incomplete final line found warning

On Sun, Dec 11, 2011 at 5:03 PM, Prof Brian Ripley wrote: > On Sun, 11 Dec 2011, David Winsemius wrote: > >> >> On Dec 10, 2011, at 10:01 PM, Xiaobo Gu wrote: > > > without following the posting guide in several respects and hence leaving us > guessing > > >>> Hi, >>> >>> I saved the followin

Re: [R] incomplete final line found warning

On Sun, 11 Dec 2011, David Winsemius wrote: On Dec 10, 2011, at 10:01 PM, Xiaobo Gu wrote: without following the posting guide in several respects and hence leaving us guessing Hi, I saved the following as a UTF-8 encoded file named amberutil.r BTW, it is hard to know how you know