[R] install package av on centos 7

2019-12-09 Thread Schneider, Dominik
I tried to install the av package on centos7 server with R 3.5.3 but am receiving an error. I installed ffmpeg and ffmpeg-devel on centos 7 as below and then followed the instructions of specifying LIB_DIR and INCLUDE_DIR for installing av. any help would be appreciated. Thanks ``` >

Re: [R] Having problems with the ifelse and negative numbers

2019-12-09 Thread Richard M. Heiberger
or even simpler sqrt(abs(A)) On Mon, Dec 9, 2019 at 8:45 AM Eric Berger wrote: > > Hi Bob, > You wrote "the following error message" - > when in fact it is a Warning and not an error message. I think your > code does what you hoped it would do, in the sense it successfully > calculates the

Re: [R] Where is the SD in output of glm with Gaussian distribution

2019-12-09 Thread Fox, John
Dear Bert, It's perhaps a bit pedantic to point it out, but the dispersion is estimated from the Pearson statistic (sum of squared residuals or weighted squared residuals) not from the residual deviance. You can see this in the code for summary.glm(). Best, John

Re: [R] how to create a txt file with parsed columns

2019-12-09 Thread Ana Marija
Thanks for getting back to me, I resolved my problem with this: library(reshape2) c=dcast(a, rs ~ GENE) d=merge(c,b,by="rs") d[is.na(d)] <- 0 On Sun, Dec 8, 2019 at 11:03 PM Jim Lemon wrote: > > Hi Ana, > Is this what you want? > > a<-read.table(text="GENErs BETA > 1

Re: [R] Where is the SD in output of glm with Gaussian distribution

2019-12-09 Thread Marc Girondot via R-help
Le 09/12/2019 à 16:45, Bert Gunter a écrit : > In addition, as John's included output shows, only 1 parameter, the > intercept, is fit. As he also said, the sd is estimated from the > residual deviance -- it is not a model parameter. > > Suggest you spend some time with a glm tutorial/text. I

Re: [R] Where is the SD in output of glm with Gaussian distribution

2019-12-09 Thread Bert Gunter
In addition, as John's included output shows, only 1 parameter, the intercept, is fit. As he also said, the sd is estimated from the residual deviance -- it is not a model parameter. Suggest you spend some time with a glm tutorial/text. Bert On Mon, Dec 9, 2019 at 7:17 AM Marc Girondot via

Re: [R] Where is the SD in output of glm with Gaussian distribution

2019-12-09 Thread Eric Berger
summary(gnul) shows the std error of the coefficient estimate On Mon, Dec 9, 2019 at 5:16 PM Marc Girondot via R-help wrote: > > Let do a simple glm: > > > y=rnorm(100) > > gnul <- glm(y ~ 1) > > gnul$coefficients > (Intercept) >0.1399966 > > The logLik shows the fit of two parameters

Re: [R] Where is the SD in output of glm with Gaussian distribution

2019-12-09 Thread Fox, John
Dear Marc, For your simple model, the standard deviation of y is the square-root of the estimated dispersion parameter: > set.seed(123) > y <- rnorm(100) > gnul <- glm(y ~ 1) > summary(gnul) Call: glm(formula = y ~ 1) Deviance Residuals: Min1QMedian3Q Max

Re: [R] A better scales::dollar() ?

2019-12-09 Thread Hadley Wickham
You can get pretty close with label_number_si(): pounds <- scales::label_number_si(prefix = "£") pounds(10 ^ (0:7)) #> [1] "£1""£10" "£100" "£1K" "£10K" "£100K" "£1M" "£10M" Created on 2019-12-09 by the [reprex package](https://reprex.tidyverse.org) (v0.3.0.9001) Hadley On Fri, Dec

[R] Where is the SD in output of glm with Gaussian distribution

2019-12-09 Thread Marc Girondot via R-help
Let do a simple glm: > y=rnorm(100) > gnul <- glm(y ~ 1) > gnul$coefficients (Intercept)   0.1399966 The logLik shows the fit of two parameters (DF=2) (intercept) and sd > logLik(gnul) 'log Lik.' -138.7902 (df=2) But where is the sd term in the glm object? If I do the same with optim, I can

Re: [R] Having problems with the ifelse and negative numbers

2019-12-09 Thread Eric Berger
Hi Bob, You wrote "the following error message" - when in fact it is a Warning and not an error message. I think your code does what you hoped it would do, in the sense it successfully calculates the sqrt(abs(negativeNumber)), where appropriate. If you want to run the code without seeing this

Re: [R] Having problems with the ifelse and negative numbers

2019-12-09 Thread Kevin Thorpe
The sqrt(-A) is evaluated for all A. The result returned is conditional on the first argument but the other two arguments are evaluated on the entire vector. Kevin -- Kevin E. Thorpe Head of Biostatistics, Applied Health Research Centre (AHRC) Li Ka Shing Knowledge Institute of St. Michael's

[R] Having problems with the ifelse and negative numbers

2019-12-09 Thread rsherry8
Please consider the following two R statements: A = runif(20, min=-1,max=1) ifelse( A < 0, sqrt(-A), A ) The second statement produces the following error message: rt(-A) : NaNs produced I understand that you cannot take the square root of a negative number but I thought the

[R] [R-pkgs] concurve v 2.3.0 - Comparing Functions, Bootstrapping, and Exporting Tables

2019-12-09 Thread Zad Chow
Pleased to announce that the next version (2.3.0) of our package “concurve” is out on CRAN (https://cran.r-project.org/package=concurve). In addition to plotting confidence (consonance) curves, it can now plot two functions next to one another to see the amount of overlap. It can also plot