search for: smoking

Displaying 20 results from an estimated 571 matches for "smoking".

2010 Jul 13
0
working out main effect variance when different parameterization is used and interaction term exists
Dear all, Apologies if this question is bit theoretical and for the longish email. I am meta-analyzing the coefficients and standard errors from multiple studies where the raw data is not available. Each study analyst runs a model that includes an interaction term for, say, between sex and smoking and age. Here is an illustrative example example for one study: set.seed(1066) status <- rbinom( 1000, 1, 0.2 ) males <- rbinom( 1000, 1, 0.6 ) smoke <- rbinom( 1000, 1, 0.3 ) age <- runif(1000, min=20, max=80) coef( summary( f1 <- glm( status ~ males*smoke + age...
2011 Dec 19
2
summary vs anova
Hi, I'm sure this is simple, but I haven't been able to find this in TFM, say I have some data in R like this (pasted here: http://pastebin.com/raw.php?i=sjS9Zkup): > head(df) gender age smokes disease Y 1 female 65 ever control 0.18 2 female 77 never control 0.12 3 male 40 state1 0.11 4 female 67 ever control 0.20 5 male 63 ever state1 0.16
2015 Jul 21
2
glm com etiquetas en las variables
Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG <- glm (low ~ X, family = "binomial", data = DATOS) Ejemplo: library(MASS) data(birthwt, package="MASS") birthwt$low <- factor(birthwt$low) birthwt$race <- factor(birthwt$smoke) REG_LOG <- glm (low ~ smoke, family =
2010 Oct 04
1
I have aproblem about nomogram--thank you for your help
...ression by SPSS 17.0. and I want to deveop the nomogram through R-Projject,just like this : > n<-100 > set.seed(10) > T.Grade<-factor(0:3,labels=c("G0", "G1", "G2","G3")) > Sex<-factor(0:1,labels=c("F","M")) > Smoking<-factor(0:1,labels=c("No","yes")) > L<-0.559*as.numeric(T.Grade)-0.896*as.numeric(Smoking)+0.92*as.numeric(Sex)-1.338 > y <- ifelse(runif(n) < plogis(L), 1, 0) > ddist <- datadist(as.numeric(T.Grade,Sex,Smoking)) load package "rms" > ddis...
2010 Oct 04
2
i have aproblem --thank you
...> ddist <- datadist(dfr) > options(datadist='ddist') > n<-100 > set.seed(10) > T.Grade<-factor(0:3,labels=c("G0", "G1", "G2","G3")) > Sex<-factor(sample(0:1, 100, replace=TRUE),labels=c("F","M")) > Smoking<-factor(sample(0:1, 100, replace=TRUE),labels=c("No","yes")) > dfr$L<-with(dfr,0.559*as.numeric(T.Grade)-0.896*as.numeric(Smoking)+0.92*as.numeric(Sex)-1.338) > dfr$y <- with(dfr, ifelse(runif(n) < plogis(L), 1, 0) ) > dfr <- data.frame(T.Grade,Sex,Smoki...
2004 Dec 29
0
GEE with own link function
Hello, I want to fit a GEE with a user-defined link function. For the user-defined link-function I still read http://finzi.psych.upenn.edu/R/Rhelp01/archive/6555.html and http://finzi.psych.upenn.edu/R/Rhelp02a/archive/25727.html. Only for testing purposes I added a new link function (corlogit) in make.link (as well as in binomial) with exactly the same code as logit before using my code.
2011 Jun 13
0
How to formulate an (effect-modifying) interaction with matching variable in a conditional logistic regression?
...0 = BWT>2500g # age Age of Mother Years # lwt Weight of Mother at Pounds # Last Menstrual Period # race Race 1 = White # 2 = Black # 3 = Other # smoke Smoking Status 0 = No,1 = Yes # During Pregnancy # ptd History of Premature Labor 0 = None,1 = Yes # ht History of Hypertension 0 = No, 1 = Yes # ui Presence of Uterine 0 = No, 1 = Yes # Irritability pair<-rep(1:56, each=2) low<-rep(c(0,1), 56...
2008 Mar 28
2
Comparing proportions between groups
Hello there, I have two groups (men and women) and I know per group how many of them smoke or don't smoke (women 40 of 200; men 100 of 300). I would like to know how I can compare in R if men and women differ significantly in their smoking. However, because there are more men in the sample than women I cannot just compare the number of smokers and non-smokers in both groups, right?! (I would not really know how to do that either, to be honest :(). I have done this: gensmoke <- table(data$gender,data$smoke) #save table old...
2008 Aug 19
1
R vs Stata on generalized linear mixed models: glmer and xtmelogit
Hello, I have compared the potentials of R and Stata about GLMM, analysing the dataset 'ohio' in the package 'faraway' (the same dataset is analysed with GEE in the book 'Extending the linear model with R' by Julian Faraway). Basically, I've tried the 2 commands 'glmmPQL' and 'glmer' of R and the command 'xtmelogit' of Stata. If I'm not
2015 Sep 21
2
CentOS 7 Xen 4.4.3 packages available for in virt7-xen-44-testing
On Mon, Sep 21, 2015 at 11:10:06AM +0300, Pasi K?rkk?inen wrote: > Hi, > > On Wed, Sep 09, 2015 at 02:24:15PM +0100, George Dunlap wrote: > > Updated Xen 4.4.3 packages have passed my local smoke tests, and are > > now available. > > > > Full C7 release is still waiting on suitable testing for the libvirt packages. > > > > If someone wants to step up
2009 Jun 30
0
quwstion about elrm package for Exact Logistic Regrassion
Hi I am a new R user and I did try elrm package to do Exact Logistic regression for my dataset ( about 170 patients with a binary dependent variable against age and 6 other binary covariates), I did not have n column because it was about individuals, I made n column with 1 for each individual, here is my code and result: > > #### Libraries needed #### > >
2007 Apr 25
3
aggregate similar to SPSS
Hi, Does anyone know if: with R can you take a set of numbers and aggregate them like you can in SPSS? For example, could you calculate the percentage of people who smoke based on a dataset like the following: smoke = 1 non-smoke = 2 variable 1 1 1 2 2 1 1 1 2 2 2 2 2 2 When aggregated, SPSS can tell you what percentage of persons are smokers based on the frequency of 1's and 2's. Can
2007 Jun 13
0
printer woes
greetings. my problem is i can't print from windows 2000 or windows xp. my solaris box can print fine. on smoking (solaris 10, samba host) smbclient '//smoking/hp2605' -U jgs%password put /etc/motd i get the printout. on windows, this works: echo hiya^L > \\smoking\hp2605 on windows, from M$ word, the print SEEMS to go into the queue. i click on the printer, get the box that supposedly shows me...
2008 Mar 10
3
A stats question -- about survival analysis and censoring
Dear UseRs, Suppose I have data regarding smoking habits of a prospective cohort and wish to determine the risk ratio of colorectal cancer in the smokers compared to the non-smokers. What do I do at the end of the study with people who die of heart disease? Can I just censor them exactly the same as people who become uncontactable or who die in a...
2006 Jul 21
0
[Fwd: Re: Parameterization puzzle]
Bother! This cold has made me accident-prone. I meant to hit Reply-all. Clarification below. -------- Original Message -------- Subject: Re: [R] Parameterization puzzle Date: Fri, 21 Jul 2006 19:10:03 +1200 From: Murray Jorgensen <maj at waikato.ac.nz> To: Prof Brian Ripley <ripley at stats.ox.ac.uk> References: <44C063E5.3020703 at waikato.ac.nz>
2011 Mar 08
1
Sorting
I apologize in advance if this is posted all ready I have not been able to find any information about it. I have this data frame and I want to sort smoking by retlevel. Age Gender BMI Calories Fat Fiber Alc retlevel Smoking 1 64 Female 18.87834 1828.0 63.4 14.7 0.0 Normal Non-Smoker 2 25 Female 20.64102 1517.4 59.1 5.9 0.0 Normal Smoker 3 50 Female 20.40345 1902.9 72.9 35.4 7.3 Norm...
2006 Jul 21
1
Parameterization puzzle
Consider the following example (based on an example in Pat Altham's GLM notes) pyears <- scan() 18793 52407 10673 43248 5710 28612 2585 12663 1462 5317 deaths <- scan() 2 32 12 104 28 206 28 186 31 102 Smoke <- gl(2,1,10,labels=c("No","Yes")) Age <- gl(5,2,10,labels=c("35-44","45-54","55-64","65-74","75-84"),
2018 Apr 02
2
Proposal to make Design Spec and Document for a feature mandatory.
Hi all, A better documentation about the feature, and also information about how to use the features are one of the major ask of the community when they want to use glusterfs, or want to contribute by helping get the features, bug fixes for features, etc. Finally, we have taken some baby steps to get that ask of having better design and documentation resolved. We had discussed this in our
2011 Jan 06
1
need help for chi-squared test
I've got a dataset which looks like this in the beginning: cbr dust smoking expo 1 0 0.20 1 5 2 0 0.25 1 4 3 0 0.25 1 8 4 0 0.25 1 4 5 0 0.25 1 4 (till no. 1240, anyway, a huge set) I have to analyse cbr and smoking, I know it works with chisq.test() for the whole set, but I only need cbr and smo...
2003 May 19
2
To update() or not to update()?
Hi, Suppose I have: # Fit a base model d1.ph <- coxph(Surv(start, stop, event)~ ejec + diavol + score + smoking + beta + surg.done, data = data.frame(foo)) summary(update(d1.ph, . ~ . + td1)) summary(update(d1.ph, . ~ . + td2)) As I have many columns in my data frame, foo, called td's. e.g. td1, td2, td3, .... And I'd like to add one column each time. What is...