similar to: Multilevel survival model

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

Hi, I am trying to fit the following model: sr.reg.s4.nore <- survreg(Surv(age_sym4,sym4), as.factor(lifedxm), data=bip.surv) Where age_sym4 is the age that a subject develops clinical thought problems; sym4 is whether they develop clinical thoughts problems (0 or 1); and lifedxm is mother's diagnosis: BIPOLAR, MAJOR DEPRESSION, or CONTROL. I am interested in whether or not

I have the following code: ----------------------------------------------------------------------------------------------------- result <- matrix(NA, nrow=1, ncol=5) for(i in 1:(nsnp-1)) { for(j in (i+1):nsnp){ tempsnp1 <- data.lme[,i] tempsnp2 <- data.lme[,j] fm1 <- lme(trait~sex+age+rmtemp.b+fc+tempsnp1+tempsnp2+tempsnp1*tempsnp2, random=~1|famid, na.action=na.omit) fm2 <-

I am involved in a study where, as in most of life, men demonstrate themselves to be recalcitrant. So while we have many probands and most of their mothers we only have about 50% of the trios being complete. I have been running tdt and trio.types. It appears as if it is ignoring the duos. Sometimes a duo can be informative. For instance Father ..missing Mother 1/2 Proband 1/1 This duo shows that

I am using dgc.genetics to perform TDT analysis on SNP data from a cohort of trios. I now have a file with about 6008 variables. The first few variables related to the pedigree data such as the pedigree ID the person ID etc. Thereafter each variable is a specific locus or marker. The variables are named by a pattern such as "Genotype.nnnnn" with nnnnn corresponding to a number which

Just some clarification here. On Monday, June 30, 2003 12:04 AM, Steven Critchfield [SMTP:critch@basesys.com] wrote: > > As for slips and bipolar violations... > T1s are just high speed serial lines. A sleep is when you loose sync > with the far side and when you see a 1 come across the line, you may > not > know which bit it was for. A slip is when one end of the circuit is

Hi, I am fitting a GEE model using gee R package, but I am not sure how to get OR and its CI? Could anyone give me some hints? Here are some output: > gee.obj <- gee(Affection~Sibsex+Probandsex,id = FAMID,family = binomial,corstr = "independence",data =seldata) Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate

Hello, I have a pedigree file such this: FAMID ID FA ID MO ID SEX STATUS 1 1 0 0 2 0 1 2 3 1 2 2 1 4 3 1 1 2 1 5 3 1 1 0 1 6 3 1 1 0 1 7 3 1 2 2 1 8 3 1 1 0 1 9 3 1 1 0 1 10 3 1 2 0 1 11 3 1 1 0 1 3 0 0 1 0 2 12 13

I'm emailing this as the customer in this case, since my carrier appears to be completely unable to solve this. A brief rundown of the problem: - We have several voice lines going through the ADIT, of course into a VoIP type of arrangement. - Voice traffic will become choppy, even drop calls completely, at random but quite often. Modem calls (which we needed unfortunately) are a joke. -

I want to classify bipolar neurons in human cochleas and have data of the following structure: Vol_Nuc Vol_Soma 1 186.23 731.96 2 204.58 4370.96 3 539.98 7344.86 4 477.71 6939.28 5 421.22 5588.53 6 276.61 1017.05 7 392.28 6392.32 8 424.43 6190.13 9 256.41 3850.51 10 249.17 3118.14 11 276.97 3037.29 12 295.30 3703.76 13 314.43 5265.97 14 301.15 5781.73 I

Hello, I am missing any description of zaptel/DAHDI alarms. The TE200 series user manual contains only a description of LEDs states. These alarms states are visible in zttool/dahditool or in astersick CLI (zap show status) and I wonder what is the real meaning of these alarms for E1 channel. Possible alarm states (based on zaptel.h 1.2): 1. No alarms 2. Recovering from alarm 3. In loopback

Hi, I'm using kinship2 to calculate heritabilty, but I would like calculate in GenABEL too. I trying the code: > require(kinship2) > require(GenABEL) > pedig = with(Dados, pedigree(id=IID, dadid=PAT, momid=MAT, sex=SEX, famid=FID, missid=0)) > kmat = kinship(pedig) > (mod1 = polygenic(altura ~ SEX + idade, data=Dados, kin=kmat)) Erro em intI(i, n = d[1],

Hello everyone, I bought 2 HFC-ISDN Cards and want to run the first card in NT-Mode an the second one in TE-Mode. Everything looks ok under SuSE 9.1, but I can't dial out. I removed one card, for testing purposes and want to run this one card in TE-Mode. I only want to make a call with my Grandstream BT-101 over Asterisk via ISDN. When I try to make a call I get: ---------------------

Hi Im using the kinship package to draw a pedigree. On my data set this works fine but when i add indivudals to the pedigree i keep getting an error i hope someone can help me! This is the code im using: Data<-read.table("Tree.txt", header=T, sep=",") attach(Data) ped<-pedigree(id, dadid, momid, sex, aff) par(xpd=T) plot.pedigree(ped) This is my data looks like

Hi! My question(s) in the end might be silly but I am no expert on this, so here it goes: Noy-Meir (1973), Pielou (1984) and a few others have pointed to non-centered PCA being in some cases useful. They clearly explain that "it is the case" when multi-dimensional data display distinct clusters (which have zero, or near-zero, projections in some subset of the axes) and the task is

Hi, I have been trying to get a Wildcard TE122 card running here the last couple of days. libpri and zaptel are all installed and configured to E1 specs. The jumper on the card is on, so configured for E1 (I'm in Norway). When running zttool, I get 'Alarms: RED' on the single card installed. Other information in zttool: Current alarms: Red Alarm Sync source: Internally clocked

Excuse my ignorance (I just got crap on the #centos IRC channel for this question), but is there a (easy!) way to have and IRC and/or Jabber server relay a login to a Microsoft Active Directory server for authentication? If there's a better question to ask this question, please point me in that direction, and I'll be happy to do so Thanks -------------- next part -------------- An HTML

I have 7965 and am trying to convert the firmware to SIP (SIP45.8-3-4SR1S). Does anyone have a valid XMLDefault.cnf.xml they could share? I have tried the version at voip-info<info.org/wiki/view/Asterisk+phone+cisco+79x1+xml+configuration+files+for+SIP&view_comment_id=14768#Troubleshooting>for the 7941/7961 but unfortunately /var/log/messages shows in.tftp stops sending after

Hi All; Suddenly the zaptel channel look like stop working and it is giving me this error when I do zap restart: [May 24 19:30:21] ERROR[2772]: chan_zap.c:7219 mkintf: Unable to open channel 1: Device or resource busy here = 0, tmp->channel = 1, channel = 1 [May 24 19:30:21] ERROR[2772]: chan_zap.c:10582 build_channels: Unable to register channel '1,2' To what it could be related?