search for: sex1

Dear R-project! How do i create 1 dummy from 2 already existing dummys. To be more precise, I want to create a dummy from a dummy called "sex" and another called "sex1" when both thoose dummys are 1 I want my created dummy "samesex" to take 1. Thanks for the help! Paulie [[alternative HTML version deleted]]

...ogramming problem for which I have found a solution, but I feel it is a rather poor one. I wonder if there's some other (more clever) solution, using (maybe?) vectorization or subscripting. A toy example would be: rel1 rel2 rel3 age0 age1 age2 age3 sex0 sex1 sex2 sex3 1 3 NA 25 23 2 NA 1 2 1 NA 4 1 3 35 67 34 10 2 2 1 2 1 4 4 39 40 59 60 1...

...message: In coxph(Surv(tstart, time, status) ~ SEX:strata(tgroup), data = data) : X matrix deemed to be singular; variable 2 4 6 > fit1 Call: coxph(formula = Surv(tstart, time, status) ~ SEX:strata(tgroup), data = data) coef exp(coef) se(coef) z p SEX1:strata(tgroup)tgroup=1 -0.0715 0.9310 0.0524 -1.36 0.17 SEX2:strata(tgroup)tgroup=1 NA NA 0.0000 NA NA SEX1:strata(tgroup)tgroup=2 -0.2624 0.7692 0.0655 -4.01 6.1e-05 SEX2:strata(tgroup)tgroup=2 NA NA 0.0000 NA NA SEX1:strata(tgroup)tgroup=3...

...Binomial distributed data. The main challenge am facing is the fact that the residual variance depends on the mean (if I am correct). I strongly feel that the script below may be wrong due to mean-variance relationship #####R code####### fit1s <-glm(mydata$vzv~mydata$age.c+mydata$age2+mydata$sex1, family=Gamma(link=log)) x.betahat1<-fit1s$fitted.values res1<-fit1s$residuals b<-1000 for (i in 1:b){ b.i <- sample(index, size=n, replace=T) res.star1=res1[b.i] bst1=x.betahat1+res.star1 mydata1 <-data.frame(age,age2,sex,bst1) ########Modeling ###########...

Hi! I have a dataframe like this: dat=data.frame(Age=c(rep(30,8),rep(40,8),rep(50,8)),Period=rep(seq(2005,2008,1),3),Rate=c(seq(1,8,1),seq(9,16,1),seq(17,24,1)),Sex=rep(c(rep(0,4),rep(1,4)),3))attach(dat)dat    Age Period Rate Sex1   30   2005    1   02   30   2006    2   03   30   2007    3   04   30   2008    4   05   30   2005    5   16   30   2006    6   17   30   2007    7   18   30   2008    8   19   40   2005    9   010  40   2006   10   011  40   2007   11   012  40   2008   12   013  40   2005   13   114  40   2006  ...

...mputation results: with.imputationList(miset, coxph(Surv(time, status) ~ age + sex + hepmeg + platelet + trt + trig)) MIcombine.default(mifit) results se (lower upper) missInfo age 0.035548792 0.0082506946 0.019373545 0.0517240397 4 % sex1 -0.070760613 0.2563372831 -0.580309741 0.4387885156 34 % hepmeg1 0.932378808 0.2026274576 0.532555416 1.3322021993 23 % platelet -0.001757899 0.0009480636 -0.003620446 0.0001046469 14 % trt2 0.137413162 0.1924230007 -0.243815288 0.5186416117 29 % trig 0.003979287 0....

Hi all, I have a question on the package « survey” I have some difficulties to use the function ‘calibrate’. Although it works well with one single factor variable, I cannot use it for 2 and get the message “Erreur dans regcalibrate.survey.design2(design, formula, population, aggregate.stage = aggregate.stage, : Population and sample totals are not the same length.” Here is the format

Hello, I'm using aov() to analyse changes in brain volume between males and females. For every subject (there are 331 in total) I have 8 volume measurements (4 different brain lobes and 2 different tissues (grey/white matter)). The data looks like this: Subject Sex Lobe Tissue Volume subect1 1 F g 262374 subect1 1 F w 173758 subect1 1 O g 67155 subect1 1 O w 30067 subect1 1 P g 117981