命令求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

function (unbalanced = FALSE)
{
    v.rates <- c(0.19, 0.13, 0.5, 0.19, 0.7, 0.66, 0.98, 0.72,
        0.71, 0.57, 0.78, 1.32, 1.47, 1.64, 1.32, 2.67, 3.16,
        2.53, 3.38, 3.93, 4.84, 5.6, 4.93, 6.05, 6.83, 6.6, 8.5,
        7.65, 10.59, 10.42, 10.36, 12, 12.68, 14.34, 17.95, 14.76,
        16.37, 18, 17.6, 23.91, 20.53, 22.6, 24.9, 24.33, 32.7,
        26.24, 27.7, 30.47, 36.94, 38.47, 33.47, 33.61, 36.77,
        43.69, 45.2)
    v.cases <- c(3, 2, 7, 3, 10, 11, 16, 11, 10, 7, 11, 22, 24,
        25, 15, 36, 44, 42, 53, 48, 77, 74, 68, 99, 88, 106,
        131, 99, 142, 134, 157, 184, 189, 180, 177, 193, 232,
        262, 249, 239, 219, 267, 323, 325, 343, 223, 250, 308,
        412, 358, 198, 214, 253, 338, 312)
    row.names <- c("25-29", "30-34", "35-39", "40-44", "45-49",
        "50-54", "55-59", "60-64", "65-69", "70-74", "75-79")
    col.names <- c("1955-1959", "1960-1964", "1965-1969", "1970-1974",
        "1975-1978")
    rates <- matrix(data = v.rates, nrow = 11, ncol = 5, byrow = TRUE,
        dimnames = list(row.names, col.names))
    cases <- matrix(data = v.cases, nrow = 11, ncol = 5, byrow = TRUE,
        dimnames = list(row.names, col.names))
    index.col <- seq(1, 4)
    if (unbalanced)
        index.col <- seq(1:5)
    return(list(rates = rates[, index.col], response = cases[,
        index.col], dose = cases[, index.col]/rates[, index.col],
        data.format = "AP", age1 = 25, per1 = 1955, coh1 = NULL,
        unit = 5, per.zero = NULL, per.max = NULL, time.adjust = 0))
}
<environment: namespace:apc>
>

谁能帮看看这是怎么建立的数据库，seq(1,4这类的说的什么

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：unbalanced unbalance function Balanced balance function

1. > bal <- function(unbalanced = FALSE){
   
2. +     v.rates <- c(0.19, 0.13, 0.5, 0.19, 0.7, 0.66, 0.98, 0.72,
   
3. +                  0.71, 0.57, 0.78, 1.32, 1.47, 1.64, 1.32, 2.67, 3.16,
   
4. +                  2.53, 3.38, 3.93, 4.84, 5.6, 4.93, 6.05, 6.83, 6.6, 8.5,
   
5. +                  7.65, 10.59, 10.42, 10.36, 12, 12.68, 14.34, 17.95, 14.76,
   
6. +                  16.37, 18, 17.6, 23.91, 20.53, 22.6, 24.9, 24.33, 32.7,
   
7. +                  26.24, 27.7, 30.47, 36.94, 38.47, 33.47, 33.61, 36.77,
   
8. +                  43.69, 45.2)
   
9. +     v.cases <- c(3, 2, 7, 3, 10, 11, 16, 11, 10, 7, 11, 22, 24,
   
10. +                  25, 15, 36, 44, 42, 53, 48, 77, 74, 68, 99, 88, 106,
    
11. +                  131, 99, 142, 134, 157, 184, 189, 180, 177, 193, 232,
    
12. +                  262, 249, 239, 219, 267, 323, 325, 343, 223, 250, 308,
    
13. +                  412, 358, 198, 214, 253, 338, 312)
    
14. +     row.names <- c("25-29", "30-34", "35-39", "40-44", "45-49",
    
15. +                    "50-54", "55-59", "60-64", "65-69", "70-74", "75-79")
    
16. +     col.names <- c("1955-1959", "1960-1964", "1965-1969", "1970-1974",
    
17. +                    "1975-1978")
    
18. +     rates <- matrix(data = v.rates, nrow = 11, ncol = 5, byrow = TRUE,
    
19. +                     dimnames = list(row.names, col.names))
    
20. +     cases <- matrix(data = v.cases, nrow = 11, ncol = 5, byrow = TRUE,
    
21. +                     dimnames = list(row.names, col.names))
    
22. +     index.col <- if(unbalanced){seq(1:5)}else{seq(1:4)}
    
23. +     return(list(rates = rates[, index.col], response = cases[, index.col],
    
24. +                 dose = cases[, index.col]/rates[, index.col],
    
25. +                 data.format = "AP", age1 = 25, per1 = 1955, coh1 = NULL,
    
26. +                 unit = 5, per.zero = NULL, per.max = NULL, time.adjust = 0)) }
    
27. > bal()
    
28. $rates
    
29.       1955-1959 1960-1964 1965-1969 1970-1974
    
30. 25-29      0.19      0.13      0.50      0.19
    
31. 30-34      0.66      0.98      0.72      0.71
    
32. 35-39      0.78      1.32      1.47      1.64
    
33. 40-44      2.67      3.16      2.53      3.38
    
34. 45-49      4.84      5.60      4.93      6.05
    
35. 50-54      6.60      8.50      7.65     10.59
    
36. 55-59     10.36     12.00     12.68     14.34
    
37. 60-64     14.76     16.37     18.00     17.60
    
38. 65-69     20.53     22.60     24.90     24.33
    
39. 70-74     26.24     27.70     30.47     36.94
    
40. 75-79     33.47     33.61     36.77     43.69
    
41. 
    
42. $response
    
43.       1955-1959 1960-1964 1965-1969 1970-1974
    
44. 25-29         3         2         7         3
    
45. 30-34        11        16        11        10
    
46. 35-39        11        22        24        25
    
47. 40-44        36        44        42        53
    
48. 45-49        77        74        68        99
    
49. 50-54       106       131        99       142
    
50. 55-59       157       184       189       180
    
51. 60-64       193       232       262       249
    
52. 65-69       219       267       323       325
    
53. 70-74       223       250       308       412
    
54. 75-79       198       214       253       338
    
55. 
    
56. $dose
    
57.       1955-1959 1960-1964 1965-1969 1970-1974
    
58. 25-29 15.789474 15.384615 14.000000 15.789474
    
59. 30-34 16.666667 16.326531 15.277778 14.084507
    
60. 35-39 14.102564 16.666667 16.326531 15.243902
    
61. 40-44 13.483146 13.924051 16.600791 15.680473
    
62. 45-49 15.909091 13.214286 13.793103 16.363636
    
63. 50-54 16.060606 15.411765 12.941176 13.408876
    
64. 55-59 15.154440 15.333333 14.905363 12.552301
    
65. 60-64 13.075881 14.172266 14.555556 14.147727
    
66. 65-69 10.667316 11.814159 12.971888 13.357994
    
67. 70-74  8.498476  9.025271 10.108303 11.153221
    
68. 75-79  5.915745  6.367153  6.880609  7.736324
    
69. 
    
70. $data.format
    
71. [1] "AP"
    
72. 
    
73. $age1
    
74. [1] 25
    
75. 
    
76. $per1
    
77. [1] 1955
    
78. 
    
79. $coh1
    
80. NULL
    
81. 
    
82. $unit
    
83. [1] 5
    
84. 
    
85. $per.zero
    
86. NULL
    
87. 
    
88. $per.max
    
89. NULL
    
90. 
    
91. $time.adjust
    
92. [1] 0
    
93. 
    
94. > bal(TRUE)
    
95. $rates
    
96.       1955-1959 1960-1964 1965-1969 1970-1974 1975-1978
    
97. 25-29      0.19      0.13      0.50      0.19      0.70
    
98. 30-34      0.66      0.98      0.72      0.71      0.57
    
99. 35-39      0.78      1.32      1.47      1.64      1.32
    
100. 40-44      2.67      3.16      2.53      3.38      3.93
     
101. 45-49      4.84      5.60      4.93      6.05      6.83
     
102. 50-54      6.60      8.50      7.65     10.59     10.42
     
103. 55-59     10.36     12.00     12.68     14.34     17.95
     
104. 60-64     14.76     16.37     18.00     17.60     23.91
     
105. 65-69     20.53     22.60     24.90     24.33     32.70
     
106. 70-74     26.24     27.70     30.47     36.94     38.47
     
107. 75-79     33.47     33.61     36.77     43.69     45.20
     
108. 
     
109. $response
     
110.       1955-1959 1960-1964 1965-1969 1970-1974 1975-1978
     
111. 25-29         3         2         7         3        10
     
112. 30-34        11        16        11        10         7
     
113. 35-39        11        22        24        25        15
     
114. 40-44        36        44        42        53        48
     
115. 45-49        77        74        68        99        88
     
116. 50-54       106       131        99       142       134
     
117. 55-59       157       184       189       180       177
     
118. 60-64       193       232       262       249       239
     
119. 65-69       219       267       323       325       343
     
120. 70-74       223       250       308       412       358
     
121. 75-79       198       214       253       338       312
     
122. 
     
123. $dose
     
124.       1955-1959 1960-1964 1965-1969 1970-1974 1975-1978
     
125. 25-29 15.789474 15.384615 14.000000 15.789474 14.285714
     
126. 30-34 16.666667 16.326531 15.277778 14.084507 12.280702
     
127. 35-39 14.102564 16.666667 16.326531 15.243902 11.363636
     
128. 40-44 13.483146 13.924051 16.600791 15.680473 12.213740
     
129. 45-49 15.909091 13.214286 13.793103 16.363636 12.884334
     
130. 50-54 16.060606 15.411765 12.941176 13.408876 12.859885
     
131. 55-59 15.154440 15.333333 14.905363 12.552301  9.860724
     
132. 60-64 13.075881 14.172266 14.555556 14.147727  9.995818
     
133. 65-69 10.667316 11.814159 12.971888 13.357994 10.489297
     
134. 70-74  8.498476  9.025271 10.108303 11.153221  9.305953
     
135. 75-79  5.915745  6.367153  6.880609  7.736324  6.902655
     
136. 
     
137. $data.format
     
138. [1] "AP"
     
139. 
     
140. $age1
     
141. [1] 25
     
142. 
     
143. $per1
     
144. [1] 1955
     
145. 
     
146. $coh1
     
147. NULL
     
148. 
     
149. $unit
     
150. [1] 5
     
151. 
     
152. $per.zero
     
153. NULL
     
154. 
     
155. $per.max
     
156. NULL
     
157. 
     
158. $time.adjust
     
159. [1] 0

复制代码

ryoeng 发表于 2015-5-22 21:59

你好，能不能加下我的QQ，我现在非常需要帮助，怕这是说不明白，这个例子是5行，11列的。但是我的数据时5行，8列的，所以我不知道下面的语句应该怎么更改，能帮帮忙不

ryoeng 发表于 2015-5-22 21:59

如果我的行和列有变化，index.col <- if(unbalanced){seq(1:5)}else{seq(1:4)}，这句话是不也发生改变呢？

咱没有使用QQ，无论几行几列都是一样的操作法，就像1+1=2, 数目加大的 10+20 就一定会是 30 而不可能是20，10+20≠20。微博是 http://www.weibo.com/5237211589

ryoeng 发表于 2015-5-25 18:14
咱没有使用QQ，无论几行几列都是一样的操作法，就像1+1=2, 数目加大的 10+20 就一定会是 30 而不可能是20， ...

if (unbalanced)
        index.col <- seq(1:5)
    return(list(rates = rates[, index.col], response = cases[,
        index.col],dose = cases[, index.col]/rates[, index.col],
        data.format = "AP", age1 = 25, per1 = 1955, coh1 = NULL,
        unit = 5, per.zero = NULL, per.max = NULL, time.adjust = 0))

你能不能告诉我下这个什么seq(1,4),然后又（1,5)啥意思。毕竟我的数据行和列和例子是同的，数据和这个seq(1,5)肯定有关的呀

ryoeng 发表于 2015-5-22 21:59

而且我想知道为什么，bal的时候就出来4列，bal(true)的时候就出来5列