similar to: Unbalanced parentheses printed by warnings() crash text editor

Hello, I have a big dataframe with *NO* na's (9 columns, 293380 rows). # doing memory.limit(size = 1000000000) cor(x) #gives Error in cor(x) : missing observations in cov/cor In addition: Warning message: NAs introduced by coercion #I found the obvious workaround: COR <- matrix(rep(0, 81),9,9) for (i in 1:9) for (j in 1:9) {if (i>j) COR[i,j] <- cor (x[,i],x[,j])} #which works fine,

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 + ......... >

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

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(

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

*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 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

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?

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

I have configured iax.conf and extensions.conf as instructed on FWD website (http://www.freeworlddialup.com/help/?p=knowledgebase&c=18&a=76) and I can successfully receive calls and make test calls to 612, 613, etc. The problem is that I can not make a call to another FWD user. Here is what asterisk says: -- Executing [393xxxxxx@default:1] Set("Zap/1-1",

I've been running a couple Centos5 and RHEL4/5 servers with samba for a while now and everything has been working great with our Windows 2003 AD. All of a sudden though I'm experience something really weird on one of the RHEL5 boxes. Whenever I try to connect as a Domain Admin from one particular Vista client, I get access denied and repeated prompts for a username/password - this has

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

we use samba as member servers with security=domain pointing towards a WinNT4 PDC. accounts that have workstation logon restrictions are unable to connect to samba shared ressources but are able to connect to windows nt shared ressources as they should. normally, workstation logon restriction should deny only interactive logons. the 2.0.5 logfile entry:

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: > > >

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: >>>>>

Hi all- Equipment: Xlite softphone Polycom 301 with SIP 2.1.1 and BootROM 3.2.3 Polycom 501 with SIP 2.1.1 and BootROM 3.2.3 Asterisk 1.4.2 SIP Trunk to FWD I wanted to post this problem as I haven't found it described in any of the past presence threads on here. I use an identical configs for a Polycom 501 and 301. (I actually unplug one when the other is in use). The

Hi all, who knows how to get an ESRI Shapefile for the NUTS-2 Regions of the enlarged European Union? Particularly I want to draw maps of Germany, Poland, Czech Republik, Hungary and Austria. I've found Shapefiles for the US, Russia and other countries elsewhere in the web, but for Europe it seems really difficult. With many thanks for any hint Albrecht

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz<-sapply(ma, myfunction) To my surprise, It almost costs me the same time to get my results if I directly use a loop

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz<-sapply(ma, myfunction) To my surprise, It almost costs me the same time to get my results if I directly use a loop