similar to: cor() fails with big dataframe

Hello everyone, Hope it is the good place for this (I discuss the question of the right place below). Most of the time, warnings are more than 1000 characters long and thus are truncated. Most of the time, this generates printouts with unbalanced parentheses. Intelligent text editors which do parentheses highlighting get very confused with this. After too many warnings, they give errors, and

Hi everyone. There's been several threads on baseline hazard in Cox model but I think they were all on cumulative baseline hazard, for instance http://tolstoy.newcastle.edu.au/R/help/01a/0464.html http://tolstoy.newcastle.edu.au/R/help/01a/0436.html "basehaz" in package survival seems to do a cumulative hazard. extract from the basehaz function: sfit <- survfit(fit) H

Hi, Can anybody tell me what the message below means and how to overcome it. Thanks, Neil Warning message: X matrix deemed to be singular; variable 2 in: coxph(Surv(age_at_death, death) ~ project$pluralgp + project$yrborn + ......... >

Hello everyone, I am searching for a covariate selection procedure in a cox model formulated as a counting process. I use intervals, my formula looks like coxph(Surv(start,stop,status)~ x1+x2+...+cluster(id),robust=T) where id is a country code (I study occurence of civil wars from 1962 to 1997). I'd like something not based on p-values, since they have several flaws for this purpose. I turned

Hello, coxph does not use any information that are in the dataset between event times (or "death times") , since computation only occurs at event times. For instance, removing observations when there is no event at that time in the whole dataset does not change the results: > set.seed(1) > data <- as.data.frame(cbind(start=c(1:5,1:5,1:4),stop=c(2:6,2:6,2:5),status=c(rep(

*Paid internship "web application development with R"* The main skills for this proposed paid internship are: R, Shiny, and components used by Shiny (plotly, javascript). Other skills (Debian, Docker, CKAN...) are useful. This internship offer is part of a partnership between the Grenoble Alpes Cybersecurity Institute and the CYBIS project. The *Cybersecurity Institute* aims to

Hi, I would like to use constrained splines in a GLM model (Poisson link) to take into account the nonlinear effect of some covariate. The constraints I need are described below. I have several variables that I need concurrently in the same model. I looked at package mgcv but I do not know if/how I can use it in GLM (not GAM) : I could not manage to adapt the mono.con(mgcv) example to GLM. The

Dear R helpers, this is a question connex to the subject of reading large data-sets - I have a large db (MySQL) ~80'000 records for each year / 120 variates from which I extract relevant info to be treated in R (no problem) Now each record in my db represents an experimental unit; moreover I will get - in a few days- a set of separate files(one for each records....~80'000 files pro

Hi, I have been dreaming about a complete R environment on my USB stick for a long time. Now I finally want to realize it but what I am missing is a good, portable editor for R which has tabs and syntax highlighting, can execute code, has bookmarks and a little project file management facility pretty much like Tinn-R has those. I like Tinn-R but it seems like there is only a very old version of

Am Mi, 1. Nov 2017, um 17:40, schrieb Duncan Murdoch: > On 01/11/2017 12:02 PM, Albrecht Kauffmann wrote: > > Dear Duncan, > > > > Many thanks! > > > > Am Mi, 1. Nov 2017, um 16:17, schrieb Duncan Murdoch: > >> On 01/11/2017 8:11 AM, Albrecht Kauffmann wrote: > >>> Hi all, > >>> > >>> after compiling R-3.4.2 on

On 01/11/2017 1:08 PM, Albrecht Kauffmann wrote: > Am Mi, 1. Nov 2017, um 17:40, schrieb Duncan Murdoch: >> On 01/11/2017 12:02 PM, Albrecht Kauffmann wrote: >>> Dear Duncan, >>> >>> Many thanks! >>> >>> Am Mi, 1. Nov 2017, um 16:17, schrieb Duncan Murdoch: >>>> On 01/11/2017 8:11 AM, Albrecht Kauffmann wrote: >>>>>

DeaRs, I have a doubt about: BIC (Bayesian Information Criterion) SBIC (Schwartz Bayesian Informarion Criterion) SIC (Schwartz Information Criterion) In many references these are know as the same (eg. stepAIC() function) but I just found a SAS8.2 output that show either the BIC and SIC values for a logistic regression.. simillary values but different. 1) question: What are the differences?

On 01/11/2017 12:02 PM, Albrecht Kauffmann wrote: > Dear Duncan, > > Many thanks! > > Am Mi, 1. Nov 2017, um 16:17, schrieb Duncan Murdoch: >> On 01/11/2017 8:11 AM, Albrecht Kauffmann wrote: >>> Hi all, >>> >>> after compiling R-3.4.2 on opensuse leap 42.3, make check failed. Until >>> R-3.4.1 I never had a problem with these tests. No,

Dear Duncan, Many thanks! Am Mi, 1. Nov 2017, um 16:17, schrieb Duncan Murdoch: > On 01/11/2017 8:11 AM, Albrecht Kauffmann wrote: > > Hi all, > > > > after compiling R-3.4.2 on opensuse leap 42.3, make check failed. Until > > R-3.4.1 I never had a problem with these tests. No, the programm stops > > at the following line of tests/reg-tests-1d.R: > > >

Dear Duncan, thank you very much, that is it! Indeed, there was activated in my internet router the "telekom navigation help". When I am (or the system is) looking for an not existent url, this "help" opens its own website. After deactivation of telekom navigation help, the normal message (server not found) comes. I never gave attention to this "help", although it is

I've found two problems in interpreting adv-r-master/book/build-book.r: 1. All pathes in build-book.r refer to the starting-directory "adv-r-master". However, the script build-book.r is located in the directory "book", which is located in directory "adv-r-master". Therefore, pathes starting at "." are wrong (should start at ".."). 2. Some

Hi all, I'm faced with a problem applying the sp package: The projection argument in readShapePoly(Shapefile,proj4string="CRS class argument") e.g.: CRS("+proj=aea +lat_1=46 +lat_2=73 +lat_0=60 +lon_0=84 +x_0=0 +y_0=0 +ellps=clrk66 +units=m +no_defs") doesn't have any impact on the plotted object. I also tested the simple example: xy = cbind(x = 2 * runif(100) -

Hi all, after compiling R-3.4.2 on opensuse leap 42.3, make check failed. Until R-3.4.1 I never had a problem with these tests. No, the programm stops at the following line of tests/reg-tests-1d.R: > ## available.packages() (not) caching in case of errors > tools::assertWarning(ap1 <- available.packages(repos = "http://foo.bar")) > tools::assertWarning(ap2 <-

Dear all, I am interested in plotting maps visualizing spatial statistics in an aggregated fashion, according to administrative boundaries. Partially, I want to visualize some spatial data for administrative units (autonomous republics, oblasts, krays) of the Russian Federation on a geographical map. I have found the maps package (and related) and would like to use this package e.g. in a kind

On 01/11/2017 8:11 AM, Albrecht Kauffmann wrote: > Hi all, > > after compiling R-3.4.2 on opensuse leap 42.3, make check failed. Until > R-3.4.1 I never had a problem with these tests. No, the programm stops > at the following line of tests/reg-tests-1d.R: > >> ## available.packages() (not) caching in case of errors >> tools::assertWarning(ap1 <-