Re: [R] [FORGED] Re: Dunnett Test in 'multicomp' package

2015-06-04 Thread Rolf Turner
On 05/06/15 11:08, Jim Lemon wrote: Hi James, You can change the order of levels like this: levels(viagraData$dose)<-c("placebo","low dose","high dose") As Richard Heiberger has pointed out, this is wrong. What *does* work is: viagraData$dose)<-factor(viagraData$dose,

Re: [R] Dunnett Test in 'multicomp' package

2015-06-04 Thread Richard M. Heiberger
Sorry, Jim that doesn't do the right thing. > tmp <- factor(c("mm", "cm", "dm", "m", "km")) > levels(tmp) [1] "cm" "dm" "km" "m" "mm" > tmp [1] mm cm dm m km Levels: cm dm km m mm > levels(tmp) <- c("mm", "cm", "dm", "m", "km") > tmp [1] km mm cm m dm Levels: mm cm dm m km > ## back to the beg

Re: [R] Dunnett Test in 'multicomp' package

2015-06-04 Thread Jim Lemon
Hi James, You can change the order of levels like this: levels(viagraData$dose)<-c("placebo","low dose","high dose") Although I don't know the exact names of the variable and its levels. Jim On Fri, Jun 5, 2015 at 7:19 AM, Richard M. Heiberger wrote: > This example is based on ?glht > >> data(

Re: [R] Finite mixture model: truncated normal and a normal distribution

2015-06-04 Thread David Winsemius
On Jun 4, 2015, at 5:38 PM, John Sorkin wrote: > I am looking for software that will take data what appears to come from a > series of low numbers coming from a left-truncated normal distribution and a > series of higher numbers coming from a normal distribution and produce the > mean, SD, and

[R] Finite mixture model: truncated normal and a normal distribution

2015-06-04 Thread John Sorkin
I am looking for software that will take data what appears to come from a series of low numbers coming from a left-truncated normal distribution and a series of higher numbers coming from a normal distribution and produce the mean, SD, and mixing proportion of the two source distributions. I hav

Re: [R] web scraping image

2015-06-04 Thread Jim Lemon
Hi Chris, I don't have the packages you are using, but tracing this indicates that the page source contains the relative path of the graphic, in this case: /nwisweb/data/img/USGS.12144500.19581112.20140309..0.peak.pres.gif and you already have the server URL: nwis.waterdata.usgs.gov getting the

Re: [R] Help to solve an Error

2015-06-04 Thread Jim Lemon
Hi all, Does this have something to do with "markdown"? I recently had to write something using this and found it much more difficult than HTML. One missed space or character, which is often hard to see, can turn a whole line to garbage. Jim On Thu, Jun 4, 2015 at 9:38 PM, John Kane wrote: > Hi

Re: [R] Decision tree in R using csv files

2015-06-04 Thread John Kane
No attachment. See http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example and http://adv-r.had.co.nz/Reproducibility.html John Kane Kingston ON Canada > -Original Message- > From: aks...@cosmicad.com > Sent: Thu, 4 Jun 2015 08:55:42 -0700 (PDT) > To: r-he

Re: [R] Dunnett Test in 'multicomp' package

2015-06-04 Thread Richard M. Heiberger
This example is based on ?glht > data(warpbreaks) > glht(amod, linfct = mcp(tension = "Dunnett")) General Linear Hypotheses Multiple Comparisons of Means: Dunnett Contrasts Linear Hypotheses: Estimate M - L == 0-10.0 H - L == 0-14.7 > levels(warpbreaks$tension) [1] "L" "M" "

Re: [R] Unwanted unicode - solved

2015-06-04 Thread Therneau, Terry M., Ph.D.
Thanks to all. The suggestion below by Marc Schwarz (the second I tried) showed the problem. One of the references in one of the files had been pasted in and had a funky dash in its "111-196" page number. It looked just fine in my emacs window so I hadn't picked it up. There are 90 .Rd files

[R] Dunnett Test in 'multicomp' package

2015-06-04 Thread James F. Henson
Greetings Below is my code. library("multcomp") viaModel1 <- aov(libido ~ dose, data=viagraData) dunnettModel <- glht(viaModel1 , linfct = mcp(dose = "Dunnett"), base = "placebo") The code base="placebo" is ignored. All treatments are compared to the first treatment in the order, which is "hi

Re: [R] Reproducing a 3d yield curve plot from New York Times

2015-06-04 Thread David Winsemius
At the bottom you asked for package recs: Karline Soetaert's (excuse any misspelling) plot3D package should be examined. I didn't find the fine-grained example that I remember seeing but this is an example in RGL code: http://markmail.org/message/wgetbk5h3f7zrfd4 -- David. On Jun 4, 2015,

[R] Reproducing a 3d yield curve plot from New York Times

2015-06-04 Thread matt
I am trying to closely reproduce the 3d yield curve chart that showed up in the New York Times in March, as seen here http://www.nytimes.com/interactive/2015/03/19/upshot/3d-yield-curve-economic-growth.html?ref=economy&abt=0002&abg=0 Working with chartSeries3d0 gets me close, but not quite th

Re: [R] Unwanted unicode

2015-06-04 Thread Marc Schwartz
On Jun 4, 2015, at 12:56 PM, Therneau, Terry M., Ph.D. wrote: > > I'm checking the survival package and get the following error. How do I find > the offending line? (There are a LOT of files in the man directory.) > > Terry T. > > -- > > * checking PDF version of manual ...

Re: [R] Unwanted unicode

2015-06-04 Thread Ista Zahn
If you use emacs you can use 'M-x find-grep-dired', select the directory, and search for '‑' (maybe won't work on Windows, I'm not sure). If you are on Linux (OS X?) you can run the equivalent find . \( -type f -exec grep -q -e \‑ \{\} \; \) -ls in a terminal. Best, Ista On Thu, Jun 4, 2015 at

[R] Unwanted unicode

2015-06-04 Thread Therneau, Terry M., Ph.D.
I'm checking the survival package and get the following error. How do I find the offending line? (There are a LOT of files in the man directory.) Terry T. -- * checking PDF version of manual ... WARNING LaTeX errors when creating PDF version. This typically indicates Rd proble

Re: [R] is.na for S4 object

2015-06-04 Thread Martin Morgan
On 06/04/2015 10:08 AM, cgenolin wrote: Hi the list, I have a variable y that is either NA or some S4 object. I would like to know in which case I am, but it seems taht is.na does not work with S4 object, I get a warnings: --- 8< setClass("myClass",slots=c(x="numeric")) if(runif(1)

[R] is.na for S4 object

2015-06-04 Thread cgenolin
Hi the list, I have a variable y that is either NA or some S4 object. I would like to know in which case I am, but it seems taht is.na does not work with S4 object, I get a warnings: --- 8< setClass("myClass",slots=c(x="numeric")) if(runif(1)>0.5){a <- new("myClass")}else{a <- NA} is

Re: [R] Error running help.search("keywords")

2015-06-04 Thread Duncan Murdoch
On 04/06/2015 6:42 AM, Yo Gmail wrote: Hi, When running help.search(“linear algebra”) or any other valid search keyword the following error is being returned: Error in help(db[i, "topic"], package = db[i, "Package"], lib.loc = lib, : 'topic' should be a name, length-one character vector or

[R] web scraping image

2015-06-04 Thread Curtis DeGasperi
I'm working on a script that downloads data from the USGS NWIS server. dataRetrieval makes it easy to quickly get the data in a neat tabular format, but I was also interested in getting the tabular text files - also fairly easy for me using download.file. However, I'm not skilled enough to work ou

[R] Decision tree in R using csv files

2015-06-04 Thread AKSHAYLAMBA
Hi Pls check the attachment. I am having some error while making a decision tree in R . Pls help. -- View this message in context: http://r.789695.n4.nabble.com/Decision-tree-in-R-using-csv-files-tp4708193.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] cannot run package anymore

2015-06-04 Thread Rich Shepard
On Thu, 4 Jun 2015, T.Riedle wrote: Error: package 'sandwich' could not be loaded I have no idea why this happens. I have not done anything before. Er, the answer is shown. What can I do to solve this problem? Install the latest sandwich (hold the mayo), then re-install the package req

[R] Error running help.search("keywords")

2015-06-04 Thread Yo Gmail
Hi, When running help.search(“linear algebra”) or any other valid search keyword the following error is being returned: Error in help(db[i, "topic"], package = db[i, "Package"], lib.loc = lib, : 'topic' should be a name, length-one character vector or reserved word I have not installed the

[R] cannot run package anymore

2015-06-04 Thread T.Riedle
Hello, I have been using the midasr package for several month. Now I want to start it again and it does not work anymore. I have installed the newest R version but I get following message when I am trying to run midasr install.packages("midasr") trying URL 'http://cran.rstudio.com/bin/windows/co

Re: [R] Splitting and arranging in ascending order

2015-06-04 Thread Amelia Marsh
Dear Sir, Thanks for your guidance and sinecerly apologize for late reply from my end. Regards Amelia On Wednesday, 3 June 2015 10:38 PM, Rui Barradas wrote: Hello, Try the following. tmp <- strsplit(as.character(my_dat$instrument), "_") tmp <- t(as.data.frame(tmp)) tmp <- data.frame(inst

Re: [R] reshape a data frame

2015-06-04 Thread David L Carlson
Here are two ways using stack() and reshape() from base R. The "col" variable indicates which column the bcX value came from. Both are easy to scale up to multiple columns: > exnew1 <- data.frame(ex$gIN, stack(ex, 2:3), ex$group) > names(exnew1) <- c("gIn", "bcX", "col", "group") > exnew1 gIn

Re: [R] R convert pdf/png to html

2015-06-04 Thread Boris Steipe
I'm glad this works for you. Myself, I would use paste() to collect the output to a single variable, then use one instance of cat() to write the whole thing to file. This saves you retyping ... file="treemap.html", append=TRUE,sep='\n' ... more than once and prevents errors such as the one you e

Re: [R] Passing strings with spaces to Python using system2 splits string at whitespace

2015-06-04 Thread John McKown
On Thu, Jun 4, 2015 at 8:47 AM, Wall, Wade A ERDC-RDE-CERL-IL < wade.a.w...@usace.army.mil> wrote: > Hi all, > > I am trying to pass arguments to a python script using R, but am running > into a problem with the string being split on the white spaces. > Investigation on the python end suggests tha

Re: [R] Passing strings with spaces to Python using system2 splits string at whitespace

2015-06-04 Thread Marc Schwartz
On Jun 4, 2015, at 8:47 AM, Wall, Wade A ERDC-RDE-CERL-IL wrote: > > Hi all, > > I am trying to pass arguments to a python script using R, but am running into > a problem with the string being split on the white spaces. Investigation on > the python end suggests that it is happening upstream

Re: [R] reshape a data frame

2015-06-04 Thread Bert Gunter
Yes. This is basic stuff, and it seems unnecessary to run to packages for it, Knowledge of base R should suffice. It would appear that the OP would benefit by going through an R tutorial or two. Slightly more economical and more general -- and trickier -- than explicit concatenation, which could g

[R] Passing strings with spaces to Python using system2 splits string at whitespace

2015-06-04 Thread Wall, Wade A ERDC-RDE-CERL-IL
Hi all, I am trying to pass arguments to a python script using R, but am running into a problem with the string being split on the white spaces. Investigation on the python end suggests that it is happening upstream from python, because other shells such as bash have generated similar errors.

Re: [R] SVM error

2015-06-04 Thread Uwe Ligges
On 04.06.2015 08:06, Pijush Das wrote: Dear Sir, I have converted the data into numeric as the function read.excel read the data file as a character (as it was said by David Meyer, maintainer of the package e1071) and after that run the function svm. Then another error is found there. The deta

Re: [R] Help to solve an Error

2015-06-04 Thread John Kane
Hi Jim, It is an problem in Gavin Simpson's code itself. Bad translation from R to HMTL it appears John Kane Kingston ON Canada > -Original Message- > From: drjimle...@gmail.com > Sent: Thu, 4 Jun 2015 21:22:37 +1000 > To: jacksonmrodrig...@gmail.com > Subject: Re: [R] Help to solve an

Re: [R] Help to solve an Error

2015-06-04 Thread John Kane
It seems to be some kind of "translation" error from R (or the text editior) to HMTL. That > does not work. Try grecent <- subset(gtemp, Year >= 1995, select = c(Year, Annual)) John Kane Kingston ON Canada > -Original Message- > From: jacksonmrodrig...@gmail.com

Re: [R] Help to solve an Error

2015-06-04 Thread Jim Lemon
Hi Jackson, It looks like you have picked up the HTML code for the right angle bracket. Try replacing this with a right angle bracket: grecent<-subset(gtemp,subset = Year>=1995,select = c(Year,Annual)) Jim On Thu, Jun 4, 2015 at 8:20 PM, Jackson Rodrigues wrote: > Hi, > > > > I want to apply t

Re: [R] Help to solve an Error

2015-06-04 Thread Sarah Goslee
The dangers of HTML: > should be > a greater than symbol. grecent <- subset(gtemp, subset = Year >= 1995,select = c(Year, Annual)) I've copied Gavin so he can decide whether to fix it. Sarah On Thursday, June 4, 2015, Jackson Rodrigues wrote: > Hi, > > > > I want to apply the codes of Gavin S

[R] Help to solve an Error

2015-06-04 Thread Jackson Rodrigues
Hi, I want to apply the codes of Gavin Simpson from http://www.fromthebottomoftheheap.net/2011/06/11/global-warming-since-1995-now-significant/ on my own data to detect trends in a subset. However, when I run: > grecent <- subset(gtemp, subset = Year >= 1995,select = c(Year, Annual)) I get err

Re: [R] Difference between 32-bit and 64-bit version

2015-06-04 Thread Thierry Onkelinx
"low probability of occurring" was just statisticians lingo for "rare" ;-) ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call

Re: [R] Save the result of map.market function to HTML file?

2015-06-04 Thread Jim Lemon
Hi valerio, This is a great coincidence. I have been working on exactly this problem in an analysis I have been conducting. When I supply a character variable to the png function rather than an explicit string, I get the same error in htmlize. Analyzing the error message, I realized that it was not

Re: [R] Difference between 32-bit and 64-bit version

2015-06-04 Thread Duncan Murdoch
On 04/06/2015 3:59 AM, Thierry Onkelinx wrote: > Dear Duncan, > > I had been thinking about FAQ 7.31. I tried to create a dummy dataset > with the same structure to replicate the problem with the need of > sending my dataset. However all of them gave identical() results between > 32-bit and 64-bit

[R] Vegan_p-values_RDA

2015-06-04 Thread Sibylle Stöckli
Dear R members, I understand the main principles why R-Vegan does not provide p-values for the biplot scores and/or canonical coefficients (see also post on stackoverflow). (i) We can obtain linear regression statistics and refit an ordination result as multiple response linear model (lm, see a

Re: [R] R convert pdf/png to html

2015-06-04 Thread valerio orfano
I forgot to use fie= treemap.html :) On 04 Jun 2015, at 10:53, valerio orfano wrote: > I have the following code: > > setwd("C:/Users/Administrator/Desktop") > cat("",file="treemap.html", append=TRUE,sep='\n') > cat("",file="treemap.html", append=TRUE,sep='\n') > cat("",file="treemap.html", ap

Re: [R] R convert pdf/png to html

2015-06-04 Thread valerio orfano
I have the following code: setwd("C:/Users/Administrator/Desktop") cat("",file="treemap.html", append=TRUE,sep='\n') cat("",file="treemap.html", append=TRUE,sep='\n') cat("",file="treemap.html", append=TRUE,sep='\n') cat("TreeMap",file="treemap.html", append=TRUE,sep='\n') cat("","treemap.html", a

Re: [R] Greek letters with subscript numbers on x-axis

2015-06-04 Thread Frank S.
Thank you very much Bill, Your answer to my question is exactly what I was trying to do in my R code. Best regards. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNS

Re: [R] Difference between 32-bit and 64-bit version

2015-06-04 Thread Thierry Onkelinx
Dear Duncan, I had been thinking about FAQ 7.31. I tried to create a dummy dataset with the same structure to replicate the problem with the need of sending my dataset. However all of them gave identical() results between 32-bit and 64-bit. Note that coef()$fRow is a 1266 x 6 data.frame. Is it cor

Re: [R] Save the result of map.market function to HTML file?

2015-06-04 Thread valerio orfano
Belowe My code: treemap.R library(portfolio) library(sendmailR) setwd("C:/Users/Administrator/Desktop") data_all <- read.csv("prova_data.txt", sep='\t', stringsAsFactors = FALSE) distinctstorage <- unique(data_all$Storage) stor_var <- "Storage001" filename <-