Re: [R] Problem with r project in ubuntu xenial

> On Nov 7, 2017, at 3:46 PM, George Balas wrote: > > For anyone who sees this conversation. > > There is a bug in installation of igraph in R language in Ubuntu. There is > a solution in stackoverflow. We have to use the devtools. Write this code: >

Re: [R] Problem with r project in ubuntu xenial

For anyone who sees this conversation. There is a bug in installation of igraph in R language in Ubuntu. There is a solution in stackoverflow. We have to use the devtools. Write this code: install.packages("devtools") library(devtools) install_github("igraph/rigraph") If there are errors

Re: [R] Missing information in source()

Hello, Try print(head(...)) Hope this helps, Rui Barradas Em 07-11-2017 20:01, Tom Backer Johnsen escreveu: Dear R-help, I am running a Mac under Sierra, with R version 3.4.2 and RStudio 1.1.383. When running head () or tail () on an object in a script using source

[R] Error when attempting to see "Corpus" metadata

R Project I receive the erro highlighted in yellow when attempting to combine two Corpus, so I'm assuming I'm not combining the two variables (nb_pos and nb_neg) in the following line nb_all <- c(nb_pos,nb_neg,recursive=TRUE) # anyone see anything wrong with this line of code

[R] Revolutions blog: October 2017 roundup

Since 2008, Microsoft (formerly Revolution Analytics) staff and guests have written about R every weekday at the Revolutions blog (http://blog.revolutionanalytics.com) and every month I post a summary of articles from the previous month of particular interest to readers of r-help. In case you

Re: [R] Pathview xml issue

Between the lack of a reproducible example and the question being about a Bioconductor package, you probably won't see much response here. Try the Bioconductor mailing list? -- Sent from my phone. Please excuse my brevity. On November 7, 2017 11:15:51 AM PST, Yin wrote:

Re: [R] Missing information in source()

On 11/7/2017 12:01 PM, Tom Backer Johnsen wrote: Dear R-help, I am running a Mac under Sierra, with R version 3.4.2 and RStudio 1.1.383. When running head () or tail () on an object in a script using source

Re: [R] Missing information in source()

On 07/11/2017 3:01 PM, Tom Backer Johnsen wrote: Dear R-help, I am running a Mac under Sierra, with R version 3.4.2 and RStudio 1.1.383. When running head () or tail () on an object in a script using source

Re: [R] Missing information in source()

Either change your script by adding print(...) calls where you want to see something printed or change the call to source() by adding print.eval=TRUE or echo=TRUE. E.g., > cat(file = tf <- tempfile(), "head(10:1)\ntail(letters)\nx <- gamma(0:4)\nx\n") > source(tf) Warning message: In gamma(0:4)

[R] Missing information in source()

Dear R-help, I am running a Mac under Sierra, with R version 3.4.2 and RStudio 1.1.383. When running head () or tail () on an object in a script using source

Re: [R] fill histogram in ggplot

Hi Elahe, You pass 'probable' to the fill aesthetic along the lines of; ggplot(hist) + aes(x=mms, fill = probable) + geom_histogram(binwidth=1) The NAs might give you three and not two colours. I'm guessing you want distinct colours. In this case 'probable' should be a factor and not an

[R] Pathview xml issue

Hi, I'm using GAGE/pathview to analyze my RNA-seq and phospho-protein data. The following error occurs after this command line below: >pv.out.list <- sapply(path.ids2[1:3], function(pid) pathview( gene.data = cnts.d, pathway.id = pid, gene.idtype="SYMBOL",kegg.native = F, same.layer = T,

[R] fill histogram in ggplot

Hi all, I have the following data and I have a histogram for mms like ggplot(hist,aes(x=hist$mms))+ geom_histogram(binwidth=1,fill="white",color="black")and then I want to fill the color of histogram by probable=1 and probable=0, could anyone help me in this? My data:

Re: [R] Fitdistrplus and Custom Probability Density

> On Nov 7, 2017, at 6:58 AM, Eric Berger wrote: > > Why not define your own functions based on d? > e.g. > myCumDist <- function(x) { integrate(d, lower=-Inf, upper=x)$value } > myQuantile <- function(x) { uniroot(f=function(y) { h(y) - x }, > interval=c(-5,5)) } #

[R] Survfit when new data has only 1 row of data

Dear R-help, I am using R version 3.4.0 within Windows, and survival 2.41-3. I have fit a Prentice Williams and Peterson-Counting Process model to my data as shown below. This is basically an extension of the Cox model for interval censored data. My dataset, bdat5 can be found here:

Re: [R] Fitdistrplus and Custom Probability Density

Why not define your own functions based on d? e.g. myCumDist <- function(x) { integrate(d, lower=-Inf, upper=x)$value } myQuantile <- function(x) { uniroot(f=function(y) { h(y) - x }, interval=c(-5,5)) } # limits -5,5 should be replaced by your own which might require some fiddling e.g. d <-

[R] Fitdistrplus and Custom Probability Density

Dear All, Apologies for not providing a reproducible example, but if I could, then I would be able to answer myself my question. Essentially, I am trying to fit a very complicated custom probability distribution to some data. Fitdistrplus does in principle everything which I need, but if require

Re: [R] Using MLE on a somewhat unusual likelihood function

> On 7 Nov 2017, at 11:38 , peter dalgaard wrote: > > Presumably you wanted -sum(lnagevec[lnwagevec>-200 & ln2wagevec<200], > na.rm=TRUE) Well, after correcting the obvious spelling error. Sorry about that... -pd -- Peter Dalgaard, Professor, Center for Statistics,

Re: [R] Using MLE on a somewhat unusual likelihood function

(inline) > On 7 Nov 2017, at 06:02 , Hugo André wrote: > > So I am trying to use the mle command (from stats4 package) to estimate a > number of parameters using data but it keeps throwing up this error message: > > Error in solve.default(oout$hessian) : > Lapack

[R] Fwd: FW: Time Series

[Please send replies to r-help, not individual responders] Emre, In R, when you call a function defined via something like f <- function( foo, bar ) then you can call it as, for, example a <- f(x,y) or a <- f(foo=x, bar=y) or even a <- f( bar=y, foo=x) # notice I switched the order! The

Re: [R] FW: Time Series

Following Erin's pointer: library(zoo) times <- seq(from=as.POSIXct("2015-12-18 00:00:00"), to=as.POSIXct("2017-10-24 23:00:00"), by="hour") mydata <- rnorm(length(times)) tseri <- zoo( x=mydata, order.by=times ) HTH, Eric On Tue, Nov 7, 2017 at 9:59 AM, Erin Hodgess

Re: [R] FW: Time Series

Hello! What is the error message, please? At first glance, you are using the "ts" function. That doesn't work for hourly frequency. You may want to create a zoo object. This is Round One. Sincerely, Erin On Tue, Nov 7, 2017 at 1:46 AM, Emre Karagülle wrote: > >