similar to: Survival analysis MLE gives NA or enormous standard errors

Hi, I am running a Cox Mixed Effects Hazard model using the library coxme. I am trying to model time to onset (age_sym1) of thought problems (e.g. hearing voices) (sym1). As I have siblings in my dataset, I have decided to account for this by including a random effect for family (famid). My covariate of interest is Mother's diagnosis where a 0 is bipolar, 1 is control, and 2 is major

Hi I am running the following code: sym <- c(sym1,sym2,sym4) lifedxm <- c("O-BD","O-WELL","O-UNI") life <- c(lifedxm,lifedxm,lifedxm) tp <- c("TP-ANY","TP-ANY", "TP-ANY", "TP-SUB", "TP-SUB", "TP-SUB", "TP-CLIN" , "TP-CLIN", "TP-CLIN") data <-

Hi, I am running the following model: fit1.full <- coxme(Surv(age_sym1, sym1) ~ sex + lifedxm*sex + (1|famid), data=bip.surv) I would like to extract the AIC from that object to calculate the AICC. However, when I look at str(fit1.full) and summary(fit1.full) (pasted below) I don't see anything that would allow me to get pull the AIC out from that object. Is there a way to retrieve the

Add a helper to initialize a bvec based of a page pointer. This will help removing various open code bvec initializations. Signed-off-by: Christoph Hellwig <hch at lst.de> --- block/bio-integrity.c | 7 +------ block/bio.c | 12 ++---------- include/linux/bvec.h | 15 +++++++++++++++ 3 files changed, 18 insertions(+), 16 deletions(-) diff --git a/block/bio-integrity.c

Dear all, I tried once to create one variable called bip such that: bip <- cip + (1/f(cip))*fi(f,cip) And this was working. But now, doing the same thing I did before, the software shows me the following message: Error in f(x, ...) : unused argument(s) (subdivision = 2000) I have the variable cip and the variable bip should be created such that: Fn <- ecdf(cip) f <- function(x) {(1 -

Hi, I want to calculate the t- and p-values for a linear model using the Newey West estimator. I tried this Code and it usually worked just fine: > oberlm <- lm(DYH ~ BIP + Infl + EOil, data=HU_H) > coeftest(oberlm, NeweyWest(oberlm, lag=2)) t test of coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.1509950 0.0743832 2.0300 0.179486 BIP

Hi All, I have a problem to create a variable that is a function of an integral of another function. The problem is the following: I have a variable called cip. I have to create another variable called bip that is a function of the former variable cip and also the cumulative distribution function of the variable cip. Consider the equation below:

At the same site, DTMF recognition is functioning badly, sometimes duplicating digits and sometimes totally missing others. We have checked already /proc/interrups, there is no interrupt being shared. Our zaptel has .. span=1,1,0,ccs,hdb3 On zapata we have ... switchtype=euroisdn signalling=pri_cpe relaxdtmf=no (yes doesn't seem to help) -----Original Message----- From:

* Please cc me if you reply as I am a digest subscriber * Hi, I am wondering how I can run a multilevel survival model in R? Below is some of my data. > head(bi0.test) childid famid lifedxm sex age delta 1 22.02 22 CONTROL MALES 21.36893 0 2 13.02 13 MAJOR MALES 21.18001 0 3 64.02 64 CONTROL MALES 20.09377 0 4 5.02 5 CONTROL FEMALES

Hello All, I have to create a variable that is a function of another one (already created), its cumulative distribution function and the integral of this cumulative distribution, with limits: 0 and the value of the variable. To be clear, I have the variable called “cip”. And its cdf called “cdfcip” I need to create the variable: bip = cip + ((1 – cdfcip)^4)*integral((1-cdf(u))^4*du, 0, value

Hi I have a problem with Asterisk 11.5.1. When I pick up an incoming phone call sometimes I need to transfer to someone else in the organization. I then dial a number on my phone, and press Xfer. Sometimes it works well, I mean, the number I dialed get the call and can chat with correspondent. Sometimes, the number I dialed get the communication and while he chats with correspondent, gets a bip

Hello all, I've played with background and play sounds apps and googled around and asked the list before to no avail. Does anyone know of a way to have tones played during the call progress stage of the call? We (especially on some international circuits) get up to 5 seconds of silence before the phone starts ringing or is busy. I don't want to force "R" on the Dial app as

Hello, I have posted this problem before, but thought I try to explain it a bit better. I'm using the function plm to create a fixed effects model for panel data, my method is therefor "within" my effect is "twoways". My Data contains unbalanced Panels due to missing Values, but contains 309 observation for 11 variables (incl. response), with no missing Values. These 309

Hello list, Hope you all doing fine! I've tried to use the 'alias' directive in the indications.conf file but apparently it doesn't work.... It looks like maybe this feature was removed, because old sample for the indications.conf file have example using the alias parameter, but newer samples don't have it anymore.... also I couldn't find any ticket saying this parameter

Dear List, I am having trouble with some R code I have written to perform Redundancy Analysis (RDA) on a matrix of species abundance data (Y) and a matrix of environmental data (X). RDA is a constrained form of PCA and can be thought of as a PCA of the fitted values of a regression of each variable in Y on all variables in X. For info, the first use of RDA is in: Rao, C.R, 1964. The use and

Hi everybody, I'm testing asterisk@home 0.4, looks great so far I was working when I have been alerted by a bip comming from the * pc... I connected a screen to it and saw that there was a message which looked like : Message from syslogd@asterisk1 at Thu Feb 10 09:01:00 2005 ... asterisk1 so I stopped asterisk, type mail and got a strange mail saying that user xxxx@yahoo.com could

Hi All, We are just starting to deploy a new PRI IVR system, and the incoming calls sometimes get random short 'bips' while navigating our IVR menu. Any hint on what this can be? Best regards, PauloHM

Hello, I have recently had 'Broadband' (Cable) Internet installed in my house and I have setup the my Linux PC to do Network Address Translation for my Windows PC. This system has worked well with the PPP modem, and I know that a full Proxy server would be better, but at the moment I haven't got the time to set one up etc. This setup worked fine: Modem[ppp0

Starting with 3.7 rc1, my workstation seems to loose ram. Up until (and including) 3.6, used-(buffers+cached) was roughly the same as sum(rss) (taking shared into account). Now there is an approx 6G gap. When the box first starts, it is clearly less swappy than with <= 3.6; I can''t tell whether that is related. The reduced swappiness persists. It seems to get worse when I update