Displaying 6 results from an estimated 6 matches for "smok".
Did you mean:
smoke
2001 Jul 28
2
Re: [S] Labels wrong with lrm
...uch for your excellent description of the
problem and the self-contained test code. This is a
problem that I've been meaning to either document better
or solve for some time. The root of the problem is with
the builtin S-Plus terms.inner function:
> attr(terms.inner(asthma ~ pol(age,kx) + smok),'variables')
expression(age, kx, smok)
You can see that terms.inner inappropriately includes kx
as an independent variable as it does not know that
the first argument to pol is the special variable.
When a constant replaces kx, all is well.
As I am relying on the C code called by terms.i...
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
...I fit a coxph model without frailty (M1) and with frailty (M2) on a number of data samples with a varying degree of heterogeneity (I'm running R 2.3.1, running takes ~1 min).
library(survival); set.seed(10000)
lambda <- 0.01 # Exp. hazard rate
# Beta coefficients for Age,TC,HDLC,SBP,Diab,Smok
beta <- c(0.0483,0.0064,-0.0270,0.0037,0.4284,0.5234)
n <- 1000; Ngrp <- 2; # Nr patients, Nr frailty groups
# Thetas for gamma-frailty
thetaset <- c(1,2,10,100); Ntheta <- length(thetaset);
# Define the simulated population
age <-rnorm(n,48.6,11.7);tc<-rnorm(n,2...
2012 Jan 26
2
extracting from data.frames for survival analysis
Hi,
I have a data frame:
> class(B27.vec)
[1] "data.frame"
> head(B27.vec)
AGE Gend B27 AgeOn DD uveitis psoriasis IBD CD UC InI BASDAI BASFI Smok UV
1 57 1 1 19 38 2 1 1 1 1 1 5.40 8.08 NA 1
2 35 1 1 33 2 2 1 1 1 1 1 1.69 2.28 NA 1
3 49 2 1 40 9 1 1 1 1 1 1 8.30 9.40 NA 0
4 32 1 1 21 11 1 1 1 1 1 1 5.10 9...
2007 Jul 06
0
svyglm
...tors
are numbered 1,2,3 etc instead of 0,1,2,etc - is this a problem?
Q14sum is Eating >= 5 portions fruit & veg per day vs eating less
Q16-q19 are dietary behaviours categorised into good, ok, bad
summary(svyglm(q14sum~q16x+q17ax+q17bx+q17cx+q17dx+q17ex+q17fx+q18x+q19x
+heavy+binge+smoksum+q41sum+q46+exsum+q53+q55+q61+ethnicgrp+pcs_sum+mcs_
sum,design=dudleyls1design,family=binomial))
Call:
svyglm(q14sum ~ q16x + q17ax + q17bx + q17cx + q17dx + q17ex +
q17fx + q18x + q19x + heavy + binge + smoksum + q41sum +
q46 + exsum + q53 + q55 + q61 + ethnicgrp + pcs_sum + mc...
2009 Apr 24
2
Error building package: LaTeX error when creating PDF version
...Overfull \hbox (0.6288pt too wide) in paragraph at lines 3263--3267
\T1/ptm/m/n/10 mortality and sur-vival. Pub-lished on the in-ter-net [][]$\T1/p
cr/m/n/10 www . health . wa . gov . au / Publications /
[75]
Overfull \hbox (1.66795pt too wide) in paragraph at lines 3290--3290
[]\T1/ptm/m/n/10 the smok-ing sta-tus of the women in the \T1/ptm/m/it/10 orig-
i-nal \T1/ptm/m/n/10 sur-vey; a fac-tor with lev-els \T1/pcr/m/n/10 NonSmoker
[76] [77] [78]
No file GLMsData-manual.ind.
[79] (./GLMsData-manual.aux) )
(see the transcript file for additional information){/usr/share/texmf-texlive/f
onts/enc/dvi...
2012 Jan 18
4
R-Help
I am trying to create a frequency distribution and I am a bit confused.
Here are the commands I have entered:
> data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",")
> NumberOfActionsByStatus = data$STATUS
> NumberOfActionsByUser = data$ETS_LOGIN
> NumberOfBidOffer = data$BID_OFFER
> NumberOfActionsByUser.freq = table(NumberOfActionsByUser)
>