Re: [R] help with hclust and cutree

2011-09-13 Thread Madeleine Seeland
Hello, yes, thanks! Also,how do I get the cutoff similarity value which was used for clustering when I cut the tree at a specific heigt? Best regards, Madeleine On Tue, 2011-09-13 at 09:39 -0700, Mikkel Grum wrote: > plot(clust) > rect.hclust(clust, h = 0.65) > > > > - Original Message

Re: [R] how to get a row with its its p value below 0.01 ??

2011-09-13 Thread Rolf Turner
On 14/09/11 17:49, anand m t wrote: Hi all, I'm analyzing micro array data.. it has produced a file (to be specific matrix) withdimension of 35556 2. first few lines of the matrix are as below.. probe_name control.fdr.pvals.present 10338001 0.000440001 10338002 0.000583093 10338003 0.000528449

[R] how to get a row with its its p value below 0.01 ??

2011-09-13 Thread anand m t
Hi all, I'm analyzing micro array data.. it has produced a file (to be specific matrix) withdimension of 35556 2. first few lines of the matrix are as below.. probe_name control.fdr.pvals.present 10338001 0.000440001 10338002 0.000583093 10338003 0.000528449 10338004 0.000610362 10338005 0.000151

[R] Difficulty with saving .Rhistory

2011-09-13 Thread T
Current systems     OS: RedHat SE Linux 5.7, fully patched     R: R-2.13.1 I am a new R user. The R installation needs to be very tight because I plan to use it after I am employed on Internet facing systems to analyse firewall and filesystem log data.  Everything to my knowledge regarding

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Martin Morgan
On 09/13/2011 10:54 AM, Joseph Park wrote: Hi, I'm looking for some guidance on whether to use S4 or Reference Classes for an analysis application I'm developing. I'm a C++/Python developer, and like to 'think' in OOD. I started my app with S4, thinking that was the best

Re: [R] Reading R Code aloud

2011-09-13 Thread cwdillon
Well, I've bee reading the R Manual, Nathan Yau's book on Flowing Data, and a few others like this to my kids to get them to go to sleep. I started off doing it the *wrong* way, just to keep it boring but slipped into doing it the *correct* way for my own sanity. It puts the kids to sleep in no tim

Re: [R] solving linear equations

2011-09-13 Thread Ben Bolker
varunshivashankar gmail.com> writes: > > I have a dataset > > X Y1 > 1200 1.375 > 4000 0.464 > 1333.33 0.148 > 444.440.047 > 148.148 0.014 > 49.3830.005 > 16.4610.004 > > I have to find a curve fit for the above dataset based on a 4-parameter > logist

Re: [R] drc results for binomial data

2011-09-13 Thread PatGauthier
I've been having trouble making sense of the drc results for my binomial response toxicity data. Firstly, the standard errors are far too large for how well the data fit the log-logistic model, particularly compared to other methods of LC50 estimation (e.g. probit and trimmed spearman-karber). Seco

Re: [R] count regression zero count comparison

2011-09-13 Thread jekang
Thank you, it helped me clear some confusions. Jamie -- View this message in context: http://r.789695.n4.nabble.com/count-regression-zero-count-comparison-tp3810907p3811277.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pro

[R] solving linear equations

2011-09-13 Thread varunshivashankar
I have a dataset X Y1 12001.375 40000.464 1333.33 0.148 444.44 0.047 148.148 0.014 49.383 0.005 16.461 0.004 I have to find a curve fit for the above dataset based on a 4-parameter logistic equation viz. Y1 = d + ((a-d)/(1+(X/cc)^b)), where X and Y1 are the values above. I nee

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Rich Shepard
On Tue, 13 Sep 2011, Gabor Grothendieck wrote: In the example data there were no duplicate dates for a given param value. Gabor, There are duplicate dates for a given param value, but each is on a different stream. It is the triplet (stream, sampdate, param) that is associated with each val

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Gabor Grothendieck
On Tue, Sep 13, 2011 at 7:36 PM, Rich Shepard wrote: > On Tue, 13 Sep 2011, Gabor Grothendieck wrote: > >> As in ?zoo a zoo object is a numeric matrix, numeric vector or factor >> together with an ordered time index which is unique. Its not clear that >> that is what you have; however, if we can a

Re: [R] stupid lm() question

2011-09-13 Thread Rolf Turner
``The only stupid question is the one that is not asked.'' --- Anon. (???) The reason that there is no value for level "A" of treatment is that level "A" is the reference level, under the default ``treatment'' contrasts. See ?contr.treatment. Some insight may be obtained by doing:

Re: [R] CMYK color space

2011-09-13 Thread ivo welch
obviously not. thank you, henrik. going back to square 101. (I just googled, because I usually find stuff faster by googling. cmyk and R brought up nothing.) apologies for the bandwidth, everyone. Ivo Welch (ivo.we...@gmail.com) On Tue, Sep 13, 2011 at 4:24 PM, Henrik Bengtsson wrote:

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Rich Shepard
On Tue, 13 Sep 2011, Gabor Grothendieck wrote: As in ?zoo a zoo object is a numeric matrix, numeric vector or factor together with an ordered time index which is unique. Its not clear that that is what you have; however, if we can assume that for each value of param we have a unique set of dates

Re: [R] CMYK color space

2011-09-13 Thread Henrik Bengtsson
Did you ask your question being aware of the following from help("pdf") in R v2.13.1: pdf(file = ifelse(onefile, "Rplots.pdf", "Rplot%03d.pdf"), width, height, onefile, family, title, fonts, version, paper, encoding, bg, fg, pointsize, pagecentre, colormodel, useDingbats, useKerning,

Re: [R] ZOO: Learning to apply it to my data

2011-09-13 Thread Gabor Grothendieck
On Tue, Sep 13, 2011 at 2:07 PM, Rich Shepard wrote: >  I have read ?zoo but am not sure how to relate the parameters (x, > order.by, frequency, and style) to my data.frame. The structure of the > data.frame is > > 'data.frame':   11169 obs. of  4 variables: >  $ stream  : Factor w/ 37 levels "Bur

Re: [R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 4:51 PM, Tianchan Niu wrote: Thank you David for your reply. But I am still confusing. I am using Wilcoxon test to compare two samples to assess whether their means differ because neither of them is normally distributed. Is it possible to compute the power of the Wilcoxon

Re: [R] function censReg in panel data setting

2011-09-13 Thread Arne Henningsen
Hi Igors On 13 September 2011 13:27, Igors wrote: > Any success in finding possible solutions for my problem? Somewhat. The calculation of the log-likelihood values is numerically much more robust/stable now. The log-likelihood contributions of some individuals became minus infinity in your mode

[R] CMYK color space

2011-09-13 Thread ivo welch
dear R experts---I am struggling with the requirements to prepare my files for my printers. I am printing in 2/2 format, which means cyan and black for me, which they take from my color-separated pdf files. R comes into play, because it produces all the figures that are embedded in my book (pdfla

[R] stupid lm() question

2011-09-13 Thread Carl Witthoft
I feel bad even asking, but: Rgames> data(OrchardSprays) Rgames> model<-lm(decrease~.,data=OrchardSprays) Rgames> model Call: lm(formula = decrease ~ ., data = OrchardSprays) Coefficients: (Intercept) rowpos colpos treatmentB treatmentC 22.705 -2.784 -1.234

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread William Dunlap
I often use the following function is.true <- function(x) !is.na(x) & x and, less often, is.false <- function(x) !is.na(x) & !x to report if elements of a logical vector are TRUE (not FALSE or NA) or FALSE (not TRUE or NA), respectively. Do your complicated logical expression and apply is.true

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Duncan Murdoch
On 11-09-13 5:17 PM, Timothy Bates wrote: Dear Duncan and Hadley, I stumbled across the NA behavior of subset a little while ago and thought it might do the trick. But my common usage case is not getting a subsetting sans NAs, but setting values in the whole dataframe. So I need T/F at each r

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Douglas Bates
On Tue, Sep 13, 2011 at 4:11 PM, Steve Lianoglou wrote: > Hi, > > On Tue, Sep 13, 2011 at 1:54 PM, Joseph Park wrote: >> >>   Hi, I'm looking for some guidance on whether to use >>   S4 or Reference Classes for an analysis application >>   I'm developing. >>   I'm a C++/Python developer, and like

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Douglas Bates
On Tue, Sep 13, 2011 at 12:54 PM, Joseph Park wrote: >   Hi, I'm looking for some guidance on whether to use >   S4 or Reference Classes for an analysis application >   I'm developing. >   I'm a C++/Python developer, and like to 'think' in OOD. >   I started my app with S4, thinking that was the

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Timothy Bates
Dear Duncan and Hadley, I stumbled across the NA behavior of subset a little while ago and thought it might do the trick. But my common usage case is not getting a subsetting sans NAs, but setting values in the whole dataframe. So I need T/F at each row, not just the list of rows that match th

Re: [R] S4 vs Reference Classes

2011-09-13 Thread Steve Lianoglou
Hi, On Tue, Sep 13, 2011 at 1:54 PM, Joseph Park wrote: > >   Hi, I'm looking for some guidance on whether to use >   S4 or Reference Classes for an analysis application >   I'm developing. >   I'm a C++/Python developer, and like to 'think' in OOD. >   I started my app with S4, thinking that was

Re: [R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread Tianchan Niu
Thank you David for your reply. But I am still confusing. I am using Wilcoxon test to compare two samples to assess whether their means differ because neither of them is normally distributed. Is it possible to compute the power of the Wilcoxon test similar to that of t test? Or is it just a wron

[R] Reading R Code aloud

2011-09-13 Thread andrewH
Dearfolks-- I have been told by an experienced R programmer and teacher whom I trust that it is easier to understand R code if you read it aloud, as the language that it is. However, she was clear that reading it aloud was not simply reading the marks on the screen: you read A.df[5,] as "the fift

Re: [R] count regression zero count comparison

2011-09-13 Thread Achim Zeileis
On Tue, 13 Sep 2011, jekang wrote: Hello, It might be more of a statistical question than an R question. I was reading http://cran.r-project.org/web/packages/pscl/vignettes/countreg.pdf, and I was wondering why the following functions were used to compare zero counts (observed and predicted),

Re: [R] space in directory name

2011-09-13 Thread Jean-Christophe BOUËTTÉ
And a final word of advice: be lazy, and use file.choose() to get the string that R expects as a filename. This will help avoid such typos! JC 2011/9/13 Mikkel Grum : > You are missing \\ between Documents and settings and Administrator. > > > - Original Message - > From: dbonneau > To:

Re: [R] import csv file into R, strange problem

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 12:58 PM, Greg Snow wrote: It would appear that your file is not a csv file, but we cannot tell for sure without seeing an example of the data that you are reading in, can you cut and paste the first few lines of the file? Or post a link to a copy of the file. If the d

[R] Online literature database search in R - RCurl?

2011-09-13 Thread Brian Oney
Hi All, theoretically, I should be able (with proper proxy setting and IP address) to connect and perform literature database searches, such as in the ISI Web of Science, with R. I can imagine this working nicely with the RCurl package. I am pestering the list to ask if I missed a certain pa

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 11:56 AM, David Winsemius wrote: On Sep 13, 2011, at 11:44 AM, David Winsemius wrote: On Sep 13, 2011, at 9:43 AM, RCulloch wrote: Dear John, Thank you for that, and for explaining why the abline() command wont/dosen't work. The approach is based on reviewers commen

Re: [R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 3:34 PM, tn85 wrote: Hello All, I perform a Wilcoxon signed rank test for two sets of data to test whether they two have significantly different means. I would also like to know the power of this test. Given that none of the various "Wilcoxon tests" are for differen

[R] count regression zero count comparison

2011-09-13 Thread jekang
Hello, It might be more of a statistical question than an R question. I was reading http://cran.r-project.org/web/packages/pscl/vignettes/countreg.pdf, and I was wondering why the following functions were used to compare zero counts (observed and predicted), instead of just using hist(fitted(fm_

[R] S4 vs Reference Classes

2011-09-13 Thread Joseph Park
Hi, I'm looking for some guidance on whether to use S4 or Reference Classes for an analysis application I'm developing. I'm a C++/Python developer, and like to 'think' in OOD. I started my app with S4, thinking that was the best set of OO features in R. However, it appears that o

Re: [R] space in directory name

2011-09-13 Thread Mikkel Grum
You are missing \\ between Documents and settings and Administrator. - Original Message - From: dbonneau To: r-help@r-project.org Cc: Sent: Tuesday, September 13, 2011 10:59 AM Subject: [R] space in directory name Hi, I am trying to read a text file located in following paths. I am get

Re: [R] writing a loop using several list()-objects

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 10:26 AM, travellina wrote: Hi there, I am looking for some help with a loop which uses three list()- objects: "sequence" has got numbers, "start" and "end" have got the begin/end time points which should be applied to the numbers in "sequence". length() of all three o

Re: [R] barplot in hexagram layout

2011-09-13 Thread Schatzi
I decided to go with circles instead of rectangles. Thank you for your help. Here is the new code: dev.new(width=2.5, height=3,mar=c(0,0,0,0)) par(mfrow=c(1,1),mar=c(0,0,0,0)) x=c(-1,1,1,-1,-3,-3,5) y=c(1.2,0.6,-.7,-1.3,-.7,0.6,5) plot(0,xlim=c(-4,2),ylim=c(-2,2),type="n",axes=FALSE,xlab="",ylab="

Re: [R] ESS: Viewing Restored Workspace At Startup

2011-09-13 Thread Rich Shepard
On Wed, 14 Sep 2011, Patrick Connolly wrote: The objects will be available already. See what you have with the ls() command. They are stored in .RData so the fact that you can't read .RData directly matters not. Patrick, I'm aware of .Rdata and use ls() regularly outside of R. What I did n

Re: [R] ESS: Viewing Restored Workspace At Startup

2011-09-13 Thread Patrick Connolly
On Tue, 13-Sep-2011 at 06:53AM -0700, Rich Shepard wrote: |> I use ESS and have no problems invoking it when I fire up emacs. When it |> loads it tells me, [Previously saved workspace restored]. However, I do not |> know how to view that history of commands so I can continue what I was |> doing.

Re: [R] envfit vector labels with ordiplot3d

2011-09-13 Thread Gavin Simpson
On Tue, 2011-09-13 at 05:18 -0700, Briony wrote: > It *really* does help to read the documentation for functions. Like most > of the functions in vegan that extract information from ordinations, > ordilabel() extracts information for "axes" 1 and 2 only, by default. It also helps if I read the doc

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Hadley Wickham
> Because in coding, I often end up with big chunks looking like this: > > ((mydataframeName$myvariableName > 2 & > !is.na(mydataframeName$myvariableName)) & > (mydataframeName$myotherVariableName == "male" & > !is.na(mydataframeName$myotherVariableName))) > > Which is much less readable/maintai

[R] How to calculate the power of Wilcoxon signed rank test

2011-09-13 Thread tn85
Hello All, I perform a Wilcoxon signed rank test for two sets of data to test whether they two have significantly different means. I would also like to know the power of this test. The third part of this tutorial is similar to what I want except the t distribution. http://www.cyclismo.org/tutori

Re: [R] space in directory name

2011-09-13 Thread Henrik Bengtsson
FYI, by using library("R.utils"); pathname <- "C:\\Documents and Settings Administrator\\My Documents\\My Dropbox\\data.txt"; pathname <- Arguments$getReadablePathname(pathname); data1 <- read.table(pathname, sep="\t", header=TRUE); you would get an explanation what is wrong and helps you troubl

Re: [R] space in directory name

2011-09-13 Thread Ben Bolker
Jean-Christophe BOUËTTÉ gmail.com> writes: > > I suggest changing > Documents and Settings Administrator > to > Documents and Settings\\Administrator > Please don't cross-post to R-help and Stack Overflow. People in both places will gladly answer well-posed questions, but it turns into a was

Re: [R] Min of

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 1:17 PM, bradford wrote: With the help of Andrie on StackOverflow.com, I was able to learn about ddply. I have another question that is more trivial and cannot seem to find help on IRC and do not want to bother Andrie again. It's doubtful that he would have considered

Re: [R] Convert "RDX" and "RDB" to ASCII format?

2011-09-13 Thread Duncan Murdoch
On 13/09/2011 2:28 PM, Henri Mone wrote: Dear R users and experts, I need to modify an exisiting R function from a package. The function/ variable is inside a namespace and I have problems accessing it. I only need to get access to to corresponding sourcecode. I localised the "RDX" and "RDB" fil

Re: [R] Selecting row indices from a data.frame by a factor and simple calculation (miniumum)

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 12:31 PM, Rawlins, Barry G. wrote: Hello I wish to extract the row indices from a data.frame in which a column contains numeric data by calculating the minimum value, but grouped on another column factor: An example data.frame: Codeabsdiff NY14/3070

Re: [R] Convert "RDX" and "RDB" to ASCII format?

2011-09-13 Thread Jean-Christophe BOUËTTÉ
Hi henri, R is open source, so you should be able to download the source code of the package from CRAN. JC 2011/9/13 Henri Mone : > Dear R users and experts, > > I need to modify an exisiting R function from a package. The function/ > variable is inside a namespace and I have problems accessing i

Re: [R] class weights with Random Forest

2011-09-13 Thread Liaw, Andy
The current "classwt" option in the randomForest package has been there since the beginning, and is different from how the official Fortran code (version 4 and later) implements class weights. It simply account for the class weights in the Gini index calculation when splitting nodes, exactly as

[R] Convert "RDX" and "RDB" to ASCII format?

2011-09-13 Thread Henri Mone
Dear R users and experts, I need to modify an exisiting R function from a package. The function/ variable is inside a namespace and I have problems accessing it. I only need to get access to to corresponding sourcecode. I localised the "RDX" and "RDB" files. But it seems to be a binary format. How

[R] ZOO: Learning to apply it to my data

2011-09-13 Thread Rich Shepard
I have read ?zoo but am not sure how to relate the parameters (x, order.by, frequency, and style) to my data.frame. The structure of the data.frame is 'data.frame': 11169 obs. of 4 variables: $ stream : Factor w/ 37 levels "Burns","CIL",..: 1 1 1 1 1 1 1 1 1 1 ... $ sampdate: Date, format

Re: [R] space in directory name

2011-09-13 Thread Jean-Christophe BOUËTTÉ
I suggest changing Documents and Settings Administrator to Documents and Settings\\Administrator ... JC 2011/9/13 dbonneau : > Hi, I am trying to read a text file located in following paths. I am getting > error if I try to read from long directory..  But if I place the file right > under C, It r

[R] Min of

2011-09-13 Thread bradford
With the help of Andrie on StackOverflow.com, I was able to learn about ddply. I have another question that is more trivial and cannot seem to find help on IRC and do not want to bother Andrie again. I can't seem to figure out what to google for, so I thought I'd ask here. I have: library(plyr)

Re: [R] Selecting row indices from a data.frame by a factor and simple calculation (miniumum)

2011-09-13 Thread Mikkel Grum
Code <- c(rep("NY14/3070", 3), rep("NY14/5459", 2)) Code <- as.factor(Code) absdiff <- c(2, 4, 1, 5, 7) df <- data.frame(Code, absdiff) which(     paste(df$Code, df$absdiff) ==      paste(         aggregate(df$absdiff, by = list(df$Code), min)$Group.1,         aggregate(df$absdiff, by = list(df$Cod

Re: [R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Duncan Murdoch
On 13/09/2011 12:42 PM, Timothy Bates wrote: Dear R cognoscenti, While having NA as a native type is nifty, it is annoying when making binary choices. Question: Is there anything bad about writing comparison functions that behavior like %in% (which I love) and ignore NAs? "%>%"<- function(ta

Re: [R] map

2011-09-13 Thread Justin Haynes
i responded offline the first time, but: google is your friend: search for R maps and you'll find what I mention below. In the future make sure to perform a thorough search of google and the help forums before you post That said... you're looking for the maps package install.packages('maps')

Re: [R] import csv file into R, strange problem

2011-09-13 Thread Greg Snow
It would appear that your file is not a csv file, but we cannot tell for sure without seeing an example of the data that you are reading in, can you cut and paste the first few lines of the file? Or post a link to a copy of the file. If the data is not able to be shared, then create a new (smal

Re: [R] Deleting Rows based on Factor and Time Period

2011-09-13 Thread Mikkel Grum
The following will get you the first stock in each week. Is that useful? install.packages("surveillance") library(surveillance) alldat$year <- isoWeekYear(alldat$mydate)$ISOYear alldat$week <- isoWeekYear(alldat$mydate)$ISOWeek alldat <- alldat[order(alldat$year, alldat$week), ] alldat[!duplicated

Re: [R] Loops on data˜1

2011-09-13 Thread Trying To learn again
Hi Michael, First of all thanks for your response. I do know that if I make my stimation on a single data the regression has no sense but it will be getting sense in the growing next estimations. I change my asking doubt. I want to use this regressions as a first filter. Only this. Can anyone

[R] x %>% y as an alternative to which( x > y)

2011-09-13 Thread Timothy Bates
Dear R cognoscenti, While having NA as a native type is nifty, it is annoying when making binary choices. Question: Is there anything bad about writing comparison functions that behavior like %in% (which I love) and ignore NAs? "%>%" <- function(table, x) { return(which(table > x)) }

[R] implicit data frame reference

2011-09-13 Thread Benjamin Polidore
If I create an aggregation like this: aggregate(lastYear[,8:10],list(Stadium=lastYear$STAD),mean) I'll get a new data frame, which I can order if I assign it like this: newFrame <- aggregate(lastYear[,8:10],list(Stadium=lastYear$STAD),mean) newFrame[order(newFrame$TEMP),] But.. if I just want t

[R] space in directory name

2011-09-13 Thread dbonneau
Hi, I am trying to read a text file located in following paths. I am getting error if I try to read from long directory.. But if I place the file right under C, It runs fine. Could anyone tell me how to read table from following directory path ? data1 <-read.table("C:\\Documents and Settings Adm

Re: [R] map

2011-09-13 Thread Batur
Adding to the previous question, I would like to map central Asia along with those five countries (Kazakhstan, Kyrgyzstan, Uzbekstan, Tajikstan and Turkmenstan). Please tell us the right data base!!! Thanks a lot!!! -- View this message in context: http://r.789695.n4.nabble.com/map-tp3810363p3810

[R] Building R package with precompiled shared library

2011-09-13 Thread Francesco Finazzi
Dear R users, we are trying to build a R package that includes a precompiled shared library, let's say mylib.so. We created the skeleton of the package and we moved the mylib.so file into the libs folder that we created at the same level of the folders man and R. Moreover we created the file NAMES

[R] Selecting row indices from a data.frame by a factor and simple calculation (miniumum)

2011-09-13 Thread Rawlins, Barry G.
Hello I wish to extract the row indices from a data.frame in which a column contains numeric data by calculating the minimum value, but grouped on another column factor: An example data.frame: Codeabsdiff NY14/3070 2 NY14/3070 4 NY14/3070 1 NY14/5459

[R] import csv file into R, strange problem

2011-09-13 Thread bonnieyuan
Hi, I used read.csv(file name, header=T, sep=",") to bring in a csv file I saved in MS Excel. The strange thing is all the data ended up in one big column. The number of rows match with the number of observations, but all the variables got squeezed into one column. Also the first row where the head

Re: [R] help with hclust and cutree

2011-09-13 Thread Mikkel Grum
plot(clust) rect.hclust(clust, h = 0.65) Is that what you wanted? - Original Message - From: Madeleine Seeland To: r-help@r-project.org Cc: Sent: Tuesday, September 13, 2011 9:25 AM Subject: [R] help with hclust and cutree Hello, I would like to cut a hclust tree into several groups

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 11:44 AM, David Winsemius wrote: On Sep 13, 2011, at 9:43 AM, RCulloch wrote: Dear John, Thank you for that, and for explaining why the abline() command wont/dosen't work. The approach is based on reviewers comments that I am a tad sceptical about myself but yet curi

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 9:43 AM, RCulloch wrote: Dear John, Thank you for that, and for explaining why the abline() command wont/ dosen't work. The approach is based on reviewers comments that I am a tad sceptical about myself but yet curious enough to test their suggestion..I don't think

[R] writing a loop using several list()-objects

2011-09-13 Thread travellina
Hi there, I am looking for some help with a loop which uses three list()-objects: "sequence" has got numbers, "start" and "end" have got the begin/end time points which should be applied to the numbers in "sequence". length() of all three objects is identical, length() of vectors stored in "sequen

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread RCulloch
David & JC, Excellent point, of course it does - and of course that is (should have been) obvious!!! That is what I get for taking a reviewers comment/suggestion as gospel without applying a bit of thought! I'm off to go and kick myself. Cheers, Ross -- View this message in context: htt

[R] ESS: Viewing Restored Workspace At Startup

2011-09-13 Thread Rich Shepard
I use ESS and have no problems invoking it when I fire up emacs. When it loads it tells me, [Previously saved workspace restored]. However, I do not know how to view that history of commands so I can continue what I was doing. Since the buffer begins with ESS, what do I do to view that restor

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread RCulloch
Dear John, Thank you for that, and for explaining why the abline() command wont/dosen't work. The approach is based on reviewers comments that I am a tad sceptical about myself but yet curious enough to test their suggestion..I don't think it is very straightforward to explain; however, it in

Re: [R] help with glmm.admb-it worked thanks!

2011-09-13 Thread eeadie
Hi Ben, Your suggestions about the offset worked. Thank you!! Lizzy -- View this message in context: http://r.789695.n4.nabble.com/help-with-glmm-admb-tp3788082p3810141.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pr

[R] Help in SVM prediction

2011-09-13 Thread Divyam
Hello, I am trying to use SVM from e1071 package for doing binary classification and I am having problems in prediction using SVM. I ran a SVM for X~as.factor(X1)+as.factor(X2)+as.factor(X3),data=data1,cross=10. str(data1) gives me data.frame':5040 obs. of 5 variables: $ X4: int 1 2 3

[R] map

2011-09-13 Thread n.via...@libero.it
Hi!I'm trying to plot the italian map with provincial administration but I don't find the right database and packages to do so.Does anyone help me and tell me how to do???Thanks for your attention [[alternative HTML version deleted]] __ R-help

Re: [R] using vif from package "car" - "aliased coefficients in the model"

2011-09-13 Thread Dimitri Liakhovitski
John, you were right. After I checked: I had MORE predictors that entered the model than predictors that were left in and that show up in in summary(myreg)$coeff). That means, some of the predictors were kicked out from the results. I have rerun my lm with all the predictors left in and then vif wo

[R] help with hclust

2011-09-13 Thread Madeleine Seeland
Hello, how can I get the similarity value (i.e., the inner cluster similarity) that was used to cut a hierarchical tree at a specific height? I would appreciate your help! Best regards, Madeleine __ R-help@r-project.org mailing list https://stat.ethz

Re: [R] Problems of metafile plots when converting word to pdf file

2011-09-13 Thread Jin Minming
Hi JC and Bill, Thanks a lot for your help. I found that Office 2007 and 2010 both have this problem when I choose to save as PDF file in Word program.  If I copy a graph as bitmap file, then everything is ok, but not look as nice as metafile.  I have tried to use other two tools (CutePDF and

Re: [R] using vif from package "car" - "aliased coefficients in the model"

2011-09-13 Thread John Fox
Dear Dimitri, > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Dimitri Liakhovitski > Sent: September-13-11 10:29 AM > To: r-help > Subject: [R] using vif from package "car" - "aliased coefficients in > the model" > > Hello! >

Re: [R] nls, the four parameter logistic equation, and prediction band

2011-09-13 Thread John C Nash
The error msg is telling you that R cannot evaluate the loss function, so you should not expect answers. You might try examining the data -- Are there NA or Inf entries? Or prepare a dataframe with just X and Y, sort by X and graph. Then check the nls computations by sampling, say, every 100 X'

[R] using vif from package "car" - "aliased coefficients in the model"

2011-09-13 Thread Dimitri Liakhovitski
Hello! I have run a simple regression - lm and created a regression object "myreg". I can see all the coefficients when I print(myreg). Then I tried to run vif(myreg) from the package "car". However, it's giving me an error: in vif.lm(regr.f) : there are aliased coefficients in the model Very sorr

[R] help with hclust and cutree

2011-09-13 Thread Madeleine Seeland
Hello, I would like to cut a hclust tree into several groups at a specific similarity. I assume this can be achieved by specifying the "h" argument with the specified similarity, e.g.: clust<-hclust(dist,"average") cut<-cutree(clust,h=0.65) Now, I would like to draw rectangles around the branche

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread John Fox
Dear Ross, But the residuals are just y - x. Best, John > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of RCulloch > Sent: September-13-11 9:19 AM > To: r-help@r-project.org > Subject: Re: [R] Force regression line to a 1:1

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread John Fox
Dear JC, > -Original Message- > From: Jean-Christophe BOUËTTÉ [mailto:jcboue...@gmail.com] > Sent: September-13-11 9:35 AM > To: John Fox > Cc: RCulloch; r-help@r-project.org > Subject: Re: [R] Force regression line to a 1:1 relationship > > And you can easily get these predictions using

[R] Deleting Rows based on Factor and Time Period

2011-09-13 Thread Anna Dunietz
Hi All! I have been messing around with this problem for about a week but to no avail! The following data has been cut down in order to make my question reproducible. The alldat data frame includes 2 columns: 1 date column and 1 factor column (equity names)). mydate<-as.Date(c("2001-07-02","2001

Re: [R] Problems of metafile plots when converting word to pdf file

2011-09-13 Thread Bill Pikounis
Jim, I have seen this outcome as well with graphs from time to time and agree with Jean-Christophe that the bug seems to creep in at the Windows / Office level. You do not mention your version of Word, but I have seen this with Word 2000 and Word 2007 (doc and docx) software and formats, for examp

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 9:34 AM, Jean-Christophe BOUËTTÉ wrote: And you can easily get these predictions using the following code : A <- B I had a similar thought in that: resA_B <- A - B would give the "residuals" from a line through the origin with a slope of 1. No regression needed

Re: [R] Mathematical expressions in the group labels in lattice

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 9:30 AM, David Winsemius wrote: On Sep 13, 2011, at 8:30 AM, Timo Schmid wrote: Hi, I am working with the lattice package and I want to label to groups in the xyplots with mathematical expressions. I short example for this library(lattice) Case<-factor(rep(1:2,each

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread Jean-Christophe BOUËTTÉ
And you can easily get these predictions using the following code : A <- B ;-) JC 2011/9/13 John Fox : > Dear Ross, > > lm(y ~ 0 + offset(x)) will do the trick, but the resulting model has no > coefficient estimates and thus can't be used with abline(). You can, e.g., > get predictions from the m

Re: [R] Matrix multiplication gives different results some of the times (NaN)

2011-09-13 Thread Xavier Fernández i Marín
Ben Bolker vas escriure el dia dt, 13 set 2011: > You mentioned that you had compiled R yourself. > Can you replicate this with a stock R binary, i.e. from a > Red Hat or Ubuntu repository? If not, can you give details of > your compilation setup, especially if you are using a specialized > or

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread RCulloch
yes, that is correct. The idea being that I want to know the residuals of the data points compared to a 1:1 line (as shown in the plot), if that makes sense? I appreciate that this might not be considered a typical approach, and it would probably take a while to explain (defend) why I am doing it!

Re: [R] Mathematical expressions in the group labels in lattice

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 8:30 AM, Timo Schmid wrote: Hi, I am working with the lattice package and I want to label to groups in the xyplots with mathematical expressions. I short example for this library(lattice) Case<-factor(rep(1:2,each=300)) xx<-rnorm(600,0,1) yy<-rnorm(600,0,1) xyplot(yy~

Re: [R] Problems of metafile plots when converting word to pdf file

2011-09-13 Thread Jean-Christophe BOUËTTÉ
Hi there, did you try saving your graph directly as pdf with R, using ?pdf What you're describing seems more a Word problem. JC 2011/9/13 Jin Minming : > Hi all, > > I copy metafile boxplot from R to Word. Then save as Pdf file. But I found > there are some unexpected black lines in some plots wit

Re: [R] Force regression line to a 1:1 relationship

2011-09-13 Thread John Fox
Dear Ross, lm(y ~ 0 + offset(x)) will do the trick, but the resulting model has no coefficient estimates and thus can't be used with abline(). You can, e.g., get predictions from the model, but I'm not sure what real use it will be to you. I hope this helps, John ---

Re: [R] How to use a variable after $ ($x) for subscripting a list, How to source from a character object (source(x))

2011-09-13 Thread David Winsemius
On Sep 13, 2011, at 4:11 AM, drflxms wrote: Dear R colleagues, as result of a function a huge list is generated. From this result list I'd like to extract information. Think of the list i.e. as an object named "listResult" with the following form: [[a]] [1] [2] [3] [4] [5] [[b]] [1] [2]

Re: [R] GO & Protein Complex Analysis for Homo sapiens

2011-09-13 Thread Martin Morgan
On 09/13/2011 05:21 AM, Sandeep Amberkar wrote: Dear All, I need to fetch GO ontologies for Homo sapiens with their mappings to corresponding Uniprot identifiers. I would be using this information to compare result from a clustering algorithm with existing protein complexes. This would be a test

  1   2   >