Displaying 7 results from an estimated 7 matches for "ci_l".
2018 Mar 06
4
Capturing warning within user-defined function
...# subset
grp1 <- subset(my.svy, grp == "Group1")
# set options and clear old warnings
options(warn=0)
assign("last.warning", NULL, envir = baseenv())
## proportions and CIs
p <- ((svyciprop(~grp, grp1, family=quasibinomial))[1])
# save warnings
wrn1 <- warnings(p)
ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1])
ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2])
## sample counts
n <- unwtd.count(~grp, grp1)[1]
## combine into table
overall <- data.frame(n, p, ci_l, ci_u)
colnames(overall) &...
2018 Mar 06
0
Capturing warning within user-defined function
...")
>
> # set options and clear old warnings
> options(warn=0)
> assign("last.warning", NULL, envir = baseenv())
>
> ## proportions and CIs
> p <- ((svyciprop(~grp, grp1, family=quasibinomial))[1])
>
> # save warnings
> wrn1 <- warnings(p)
>
> ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1])
> ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2])
>
> ## sample counts
> n <- unwtd.count(~grp, grp1)[1]
>
> ## combine into table
> overall <- data.frame(n, p...
2018 Mar 06
0
Capturing warning within user-defined function
...")
>
> # set options and clear old warnings
> options(warn=0)
> assign("last.warning", NULL, envir = baseenv())
>
> ## proportions and CIs
> p <- ((svyciprop(~grp, grp1, family=quasibinomial))[1])
>
> # save warnings
> wrn1 <- warnings(p)
>
> ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1])
> ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2])
>
> ## sample counts
> n <- unwtd.count(~grp, grp1)[1]
>
> ## combine into table
> overall <- data.frame(n, p...
2018 Mar 06
1
Capturing warning within user-defined function
...ings
> > options(warn=0)
> > assign("last.warning", NULL, envir = baseenv())
> >
> > ## proportions and CIs
> > p <- ((svyciprop(~grp, grp1, family=quasibinomial))[1])
> >
> > # save warnings
> > wrn1 <- warnings(p)
> >
> > ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1])
> > ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2])
> >
> > ## sample counts
> > n <- unwtd.count(~grp, grp1)[1]
> >
> > ## combine into table
>...
2009 Nov 16
1
lmomco package and confidence limits?
Hello,
I am using the lmomco package (lmom.ub and pargev) to compute the GEV
parameters (location, scale, and shape), which are used to estimate
return values. I was wondering how/if I can calculate upper and lower
confidence (CI_u, CI_l) intervals for each return frequency using the
GEV parameters to fill-in the table below?
Xi (location) = 35.396
Alpha (scale) = 1.726
Kappa (shape) = 0.397
Year, Return value, CI_u, CL_l
2, 36.0,
5, 37.3,
10, 38.0,
25, 38.5,
50, 38.8,
100, 39.0,
Thank you for your help/suggestions,
Doug...
2013 Jan 03
1
Survey package help with svystandardize
...2.5%
97.5%
I(rfsmok == "At risk") 0.187685 0.176309 0.19962
Estimates by sex are a match.
> svyby(~I(rfsmok=="At risk"), ~sex, stdes, svyciprop, na.rm=TRUE, ci=TRUE,
vartype="ci")
sex I(rfsmok == "At risk") ci_l
ci_u
Male Male 0.2277632725 0.2097154285 0.2468790907
Female Female 0.1486524831 0.1349832066 0.1634444481
-------------------------
R Code:
sample.brfss <- svydesign(id=~brfss.2011$SEQNO, strata=~brfss.2011$ststr,
weights=~brfss.2011$LLCPWT, data=brfss.2011,...
2007 Jun 27
2
Meta-Analysis of proportions
Dear colleagues,
I'm conducting a meta-analysis of studies evaluating adherence of HIV-positive drug users into AIDS treatment, therefore I'm looking for some advice and syntax suggestion for running the meta-regression using proportions, not the usual OR/RR frequently used on RCT studies.
Have already searched already several handbooks, R-manuals, mailing lists, professors, but... not