[R] Dose-response relationship using metafor?

2013-09-18 Thread Ross, Stephanie
I am currently working a meta-analysis that is exploring the effect of a drug 
on plasma lipid levels. However, I am primarily interested in assessing the 
overall dose-response relationship with mean change in lipids. I have used the 
metafor package to conduct a meta-regression but is there a way to determine if 
this relationship is linear? Something like a p-value for trend?
Thanks so much!


PHRI DISCLAIMER This information is directed in confidence solely to the person 
named above and may not otherwise be distributed, copied or disclosed. 
Therefore, this information should be considered strictly confidential. If you 
have received this email in error, please notify the sender immediately via a 
return email for further direction. Thank you for your assistance.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
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.


[R] Contingency table and mean(sd)

2012-05-19 Thread Ross, Stephanie
Hi All,

I have a question regarding contingency tables. I would like to calculate the 
mean and standard deviation of a continuous variable from my own dataset based 
on the percentages of a contingency table I obtained from a scientific article.

dataset- data.frame(cbind(case=rep(0:1,5), x1=c(1:10), x2=c(0:9)))
contingency_table- matrix(c(100, 75, 65, 85, 90, 87), nrow=3, ncol=3)

In other words I would like to calculate the mean and standard deviation of 
X1 based on the frequencies of the contingency table.

Thank you so much for all your help!


PHRI DISCLAIMER
This information is directed in confidence solely to the person named above and 
may not otherwise be distributed, copied or disclosed. Therefore, this 
information should be considered strictly confidential. If you have received 
this email in error, please notify the sender immediately via a return email 
for further direction. Thank you for your assistance.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
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.


[R] Forest Plot Text

2011-01-31 Thread Ross, Stephanie
Hello All,

I am having a problem creating text for my forest plot using the forestplot 
function in the rmeta package.  I think my problem is that I have too many 
columns of text, so I was wondering if there is a way to shrink my text or to 
change the default setting.

Here is an example:

forest_mean- as.numeric(c(NA, 1.00, 2.00, NA, 0.75, 1.10))
forest_lower- as.numeric(c(NA, 0.85, 1.35, NA, 0.70, 0.80))
forest_upper- as.numeric(c(NA, 1.15, 2.65, NA, 0.80, 1.40))
forest_pval- as.numeric(c(NA, 0.004, 0.100,NA, 0.360, 0.080))
forest_int- as.numeric(c(NA,NA,NA,0.00, NA, NA))
forest_number_clop- c(NA, 10.0% (10/100), 5% (5/100), NA, 9% (9/100), 
4% (4/100))
forest_number_plac- c(NA, 11% (11/100),  30% (30/100), NA, 15% (15/100), 
20% (20/100))
forest_text- c(NA, 1.00 (0.85-1.15), 2.00 (1.35-2.65), NA, 0.75 
(0.70-0.80), 1.10 (0.80-1.40))
Label_text_forest- c(Outcome A, Exposure 1, Exposure 2, Outcome B, 
Exposure 1, Exposure 2)

forest_text_table- cbind(Label_text_forest, forest_number_clop, 
forest_number_plac, forest_text, forest_pval, forest_int)

plot1- forestplot(forest_text_table, forest_mean, forest_lower, forest_upper, 
digits=2,graphwidth = unit(2,inches), zero=1, clip = c(0.1, 2.50), 
title(main=Forest Plot))

I get the error message: Error in unit(rep(1, sum(widthcolumn)), grobwidth, 
labels[[1]][widthcolumn]) :   'x' and 'units' must have length  0.

Does anyone know how to fix this?

S



PHRI DISCLAIMER
This information is directed in confidence solely to the person named above and 
may not otherwise be distributed, copied or disclosed. Therefore, this 
information should be considered strictly confidential. If you have received 
this email in error, please notify the sender immediately via a return email 
for further direction. Thank you for your assistance.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
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.


[R] Combing forest plots

2011-01-26 Thread Ross, Stephanie
Hi All,

I am trying to combine two forest plots on the same page using the forestplot 
function in the rmeta package. Once I use the par() function to combine my 
plots on the same page, I find that my two plots are overlaying each other. 
Does anyone have any suggestions on how to  fix this?

Thanks!



PHRI DISCLAIMER
This information is directed in confidence solely to the person named above and 
may not otherwise be distributed, copied or disclosed. Therefore, this 
information should be considered strictly confidential. If you have received 
this email in error, please notify the sender immediately via a return email 
for further direction. Thank you for your assistance.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
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.


[R] Three-way interaction plot

2011-01-19 Thread Ross, Stephanie
Hello All,

I need to create a  three-way interaction plot and I was wondering if anyone 
could recommend a package I could use to design these plots.

Thank you!


PHRI DISCLAIMER
This information is directed in confidence solely to the person named above and 
may not otherwise be distributed, copied or disclosed. Therefore, this 
information should be considered strictly confidential. If you have received 
this email in error, please notify the sender immediately via a return email 
for further direction. Thank you for your assistance.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
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.


[R] Reshape Columns

2010-12-09 Thread Ross, Stephanie
Hello,

I have a general formatting question. I have two columns of data:

ColA - c(m, m, m, m)
ColB- c(d,d,d,d)

And I would like to reorder them into a new column that looks like this:

ColC- c(m,d,m,d,m,d,m,d)

Thank you!

Stephanie


PHRI DISCLAIMER
This information is directed in confidence solely to the person named above and 
may not otherwise be distributed, copied or disclosed. Therefore, this 
information should be considered strictly confidential. If you have received 
this email in error, please notify the sender immediately via a return email 
for further direction. Thank you for your assistance.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
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.