[R] trouble double looping to generate data for a meta-analysis

Federman, Douglas Douglas.Federman at utoledo.edu
Fri Jul 1 21:06:36 CEST 2016


You might look at the package 
     Wakefield
For data generation

-----Original Message-----
From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Marietta Suarez
Sent: Friday, July 01, 2016 1:28 PM
To: r-help at r-project.org
Subject: [R] trouble double looping to generate data for a meta-analysis

i'm trying to generate data for a meta analysis. 1- generate data following
a normal distribution, 2- generate data following a skewed distribution, 3-
generate data following a logistic distribution. i need to loop this
because the # of studies in each meta will be either 10 or 15. k or total
number of studies in the meta will be 5. i need to loop twice to repeat
this process 10 times. database should be 3 columns (distributions) by 65
rows x 10 reps


here's my code, not sure what's not working:
library(fGarch)

#n reps =10
rep=10

#begin function here, need to vary n and k, when k=2 n=10, when k3 n=15
fun=function(n, k){

  #prepare to store data
  data=matrix(0,nrow=10*k, ncol=3)
  db=matrix(0,nrow=650, ncol=3)

  for (j in 1:rep)
  {
    for (i in 1:k)
    {
      #generate data under normal, skewed, and logistic distributions here

      data[,1]=rnorm(n, 100, 15)
      data[,2]=rsnorm(n, 100, 15, 1)
      data[,3]=rlogis(n, 100, 15)
    }
      [j]=db
  }
}

save=fun(10,2)

Please help!!!

	[[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.



More information about the R-help mailing list