全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
1393 2
2013-05-03
Could anybody please help me to create a longitudinal concentration dataset including:

100 subjects (ID)
51% female (randomly assigned)
every subject has 3 study days (study day 1, study day 30 and study day 60)
ID from 1 to 50 was given 50mg dose and rest was given 150mg; value of 50 and 150 are under colume of dose

75% sparse data (group = 1) which means on each study day, time points are 0 and random number between 1 and 3 hours (total 2 time points)

25% rich data (group = 0)which means on each study day, time points are 0,1,2,3,4,5,6,8,10,24 hours (total 10 time points)

Appreciate for teaching.

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-5-3 20:30:52
不知道有没有理解对你的意思

ID <- 1:100
sex <- rep("Male", 100)
sex[sample(1:100, 51)] <- "Female"
ID <- rep(ID, each = 3)
sex <- rep(sex, each = 3)
day <- rep(c(1, 30, 60), times=100)
dose <- rep(50, 300)
dose[ID>50] <- 150

data1 <- data.frame(ID, sex, day, dose)

sparseid <- sample(1:100, 75)
sparsetime1 <- cbind(sparseid, rep(0, 75))
sparsetime2 <- cbind(sparseid, sample(1:3, 75, replace=T))
sparsetime <- rbind(sparsetime1, sparsetime2)
sparsetime <- sparsetime[order(sparsetime[, 1]), ]

richid <- rep(setdiff(1:100, sparseid), each=10)
richtime <- rep(c(0,1,2,3,4,5,6,8,10,24), 25)
richtime <- cbind(richid, richtime)

times <- rbind(sparsetime, richtime)
times <- times[order(times[, 1]), ]

data2 <- data.frame(ID = rep(times[, 1], each=3), time = rep(times[, 2], each=3), day = rep(c(1, 30, 60), times=nrow(times)))

data <- merge(data1, data2, all=T, sort=T)

data <- data[order(data$ID, data$day, data$time), ]
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-4 04:13:40
Thank you very much! it is really helpful! Thanks again!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群