#7 BONUS: Read file from github
# Link: https://github.com/HildaRamirez/HW2/blob/master/pbc.csv
# HMTL link didn't work, had to use raw file
theURL <- "https://raw.githubusercontent.com/HildaRamirez/HW2/master/pbc.csv"
PBC <- read.table(file=theURL, header = TRUE, sep = ",")
head(PBC)
## X id time status trt age sex ascites hepato spiders edema bili chol
## 1 1 1 400 2 1 58.76523 f 1 1 1 1.0 14.5 261
## 2 2 2 4500 0 1 56.44627 f 0 1 1 0.0 1.1 302
## 3 3 3 1012 2 1 70.07255 m 0 0 0 0.5 1.4 176
## 4 4 4 1925 2 1 54.74059 f 0 1 1 0.5 1.8 244
## 5 5 5 1504 1 2 38.10541 f 0 1 1 0.0 3.4 279
## 6 6 6 2503 2 2 66.25873 f 0 1 0 0.0 0.8 248
## albumin copper alk.phos ast trig platelet protime stage
## 1 2.60 156 1718.0 137.95 172 190 12.2 4
## 2 4.14 54 7394.8 113.52 88 221 10.6 3
## 3 3.48 210 516.0 96.10 55 151 12.0 4
## 4 2.54 64 6121.8 60.63 92 183 10.3 4
## 5 3.53 143 671.0 113.15 72 136 10.9 3
## 6 3.98 50 944.0 93.00 63 NA 11.0 3
#1 Use summary function to gain overview.
summary(PBC)
## X id time status
## Min. : 1.0 Min. : 1.0 Min. : 41 Min. :0.0000
## 1st Qu.:105.2 1st Qu.:105.2 1st Qu.:1093 1st Qu.:0.0000
## Median :209.5 Median :209.5 Median :1730 Median :0.0000
## Mean :209.5 Mean :209.5 Mean :1918 Mean :0.8301
## 3rd Qu.:313.8 3rd Qu.:313.8 3rd Qu.:2614 3rd Qu.:2.0000
## Max. :418.0 Max. :418.0 Max. :4795 Max. :2.0000
##
## trt age sex ascites
## Min. :1.000 Min. :26.28 Length:418 Min. :0.00000
## 1st Qu.:1.000 1st Qu.:42.83 Class :character 1st Qu.:0.00000
## Median :1.000 Median :51.00 Mode :character Median :0.00000
## Mean :1.494 Mean :50.74 Mean :0.07692
## 3rd Qu.:2.000 3rd Qu.:58.24 3rd Qu.:0.00000
## Max. :2.000 Max. :78.44 Max. :1.00000
## NA's :106 NA's :106
## hepato spiders edema bili
## Min. :0.0000 Min. :0.0000 Min. :0.0000 Min. : 0.300
## 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.: 0.800
## Median :1.0000 Median :0.0000 Median :0.0000 Median : 1.400
## Mean :0.5128 Mean :0.2885 Mean :0.1005 Mean : 3.221
## 3rd Qu.:1.0000 3rd Qu.:1.0000 3rd Qu.:0.0000 3rd Qu.: 3.400
## Max. :1.0000 Max. :1.0000 Max. :1.0000 Max. :28.000
## NA's :106 NA's :106
## chol albumin copper alk.phos
## Min. : 120.0 Min. :1.960 Min. : 4.00 Min. : 289.0
## 1st Qu.: 249.5 1st Qu.:3.243 1st Qu.: 41.25 1st Qu.: 871.5
## Median : 309.5 Median :3.530 Median : 73.00 Median : 1259.0
## Mean : 369.5 Mean :3.497 Mean : 97.65 Mean : 1982.7
## 3rd Qu.: 400.0 3rd Qu.:3.770 3rd Qu.:123.00 3rd Qu.: 1980.0
## Max. :1775.0 Max. :4.640 Max. :588.00 Max. :13862.4
## NA's :134 NA's :108 NA's :106
## ast trig platelet protime
## Min. : 26.35 Min. : 33.00 Min. : 62.0 Min. : 9.00
## 1st Qu.: 80.60 1st Qu.: 84.25 1st Qu.:188.5 1st Qu.:10.00
## Median :114.70 Median :108.00 Median :251.0 Median :10.60
## Mean :122.56 Mean :124.70 Mean :257.0 Mean :10.73
## 3rd Qu.:151.90 3rd Qu.:151.00 3rd Qu.:318.0 3rd Qu.:11.10
## Max. :457.25 Max. :598.00 Max. :721.0 Max. :18.00
## NA's :106 NA's :136 NA's :11 NA's :2
## stage
## Min. :1.000
## 1st Qu.:2.000
## Median :3.000
## Mean :3.024
## 3rd Qu.:4.000
## Max. :4.000
## NA's :6
#1 Mean and median for two liver enzymes: alkaline phosphotase and aspartate transaminase (AST). Without the na.rm=TRUE (na remove) they all returned NA.
a <- mean(PBC$alk.phos, na.rm=TRUE)
b <- median(PBC$alk.phos, na.rm=TRUE)
c <- mean(PBC$ast, na.rm=TRUE)
d <- median(PBC$ast, na.rm=TRUE)
chart <- data.frame(enzyme=c("AlkPhos Mean", "AlkPhos Median", "AST Mean", "AST Median"),
allpts=c(a,b,c,d))
chart
## enzyme allpts
## 1 AlkPhos Mean 1982.6558
## 2 AlkPhos Median 1259.0000
## 3 AST Mean 122.5563
## 4 AST Median 114.7000
#2 Subset of PBC data. Use only the first 312 rows since the last ones did not actually participate in the trial and have limited data. Copy over only a selection of columns. Create subset for patients in stages 0-2. This yielded 83 results.
earlyPBC <- PBC[1:312, c(2,5:9,12:14,16,17,20,21)]
tail(earlyPBC)
## id trt age sex ascites hepato bili chol albumin alk.phos ast protime
## 307 307 2 30.57358 f 0 0 0.8 273 3.56 1282 130 10.5
## 308 308 1 61.18275 f 0 1 0.4 246 3.58 797 91 10.4
## 309 309 2 58.29979 f 0 0 0.4 260 2.75 1166 70 10.8
## 310 310 1 62.33265 f 0 0 1.7 434 3.35 1713 171 10.2
## 311 311 1 37.99863 f 0 0 2.0 247 3.16 1050 117 10.5
## 312 312 2 33.15264 f 0 0 6.4 576 3.79 2115 136 10.8
## stage
## 307 2
## 308 2
## 309 2
## 310 2
## 311 2
## 312 2
nrow(earlyPBC)
## [1] 312
earlyPBC <- subset(earlyPBC, stage < 3)
nrow(earlyPBC)
## [1] 83
#3 Rename columns
library(plyr)
rename(earlyPBC,c("trt"="treatment", "bili"="bilirubin", chol="cholesterol" ))
## id treatment age sex ascites hepato bilirubin cholesterol albumin
## 9 9 1 42.50787 f 0 0 3.2 562 3.08
## 24 24 1 44.52019 m 0 1 2.1 456 4.00
## 25 25 2 45.07324 f 0 0 0.7 298 4.10
## 29 29 2 63.87680 f 0 0 0.7 370 3.78
## 31 31 2 41.55236 f 0 1 4.7 296 3.44
## 34 34 1 52.06023 f 0 0 0.8 364 3.70
## 36 36 2 56.41068 f 0 0 0.3 172 3.39
## 43 43 1 48.87064 f 0 0 1.1 361 3.64
## 45 45 2 41.79329 f 0 0 0.6 NA 3.93
## 51 51 2 52.08761 f 0 0 0.8 276 3.60
## 52 52 1 50.54073 m 0 0 6.0 614 3.70
## 56 56 2 33.61807 f 0 1 1.1 498 3.80
## 58 58 1 44.56947 m 0 0 0.7 242 4.08
## 61 61 2 43.89870 m 0 0 0.6 216 3.94
## 65 65 1 40.20260 f 0 0 1.2 256 3.60
## 70 70 1 56.39973 f 0 0 0.6 NA 4.64
## 73 73 2 38.49418 f 0 0 0.7 132 3.60
## 89 89 1 52.44353 f 0 1 2.0 408 3.65
## 90 90 1 33.47570 m 0 0 1.6 660 4.22
## 93 93 1 36.53388 f 0 0 1.3 353 3.67
## 95 95 2 46.39014 f 1 1 17.4 NA 2.64
## 98 98 1 28.88433 f 0 0 1.0 239 3.77
## 99 99 2 48.46817 m 0 0 1.8 460 3.35
## 102 102 1 56.56947 f 0 0 0.9 248 3.97
## 104 104 1 43.01711 f 0 0 1.1 303 3.64
## 106 106 1 68.50924 f 1 1 2.1 NA 3.90
## 107 107 2 62.52156 f 0 0 0.6 212 4.03
## 108 108 1 50.35729 f 0 0 0.4 127 3.50
## 127 127 2 44.10404 f 0 0 0.5 268 4.08
## 134 134 2 42.33539 f 0 0 0.7 578 3.67
## 135 135 1 42.96783 f 0 0 0.4 263 3.57
## 136 136 1 55.96167 f 0 0 0.8 263 3.35
## 137 137 1 62.86105 f 0 0 1.1 399 3.60
## 141 141 1 47.03628 f 0 0 0.9 346 3.77
## 142 142 2 55.72621 f 0 1 1.0 364 3.48
## 147 147 1 75.01164 f 0 0 1.2 288 3.37
## 150 150 2 34.98700 f 0 0 1.0 NA 3.63
## 151 151 1 55.04175 f 0 0 0.9 460 3.03
## 153 153 1 49.60438 f 0 0 0.5 217 3.85
## 156 156 2 59.40862 f 0 1 25.5 358 3.52
## 157 157 2 48.75838 f 0 0 0.6 286 3.42
## 158 158 1 36.49281 f 0 0 3.4 450 3.37
## 161 161 2 42.74333 f 0 0 2.3 502 3.56
## 166 166 2 41.35524 f 0 1 5.7 1480 3.26
## 168 168 2 47.75359 f 0 1 0.4 257 3.80
## 173 173 1 30.27515 f 0 1 3.0 236 3.42
## 174 174 1 55.56742 f 0 0 0.5 NA 3.85
## 177 177 2 55.45243 f 0 0 0.9 258 4.01
## 182 182 1 56.94456 f 0 1 0.6 NA 3.69
## 183 183 2 40.26010 f 0 0 0.5 201 3.73
## 190 190 1 54.52704 f 0 0 2.3 NA 3.93
## 194 194 1 44.39699 f 0 0 1.5 293 4.30
## 197 197 1 44.62697 f 0 0 0.6 266 3.97
## 198 198 2 35.79740 f 0 0 0.7 286 2.90
## 202 202 1 61.63997 f 0 0 0.5 223 3.80
## 203 203 2 37.05681 f 0 1 0.5 149 4.04
## 206 206 1 61.99042 f 0 0 0.6 213 4.07
## 209 209 2 52.62423 f 0 1 0.7 252 4.01
## 210 210 2 49.76318 m 0 1 0.9 346 3.37
## 218 218 1 34.59548 f 0 0 0.5 NA 3.89
## 220 220 1 63.63039 f 0 0 3.8 426 3.22
## 225 225 1 38.77070 f 0 0 0.7 274 3.66
## 227 227 1 58.95140 m 0 0 2.3 316 3.35
## 234 234 2 34.37919 f 0 0 0.9 578 3.35
## 238 238 2 66.41205 f 0 0 0.6 280 3.35
## 251 251 1 56.97467 f 0 0 0.5 227 3.61
## 257 257 1 58.26420 f 0 0 0.5 280 4.23
## 258 258 1 51.48802 f 0 0 1.1 414 3.44
## 269 269 2 45.08419 f 0 0 3.6 374 3.50
## 272 272 1 38.39836 f 0 0 0.5 226 2.93
## 276 276 1 50.10815 f 0 0 1.0 219 3.85
## 279 279 2 56.15332 f 0 0 0.5 198 3.77
## 285 285 1 46.34908 f 0 0 0.8 253 3.48
## 288 288 2 48.93634 f 0 1 8.7 310 3.89
## 289 289 1 67.57290 m 0 1 4.0 416 3.99
## 302 302 1 37.79877 f 0 0 0.7 335 3.95
## 306 306 2 56.39151 f 0 1 0.6 239 3.45
## 307 307 2 30.57358 f 0 0 0.8 273 3.56
## 308 308 1 61.18275 f 0 1 0.4 246 3.58
## 309 309 2 58.29979 f 0 0 0.4 260 2.75
## 310 310 1 62.33265 f 0 0 1.7 434 3.35
## 311 311 1 37.99863 f 0 0 2.0 247 3.16
## 312 312 2 33.15264 f 0 0 6.4 576 3.79
## alk.phos ast protime stage
## 9 2276.0 144.15 11.0 2
## 24 5719.0 221.88 9.9 2
## 25 661.0 106.95 11.3 2
## 29 5833.0 73.53 10.6 2
## 31 9933.2 206.40 10.3 2
## 34 1840.0 170.50 10.5 2
## 36 558.0 71.30 10.6 2
## 43 5430.2 67.08 10.6 2
## 45 1826.0 71.30 10.9 2
## 51 4332.0 99.33 10.6 2
## 52 5084.4 206.40 10.6 1
## 56 13862.4 95.46 10.6 2
## 58 5890.0 56.76 10.6 1
## 61 601.0 60.45 13.0 1
## 65 724.0 141.05 10.0 1
## 70 666.0 54.25 10.6 2
## 73 423.0 49.60 11.0 1
## 89 1083.0 110.05 11.4 2
## 90 1857.0 151.90 11.0 2
## 93 2039.0 232.50 11.1 2
## 95 559.0 119.35 11.7 2
## 98 1877.0 97.65 10.2 1
## 99 1472.0 108.50 10.2 2
## 102 646.0 62.00 10.1 1
## 104 2108.0 128.65 11.1 2
## 106 1087.0 103.85 10.6 2
## 107 648.0 71.30 17.1 1
## 108 1062.0 49.60 10.3 2
## 127 1174.0 86.80 10.0 2
## 134 1353.0 127.10 10.7 2
## 135 836.0 74.40 11.0 2
## 136 1636.0 116.25 9.8 2
## 137 3472.0 155.00 10.1 2
## 141 794.0 125.55 10.6 2
## 142 720.0 134.85 9.9 2
## 147 791.0 57.35 10.7 2
## 150 1536.0 134.85 10.0 1
## 151 721.0 85.25 9.4 2
## 153 453.0 54.25 11.1 1
## 156 2468.0 201.50 11.5 2
## 157 1868.0 77.50 10.0 2
## 158 1408.0 116.25 11.2 2
## 161 964.0 120.90 9.6 2
## 166 1960.0 457.25 9.5 2
## 168 842.0 97.65 9.2 2
## 173 1403.0 89.90 9.8 2
## 174 663.0 79.05 9.7 1
## 177 559.0 43.40 10.4 2
## 182 674.0 26.35 9.9 2
## 183 1345.0 54.25 10.1 2
## 190 1828.0 133.30 10.2 2
## 194 975.0 125.55 9.1 2
## 197 1164.0 102.30 10.1 2
## 198 1692.0 141.05 9.6 2
## 202 1044.0 80.60 10.0 2
## 203 598.0 52.70 9.9 2
## 206 5300.0 57.35 11.0 1
## 209 1210.0 72.85 9.5 2
## 210 1098.0 122.45 10.0 2
## 218 897.0 66.65 10.1 1
## 220 2716.0 210.80 10.6 2
## 225 1065.0 88.35 10.1 2
## 227 1601.0 179.80 9.7 2
## 234 976.0 116.25 11.2 2
## 238 1093.0 128.65 9.8 2
## 251 676.0 83.00 9.9 2
## 257 377.0 56.00 10.6 2
## 258 1003.0 99.00 9.6 1
## 269 1428.0 188.00 10.1 2
## 272 674.0 58.00 9.8 1
## 276 640.0 145.00 10.7 2
## 279 911.0 57.00 9.8 2
## 285 688.0 57.00 10.0 1
## 288 637.0 117.00 9.6 2
## 289 960.0 86.00 9.8 2
## 302 657.0 52.00 10.6 2
## 306 1072.0 55.00 10.7 2
## 307 1282.0 130.00 10.5 2
## 308 797.0 91.00 10.4 2
## 309 1166.0 70.00 10.8 2
## 310 1713.0 171.00 10.2 2
## 311 1050.0 117.00 10.5 2
## 312 2115.0 136.00 10.8 2
#4 Summary of new dataframe.
summary(earlyPBC)
## id trt age sex
## Min. : 9.0 Min. :1.000 Min. :28.88 Length:83
## 1st Qu.: 96.5 1st Qu.:1.000 1st Qu.:41.45 Class :character
## Median :158.0 Median :1.000 Median :48.76 Mode :character
## Mean :164.1 Mean :1.434 Mean :48.68
## 3rd Qu.:226.0 3rd Qu.:2.000 3rd Qu.:56.40
## Max. :312.0 Max. :2.000 Max. :75.01
##
## ascites hepato bili chol
## Min. :0.0000 Min. :0.0000 Min. : 0.300 Min. : 127.0
## 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.: 0.600 1st Qu.: 246.2
## Median :0.0000 Median :0.0000 Median : 0.900 Median : 286.0
## Mean :0.0241 Mean :0.2289 Mean : 1.925 Mean : 338.2
## 3rd Qu.:0.0000 3rd Qu.:0.0000 3rd Qu.: 1.750 3rd Qu.: 392.8
## Max. :1.0000 Max. :1.0000 Max. :25.500 Max. :1480.0
## NA's :9
## albumin alk.phos ast protime
## Min. :2.640 Min. : 377.0 Min. : 26.35 Min. : 9.10
## 1st Qu.:3.430 1st Qu.: 720.5 1st Qu.: 68.54 1st Qu.: 9.90
## Median :3.650 Median : 1087.0 Median : 99.00 Median :10.30
## Mean :3.649 Mean : 1793.3 Mean :109.57 Mean :10.45
## 3rd Qu.:3.895 3rd Qu.: 1834.0 3rd Qu.:131.65 3rd Qu.:10.70
## Max. :4.640 Max. :13862.4 Max. :457.25 Max. :17.10
##
## stage
## Min. :1.000
## 1st Qu.:2.000
## Median :2.000
## Mean :1.807
## 3rd Qu.:2.000
## Max. :2.000
##
#4 Mean and median of AP & AST for patients with earlier stage PBC. Expect these to both be lower and they are.
e <- mean(earlyPBC$alk.phos, na.rm=TRUE)
f <- median(earlyPBC$alk.phos, na.rm=TRUE)
g <- mean(earlyPBC$ast, na.rm=TRUE)
h <- median(earlyPBC$ast, na.rm=TRUE)
earlypts <- c(e,f,g,h)
chart <- cbind(chart, earlypts)
chart
## enzyme allpts earlypts
## 1 AlkPhos Mean 1982.6558 1793.2554
## 2 AlkPhos Median 1259.0000 1087.0000
## 3 AST Mean 122.5563 109.5667
## 4 AST Median 114.7000 99.0000
#5 Rename values in the following columns
# ascites to be yes or none, sex to be spelled out male and female, treatment to be D-penicill, Placebo, or Not Random
earlyPBC$sex <- gsub("m","male", earlyPBC$sex)
earlyPBC$sex <- gsub("f","female", earlyPBC$sex)
earlyPBC$ascites <- gsub(0,"none", earlyPBC$ascites)
earlyPBC$ascites <- gsub(1,"yes", earlyPBC$ascites)
earlyPBC$trt <- gsub(1,"D-penicill", earlyPBC$trt)
earlyPBC$trt <- gsub(2,"placebo", earlyPBC$trt)
rename(earlyPBC,c("trt"="treatment", "bili"="bilirubin", chol="cholesterol" ))
## id treatment age sex ascites hepato bilirubin cholesterol albumin
## 9 9 D-penicill 42.50787 female none 0 3.2 562 3.08
## 24 24 D-penicill 44.52019 male none 1 2.1 456 4.00
## 25 25 placebo 45.07324 female none 0 0.7 298 4.10
## 29 29 placebo 63.87680 female none 0 0.7 370 3.78
## 31 31 placebo 41.55236 female none 1 4.7 296 3.44
## 34 34 D-penicill 52.06023 female none 0 0.8 364 3.70
## 36 36 placebo 56.41068 female none 0 0.3 172 3.39
## 43 43 D-penicill 48.87064 female none 0 1.1 361 3.64
## 45 45 placebo 41.79329 female none 0 0.6 NA 3.93
## 51 51 placebo 52.08761 female none 0 0.8 276 3.60
## 52 52 D-penicill 50.54073 male none 0 6.0 614 3.70
## 56 56 placebo 33.61807 female none 1 1.1 498 3.80
## 58 58 D-penicill 44.56947 male none 0 0.7 242 4.08
## 61 61 placebo 43.89870 male none 0 0.6 216 3.94
## 65 65 D-penicill 40.20260 female none 0 1.2 256 3.60
## 70 70 D-penicill 56.39973 female none 0 0.6 NA 4.64
## 73 73 placebo 38.49418 female none 0 0.7 132 3.60
## 89 89 D-penicill 52.44353 female none 1 2.0 408 3.65
## 90 90 D-penicill 33.47570 male none 0 1.6 660 4.22
## 93 93 D-penicill 36.53388 female none 0 1.3 353 3.67
## 95 95 placebo 46.39014 female yes 1 17.4 NA 2.64
## 98 98 D-penicill 28.88433 female none 0 1.0 239 3.77
## 99 99 placebo 48.46817 male none 0 1.8 460 3.35
## 102 102 D-penicill 56.56947 female none 0 0.9 248 3.97
## 104 104 D-penicill 43.01711 female none 0 1.1 303 3.64
## 106 106 D-penicill 68.50924 female yes 1 2.1 NA 3.90
## 107 107 placebo 62.52156 female none 0 0.6 212 4.03
## 108 108 D-penicill 50.35729 female none 0 0.4 127 3.50
## 127 127 placebo 44.10404 female none 0 0.5 268 4.08
## 134 134 placebo 42.33539 female none 0 0.7 578 3.67
## 135 135 D-penicill 42.96783 female none 0 0.4 263 3.57
## 136 136 D-penicill 55.96167 female none 0 0.8 263 3.35
## 137 137 D-penicill 62.86105 female none 0 1.1 399 3.60
## 141 141 D-penicill 47.03628 female none 0 0.9 346 3.77
## 142 142 placebo 55.72621 female none 1 1.0 364 3.48
## 147 147 D-penicill 75.01164 female none 0 1.2 288 3.37
## 150 150 placebo 34.98700 female none 0 1.0 NA 3.63
## 151 151 D-penicill 55.04175 female none 0 0.9 460 3.03
## 153 153 D-penicill 49.60438 female none 0 0.5 217 3.85
## 156 156 placebo 59.40862 female none 1 25.5 358 3.52
## 157 157 placebo 48.75838 female none 0 0.6 286 3.42
## 158 158 D-penicill 36.49281 female none 0 3.4 450 3.37
## 161 161 placebo 42.74333 female none 0 2.3 502 3.56
## 166 166 placebo 41.35524 female none 1 5.7 1480 3.26
## 168 168 placebo 47.75359 female none 1 0.4 257 3.80
## 173 173 D-penicill 30.27515 female none 1 3.0 236 3.42
## 174 174 D-penicill 55.56742 female none 0 0.5 NA 3.85
## 177 177 placebo 55.45243 female none 0 0.9 258 4.01
## 182 182 D-penicill 56.94456 female none 1 0.6 NA 3.69
## 183 183 placebo 40.26010 female none 0 0.5 201 3.73
## 190 190 D-penicill 54.52704 female none 0 2.3 NA 3.93
## 194 194 D-penicill 44.39699 female none 0 1.5 293 4.30
## 197 197 D-penicill 44.62697 female none 0 0.6 266 3.97
## 198 198 placebo 35.79740 female none 0 0.7 286 2.90
## 202 202 D-penicill 61.63997 female none 0 0.5 223 3.80
## 203 203 placebo 37.05681 female none 1 0.5 149 4.04
## 206 206 D-penicill 61.99042 female none 0 0.6 213 4.07
## 209 209 placebo 52.62423 female none 1 0.7 252 4.01
## 210 210 placebo 49.76318 male none 1 0.9 346 3.37
## 218 218 D-penicill 34.59548 female none 0 0.5 NA 3.89
## 220 220 D-penicill 63.63039 female none 0 3.8 426 3.22
## 225 225 D-penicill 38.77070 female none 0 0.7 274 3.66
## 227 227 D-penicill 58.95140 male none 0 2.3 316 3.35
## 234 234 placebo 34.37919 female none 0 0.9 578 3.35
## 238 238 placebo 66.41205 female none 0 0.6 280 3.35
## 251 251 D-penicill 56.97467 female none 0 0.5 227 3.61
## 257 257 D-penicill 58.26420 female none 0 0.5 280 4.23
## 258 258 D-penicill 51.48802 female none 0 1.1 414 3.44
## 269 269 placebo 45.08419 female none 0 3.6 374 3.50
## 272 272 D-penicill 38.39836 female none 0 0.5 226 2.93
## 276 276 D-penicill 50.10815 female none 0 1.0 219 3.85
## 279 279 placebo 56.15332 female none 0 0.5 198 3.77
## 285 285 D-penicill 46.34908 female none 0 0.8 253 3.48
## 288 288 placebo 48.93634 female none 1 8.7 310 3.89
## 289 289 D-penicill 67.57290 male none 1 4.0 416 3.99
## 302 302 D-penicill 37.79877 female none 0 0.7 335 3.95
## 306 306 placebo 56.39151 female none 1 0.6 239 3.45
## 307 307 placebo 30.57358 female none 0 0.8 273 3.56
## 308 308 D-penicill 61.18275 female none 1 0.4 246 3.58
## 309 309 placebo 58.29979 female none 0 0.4 260 2.75
## 310 310 D-penicill 62.33265 female none 0 1.7 434 3.35
## 311 311 D-penicill 37.99863 female none 0 2.0 247 3.16
## 312 312 placebo 33.15264 female none 0 6.4 576 3.79
## alk.phos ast protime stage
## 9 2276.0 144.15 11.0 2
## 24 5719.0 221.88 9.9 2
## 25 661.0 106.95 11.3 2
## 29 5833.0 73.53 10.6 2
## 31 9933.2 206.40 10.3 2
## 34 1840.0 170.50 10.5 2
## 36 558.0 71.30 10.6 2
## 43 5430.2 67.08 10.6 2
## 45 1826.0 71.30 10.9 2
## 51 4332.0 99.33 10.6 2
## 52 5084.4 206.40 10.6 1
## 56 13862.4 95.46 10.6 2
## 58 5890.0 56.76 10.6 1
## 61 601.0 60.45 13.0 1
## 65 724.0 141.05 10.0 1
## 70 666.0 54.25 10.6 2
## 73 423.0 49.60 11.0 1
## 89 1083.0 110.05 11.4 2
## 90 1857.0 151.90 11.0 2
## 93 2039.0 232.50 11.1 2
## 95 559.0 119.35 11.7 2
## 98 1877.0 97.65 10.2 1
## 99 1472.0 108.50 10.2 2
## 102 646.0 62.00 10.1 1
## 104 2108.0 128.65 11.1 2
## 106 1087.0 103.85 10.6 2
## 107 648.0 71.30 17.1 1
## 108 1062.0 49.60 10.3 2
## 127 1174.0 86.80 10.0 2
## 134 1353.0 127.10 10.7 2
## 135 836.0 74.40 11.0 2
## 136 1636.0 116.25 9.8 2
## 137 3472.0 155.00 10.1 2
## 141 794.0 125.55 10.6 2
## 142 720.0 134.85 9.9 2
## 147 791.0 57.35 10.7 2
## 150 1536.0 134.85 10.0 1
## 151 721.0 85.25 9.4 2
## 153 453.0 54.25 11.1 1
## 156 2468.0 201.50 11.5 2
## 157 1868.0 77.50 10.0 2
## 158 1408.0 116.25 11.2 2
## 161 964.0 120.90 9.6 2
## 166 1960.0 457.25 9.5 2
## 168 842.0 97.65 9.2 2
## 173 1403.0 89.90 9.8 2
## 174 663.0 79.05 9.7 1
## 177 559.0 43.40 10.4 2
## 182 674.0 26.35 9.9 2
## 183 1345.0 54.25 10.1 2
## 190 1828.0 133.30 10.2 2
## 194 975.0 125.55 9.1 2
## 197 1164.0 102.30 10.1 2
## 198 1692.0 141.05 9.6 2
## 202 1044.0 80.60 10.0 2
## 203 598.0 52.70 9.9 2
## 206 5300.0 57.35 11.0 1
## 209 1210.0 72.85 9.5 2
## 210 1098.0 122.45 10.0 2
## 218 897.0 66.65 10.1 1
## 220 2716.0 210.80 10.6 2
## 225 1065.0 88.35 10.1 2
## 227 1601.0 179.80 9.7 2
## 234 976.0 116.25 11.2 2
## 238 1093.0 128.65 9.8 2
## 251 676.0 83.00 9.9 2
## 257 377.0 56.00 10.6 2
## 258 1003.0 99.00 9.6 1
## 269 1428.0 188.00 10.1 2
## 272 674.0 58.00 9.8 1
## 276 640.0 145.00 10.7 2
## 279 911.0 57.00 9.8 2
## 285 688.0 57.00 10.0 1
## 288 637.0 117.00 9.6 2
## 289 960.0 86.00 9.8 2
## 302 657.0 52.00 10.6 2
## 306 1072.0 55.00 10.7 2
## 307 1282.0 130.00 10.5 2
## 308 797.0 91.00 10.4 2
## 309 1166.0 70.00 10.8 2
## 310 1713.0 171.00 10.2 2
## 311 1050.0 117.00 10.5 2
## 312 2115.0 136.00 10.8 2
head(earlyPBC)
## id trt age sex ascites hepato bili chol albumin alk.phos
## 9 9 D-penicill 42.50787 female none 0 3.2 562 3.08 2276.0
## 24 24 D-penicill 44.52019 male none 1 2.1 456 4.00 5719.0
## 25 25 placebo 45.07324 female none 0 0.7 298 4.10 661.0
## 29 29 placebo 63.87680 female none 0 0.7 370 3.78 5833.0
## 31 31 placebo 41.55236 female none 1 4.7 296 3.44 9933.2
## 34 34 D-penicill 52.06023 female none 0 0.8 364 3.70 1840.0
## ast protime stage
## 9 144.15 11.0 2
## 24 221.88 9.9 2
## 25 106.95 11.3 2
## 29 73.53 10.6 2
## 31 206.40 10.3 2
## 34 170.50 10.5 2