楼主: verbatim
1209 2

[问答] help to add error to the time points, thanks [推广有奖]

  • 1关注
  • 3粉丝

已卖：25份资源

博士生

5%

还不是VIP/贵宾

-

0%
威望
0 级
论坛币
513 个
通用积分
2.0505
学术水平
0 点
热心指数
2 点
信用等级
0 点
经验
3404 点
帖子
139
精华
0
在线时间
143 小时
注册时间
2006-12-1
最后登录
2020-1-30
楼主
verbatim 发表于 2013-5-24 02:06:02 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币
  • First, I created the below dataset:

    allData = structure(list(C = structure(c(1L, NA, NA, NA, NA, NA, NA, NA,
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
    NA, NA, NA, NA, NA, NA, NA, NA, NA), .Label = "C", class = "factor"),
    TRL = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L), ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), DOSE = c(50, 50,
    50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50,
    50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50,
    50), SDAY = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 30, 30, 30,
    30, 30, 30, 30, 30, 30, 30, 30, 60, 60, 60, 60, 60, 60, 60,
    60, 60, 60, 60), TIME = c(0, 0, 1, 2, 3, 4, 5, 6, 8, 10,
    24, 0, 0, 1, 2, 3, 4, 5, 6, 8, 10, 24, 0, 0, 1, 2, 3, 4,
    5, 6, 8, 10, 24), GEN = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), GROUP = c(0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L), TSFD = c(0, 0, 1, 2, 3, 4, 5, 6, 8, 10, 24, 696, 696,
    697, 698, 699, 700, 701, 702, 704, 706, 720, 1416, 1416,
    1417, 1418, 1419, 1420, 1421, 1422, 1424, 1426, 1440), OBS = 1:33,
    CL = c(33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
    33.8712839749554, 33.8712839749554), V = c(320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
    320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579
    ), KA = c(0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293,
    0.995274973731293, 0.995274973731293, 0.995274973731293),
    KEL = c(0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526,
    0.105835444986526, 0.105835444986526, 0.105835444986526),
    conc = c(0, 0, 92.647083382467, 117.587164906218, 118.435668936488,
    111.21996697943, 101.779709034282, 92.1974188524962, 74.9088584827009,
    60.6595346538829, 13.7868848267872, 14.9672371118036, 14.9672371118036,
    106.111200764071, 129.699116863709, 129.331249838855, 121.021333729664,
    110.59675064338, 100.128988471021, 81.3273303362704, 65.8535609206108,
    14.9672371118036, 14.9672371118036, 14.9672371118036, 106.111200764071,
    129.699116863709, 129.331249838855, 121.021333729664, 110.59675064338,
    100.128988471021, 81.3273303362704, 65.8535609206108, 14.9672371118036
    ), obs = c(0, 0, 93.0705373080666, 119.158931115447, 101.119131737192,
    98.3850237974117, 93.3146423825588, 94.4074655857413, 85.4798011942845,
    64.3246842238628, 12.3877936201801, 12.9611327916651, 12.9611327916651,
    94.4599756673357, 128.418170457443, 127.880610571531, 112.014925161075,
    117.557276099122, 119.781561548053, 78.4235789891384, 61.6785269162022,
    13.3827881024542, 15.3454182191101, 15.3454182191101, 91.6833753586737,
    133.544835331848, 139.86721964799, 124.664449544007, 102.408937731944,
    93.5604355646615, 83.1683032098305, 76.5718951980788, 16.53624694588
    ), CMT = c(2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2,
    2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2), MDV = c(0,
    1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), EVID = c(0, 1, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), AMT = c(NA, 50000, NA,
    NA, NA, NA, NA, NA, NA, NA, NA, NA, 50000, NA, NA, NA, NA,
    NA, NA, NA, NA, NA, NA, 50000, NA, NA, NA, NA, NA, NA, NA,
    NA, NA), DV = c(0, 0, 93.0705373080666, 119.158931115447,
    101.119131737192, 98.3850237974117, 93.3146423825588, 94.4074655857413,
    85.4798011942845, 64.3246842238628, 12.3877936201801, 12.9611327916651,
    0, 94.4599756673357, 128.418170457443, 127.880610571531,
    112.014925161075, 117.557276099122, 119.781561548053, 78.4235789891384,
    61.6785269162022, 13.3827881024542, 15.3454182191101, 0,
    91.6833753586737, 133.544835331848, 139.86721964799, 124.664449544007,
    102.408937731944, 93.5604355646615, 83.1683032098305, 76.5718951980788,
    16.53624694588), II = c(NA, NA, NA, NA, NA, NA, NA, NA, NA,
    NA, NA, NA, 24, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 24,
    NA, NA, NA, NA, NA, NA, NA, NA, NA), SS = c(NA, NA, NA, NA,
    NA, NA, NA, NA, NA, NA, NA, NA, 2, NA, NA, NA, NA, NA, NA,
    NA, NA, NA, NA, 2, NA, NA, NA, NA, NA, NA, NA, NA, NA)), .Names = c("C",
    "TRL", "ID", "DOSE", "SDAY", "TIME", "GEN", "GROUP", "TSFD",
    "OBS", "CL", "V", "KA", "KEL", "conc", "obs", "CMT", "MDV", "EVID",
    "AMT", "DV", "II", "SS"), row.names = c("1201", "1", "4", "7",
    "10", "9", "3", "6", "5", "8", "2", "1510", "15", "13", "12",
    "11", "14", "17", "20", "19", "18", "16", "2910", "29", "23",
    "26", "25", "28", "22", "21", "24", "27", "30"), class = "data.frame")

    Then to add error the sampling time points in TSFD (EVID=0)

    allData$TSFD[allData$EVID == 0 & allData$TSFD >0] = round(jitter(allData$TSFD[allData$EVID == 0 & allData$TSFD >0]),2)




then, I got the below dataset:

allData = structure(list(C = structure(c(1L, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA), .Label = "C", class = "factor"),
TRL = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L), ID = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), DOSE = c(50, 50,
50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50,
50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50,
50), SDAY = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 30, 30, 30,
30, 30, 30, 30, 30, 30, 30, 30, 60, 60, 60, 60, 60, 60, 60,
60, 60, 60, 60), TIME = c(0, 0, 1, 2, 3, 4, 5, 6, 8, 10,
24, 0, 0, 1, 2, 3, 4, 5, 6, 8, 10, 24, 0, 0, 1, 2, 3, 4,
5, 6, 8, 10, 24), GEN = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), GROUP = c(0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L), TSFD = c(0, 0, 1.07, 2.19, 3, 4.15, 4.82, 6.06, 7.96,
10.02, 23.81, 696.1, 696, 696.9, 697.92, 698.87, 699.94,
701.11, 701.96, 704.19, 705.89, 720.06, 1416, 1416, 1417.11,
1417.81, 1419, 1419.84, 1421.06, 1422.15, 1424.11, 1426.04,
1439.83), OBS = 1:33, CL = c(33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554, 33.8712839749554,
33.8712839749554, 33.8712839749554, 33.8712839749554), V = c(320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579,
320.037242525579, 320.037242525579, 320.037242525579, 320.037242525579
), KA = c(0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293,
0.995274973731293, 0.995274973731293, 0.995274973731293),
KEL = c(0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526,
0.105835444986526, 0.105835444986526, 0.105835444986526),
conc = c(0, 0, 92.647083382467, 117.587164906218, 118.435668936488,
111.21996697943, 101.779709034282, 92.1974188524962, 74.9088584827009,
60.6595346538829, 13.7868848267872, 14.9672371118036, 14.9672371118036,
106.111200764071, 129.699116863709, 129.331249838855, 121.021333729664,
110.59675064338, 100.128988471021, 81.3273303362704, 65.8535609206108,
14.9672371118036, 14.9672371118036, 14.9672371118036, 106.111200764071,
129.699116863709, 129.331249838855, 121.021333729664, 110.59675064338,
100.128988471021, 81.3273303362704, 65.8535609206108, 14.9672371118036
), obs = c(0, 0, 93.0705373080666, 119.158931115447, 101.119131737192,
98.3850237974117, 93.3146423825588, 94.4074655857413, 85.4798011942845,
64.3246842238628, 12.3877936201801, 12.9611327916651, 12.9611327916651,
94.4599756673357, 128.418170457443, 127.880610571531, 112.014925161075,
117.557276099122, 119.781561548053, 78.4235789891384, 61.6785269162022,
13.3827881024542, 15.3454182191101, 15.3454182191101, 91.6833753586737,
133.544835331848, 139.86721964799, 124.664449544007, 102.408937731944,
93.5604355646615, 83.1683032098305, 76.5718951980788, 16.53624694588
), CMT = c(2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2), MDV = c(0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), EVID = c(0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), AMT = c(NA, 50000, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, 50000, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, 50000, NA, NA, NA, NA, NA, NA, NA,
NA, NA), DV = c(0, 0, 93.0705373080666, 119.158931115447,
101.119131737192, 98.3850237974117, 93.3146423825588, 94.4074655857413,
85.4798011942845, 64.3246842238628, 12.3877936201801, 12.9611327916651,
0, 94.4599756673357, 128.418170457443, 127.880610571531,
112.014925161075, 117.557276099122, 119.781561548053, 78.4235789891384,
61.6785269162022, 13.3827881024542, 15.3454182191101, 0,
91.6833753586737, 133.544835331848, 139.86721964799, 124.664449544007,
102.408937731944, 93.5604355646615, 83.1683032098305, 76.5718951980788,
16.53624694588), II = c(NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, 24, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 24,
NA, NA, NA, NA, NA, NA, NA, NA, NA), SS = c(NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, 2, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, 2, NA, NA, NA, NA, NA, NA, NA, NA, NA)), .Names = c("C",
"TRL", "ID", "DOSE", "SDAY", "TIME", "GEN", "GROUP", "TSFD",
"OBS", "CL", "V", "KA", "KEL", "conc", "obs", "CMT", "MDV", "EVID",
"AMT", "DV", "II", "SS"), row.names = c("1201", "1", "4", "7",
"10", "9", "3", "6", "5", "8", "2", "1510", "15", "13", "12",
"11", "14", "17", "20", "19", "18", "16", "2910", "29", "23",
"26", "25", "28", "22", "21", "24", "27", "30"), class = "data.frame")

However, the sequence of TSFD is not correct, if you look at ID 1, for example, first two records at SDAY=30, first record is pre-dose record, second record is doing record (EVID=1), but the sequence of time is wrong (equal is acceptable). how can I add error by jitter() to have correct time sequence?




二维码

扫码加我 拉你入群

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

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

分享0
关键词：thanks Points Error Point Thank structure created points factor error

相关帖子

沙发
qoiqpwqr 发表于 2013-5-24 02:56:30
你是想TSFD是从小到大的顺序？
藤椅
verbatim 发表于 2013-5-24 03:30:13
yes. But pre-dose record must be ahead of doing record (EVID=1).
返回列表
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2025-12-21 08:17