Hello Im analyzing a survey and i need to obtain some statistics per groups. Im able to create a table with sex and age. However, if i want to know how much income earns the population by sex and age, i can't. Im loading the dataset as describe the line below NN <- svydesign(ids = ~1, data = encuesta, weights = fact) Some simple table i can create table(svytable(~age+sex,design=NN)) But im not able to handle the same tabulate referencing a income variable, in this case, wage. Can you help me? Thanks for your replies and time.
svytable: How do i create a table informing a third variable?
5 messages · Anthony Damico, Juan Ceccarelli Arias, David Winsemius
# mean svymean( ~ income_variable , NN ) svyby( ~ income_variable , ~ age + sex , NN , svymean ) # median svyquantile( ~ income_variable , NN ) svyby( ~ income_variable , ~ age + sex , NN , svyquantile , 0.5 ) On Fri, Sep 2, 2016 at 3:04 PM, Juan Ceccarelli Arias <jfca283 at gmail.com> wrote:
Hello
Im analyzing a survey and i need to obtain some statistics per groups.
Im able to create a table with sex and age. However, if i want to know how
much income earns the population by sex and age, i can't.
Im loading the dataset as describe the line below
NN <- svydesign(ids = ~1, data = encuesta, weights = fact)
Some simple table i can create
table(svytable(~age+sex,design=NN))
But im not able to handle the same tabulate referencing a income variable,
in this case, wage.
Can you help me?
Thanks for your replies and time.
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/ posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Thanks a lot. Your code does the trick.
One last question:
The tabulate produced is showing every cross in just one column.
I mean, it presents the region by order and sex=1, and then again the
region but by sex==2.
Can i list or present as this:
sex1 sex2
region1 323. 3434..
...
regionN 123.. 432..
and ignoring the remaining info (standar errors or se in this case)?
Again, thanks Anthony. Really.
On Fri, Sep 2, 2016 at 8:24 PM, Anthony Damico <ajdamico at gmail.com> wrote:
# mean svymean( ~ income_variable , NN ) svyby( ~ income_variable , ~ age + sex , NN , svymean ) # median svyquantile( ~ income_variable , NN ) svyby( ~ income_variable , ~ age + sex , NN , svyquantile , 0.5 ) On Fri, Sep 2, 2016 at 3:04 PM, Juan Ceccarelli Arias <jfca283 at gmail.com> wrote:
Hello
Im analyzing a survey and i need to obtain some statistics per groups.
Im able to create a table with sex and age. However, if i want to know how
much income earns the population by sex and age, i can't.
Im loading the dataset as describe the line below
NN <- svydesign(ids = ~1, data = encuesta, weights = fact)
Some simple table i can create
table(svytable(~age+sex,design=NN))
But im not able to handle the same tabulate referencing a income variable,
in this case, wage.
Can you help me?
Thanks for your replies and time.
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posti ng-guide.html and provide commented, minimal, self-contained, reproducible code.
On Sep 2, 2016, at 6:08 PM, Juan Ceccarelli Arias <jfca283 at gmail.com> wrote:
Thanks a lot. Your code does the trick.
One last question:
The tabulate produced is showing every cross in just one column.
I mean, it presents the region by order and sex=1, and then again the
region but by sex==2.
Can i list or present as this:
sex1 sex2
region1 323. 3434..
...
regionN 123.. 432..
and ignoring the remaining info (standar errors or se in this case)?
Again, thanks Anthony. Really.
(Anthony's probably asleep.) This doesn't ignore the se's but that could be easily done by omitting that column from the data argument:
Your help was everything i needed it. Please, declare this topic as solved. And thanks again. On Sat, Sep 3, 2016 at 10:06 PM, David Winsemius <dwinsemius at comcast.net> wrote:
On Sep 2, 2016, at 6:08 PM, Juan Ceccarelli Arias <jfca283 at gmail.com>
wrote:
Thanks a lot. Your code does the trick.
One last question:
The tabulate produced is showing every cross in just one column.
I mean, it presents the region by order and sex=1, and then again the
region but by sex==2.
Can i list or present as this:
sex1 sex2
region1 323. 3434..
...
regionN 123.. 432..
and ignoring the remaining info (standar errors or se in this case)?
Again, thanks Anthony. Really.
(Anthony's probably asleep.) This doesn't ignore the se's but that could be easily done by omitting that column from the data argument: From the examples on the help page for svymean:
svyby( ~ mobility , ~ stype + comp.imp , dclus1 , svymean )
stype comp.imp mobility se E.No E No 19.71875 1.347583 H.No H No 13.14286 0.740017 M.No M No 14.81818 2.960618 E.Yes E Yes 17.28571 1.536158 H.Yes H Yes 35.14286 16.570001 M.Yes M Yes 13.71429 2.628573 apimeans1 <- svyby( ~ mobility , ~ stype + comp.imp , dclus1 , svymean )
reshape(apimeans1, idvar='stype', direction="wide", timevar="comp.imp")
stype mobility.No se.No mobility.Yes se.Yes E.No E 19.71875 1.347583 17.28571 1.536158 H.No H 13.14286 0.740017 35.14286 16.570001 M.No M 14.81818 2.960618 13.71429 2.628573 -- David.
On Fri, Sep 2, 2016 at 8:24 PM, Anthony Damico <ajdamico at gmail.com>
wrote:
# mean svymean( ~ income_variable , NN ) svyby( ~ income_variable , ~ age + sex , NN , svymean ) # median svyquantile( ~ income_variable , NN ) svyby( ~ income_variable , ~ age + sex , NN , svyquantile , 0.5 ) On Fri, Sep 2, 2016 at 3:04 PM, Juan Ceccarelli Arias <
jfca283 at gmail.com>
wrote:
Hello Im analyzing a survey and i need to obtain some statistics per groups. Im able to create a table with sex and age. However, if i want to know
how
much income earns the population by sex and age, i can't. Im loading the dataset as describe the line below NN <- svydesign(ids = ~1, data = encuesta, weights = fact) Some simple table i can create table(svytable(~age+sex,design=NN)) But im not able to handle the same tabulate referencing a income
variable,
in this case, wage.
Can you help me?
Thanks for your replies and time.
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posti ng-guide.html and provide commented, minimal, self-contained, reproducible code.
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/
posting-guide.html
and provide commented, minimal, self-contained, reproducible code.