Skip to content

Commit

Permalink
update: tidy
Browse files Browse the repository at this point in the history
  • Loading branch information
@shanmdphd
shanmdphd committed Dec 27, 2018
1 parent f1efd01 commit c57b249
Show file tree
Hide file tree
Showing 12,141 changed files with 1,281 additions and 2,386,612 deletions.
The diff you're trying to view is too large. We only load the first 3000 changed files.
2 changes: 1 addition & 1 deletion .gitignore
View file
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
lib.zip
result

# History files
.Rhistory
Expand Down
Binary file added Edison 발표 포스터.pdf
View file
Binary file not shown.
Binary file added Edison을 통한 비선형 약동학.pptx
View file
Binary file not shown.
8 changes: 1 addition & 7 deletions PK02-report.Rmd
View file
Original file line number Diff line number Diff line change
Expand Up @@ -16,11 +16,5 @@ knitr::opts_chunk$set(error = TRUE)
소스 코드는 [깃헙](https://github.com/asancpt/edison-gab)에 올라와 있습니다.
에디슨 앱은 <https://www.edison.re.kr/simulation> 에서 확인할 수 있습니다.

```{r, code = readLines('R/PK02.R')}
```{r pk02-figures, code = readLines('PK02.R')}
```

## Figures

![](figure1.png)

![](figure2.png)
267 changes: 244 additions & 23 deletions PK02-report.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -53,13 +53,16 @@ plot(DV2 ~ TIME, data=PK02, type="o")
lines(DV ~ TIME, data=PK02, type="o", col="red")
```

![plot of chunk unnamed-chunk-1](figure/unnamed-chunk-1-1.svg)
![plot of chunk pk02-figures](figure/pk02-figures-1.png)

```r
dev.new()
plot(log(DV2) ~ TIME, data=PK02, type="o")
lines(log(DV) ~ TIME, data=PK02, type="o", col="red")
```

![plot of chunk pk02-figures](figure/pk02-figures-2.png)

```r
## NCA
library(NonCompart)
R1 = sNCA(PK02$TIME, PK02$DV2, dose=100, adm="Bolus", doseUnit="ug", timeUnit="min") ; R1
Expand Down Expand Up @@ -93,7 +96,9 @@ R1 = sNCA(PK02$TIME, PK02$DV2, dose=100, adm="Bolus", doseUnit="ug", timeUnit="m
## [26] "min2*ug/L" "min2*ug/L" "%" "%" "ug/L"
## [31] "%" "%" "L" "L" "L/min"
## [36] "L/min" "min" "min" "min" "L"
## [41] "L"
## [41] "L"
## attr(,"UsedPoints")
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
```

```r
Expand Down Expand Up @@ -125,7 +130,9 @@ R2 = sNCA(PK02$TIME, PK02$DV, dose=100, adm="Extravascular", doseUnit="ug", time
## [21] "min*ug/L" "min*ug/L/ug" "%" "%" "min2*ug/L"
## [26] "min2*ug/L" "min2*ug/L" "%" "%" "L"
## [31] "L" "L/min" "L/min" "min" "min"
## [36] "min"
## [36] "min"
## attr(,"UsedPoints")
## [1] 8 9 10 11 12 13 14 15 16 17
```

```r
Expand Down Expand Up @@ -157,19 +164,19 @@ fPK01 = function(THETA) # Prediction function
return(cbind(F, H1))
}

nlr(fPK01, dPK02, pNames=c("Ka", "Ka", "V"), IE=c(0.05, 0.1, 30))
nlr(fPK01, dPK02, pNames=c("k", "Ka", "V"), IE=c(0.05, 0.1, 30))
```

```
## $Est
## Ka Ka V AddErrVar AddErrSD
## k Ka V AddErrVar AddErrSD
## PE 0.013698916 0.013699384 12.512458 0.27325141 0.52273455
## SE 0.009900321 0.009900641 8.994196 0.06627321 0.06339088
## RSE 72.270834504 72.270699963 71.881932 24.25356542 12.12678271
##
## $Cov
## Ka Ka V AddErrVar
## Ka 9.801636e-05 -9.172333e-05 -8.727111e-02 -3.783325e-09
## k Ka V AddErrVar
## k 9.801636e-05 -9.172333e-05 -8.727111e-02 -3.783325e-09
## Ka -9.172333e-05 9.802268e-05 8.603379e-02 3.734622e-09
## V -8.727111e-02 8.603379e-02 8.089557e+01 3.532174e-06
## AddErrVar -3.783325e-09 3.734622e-09 3.532174e-06 4.392138e-03
Expand All @@ -188,20 +195,78 @@ nlr(fPK01, dPK02, pNames=c("Ka", "Ka", "V"), IE=c(0.05, 0.1, 30))
## [1] 7.541375e-05
##
##
## $`Objective Function Value`
## [1] -10.11035
##
## $`-2LL`
## [1] 21.13356
##
## $AIC
## [1] 29.13356
##
## $AICc
## [1] 32.4669
##
## $BIC
## [1] 32.46642
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## F H1
## [1,] 0.5861090 1
## [2,] 0.8209607 1
## [3,] 1.0221479 1
## [4,] 1.3369352 1
## [5,] 1.5543702 1
## [6,] 1.6942200 1
## [7,] 1.7727864 1
## [8,] 1.7630462 1
## [9,] 1.7081521 1
## [10,] 1.5585419 1
## [11,] 1.3825345 1
## [12,] 1.2013751 1
## [13,] 1.0276420 1
## [14,] 0.7948857 1
## [15,] 0.6022997 1
## [16,] 0.3309442 1
## [17,] 0.1745694 1
##
## $Residual
## F H1
## [1,] -0.586109032 -1.00
## [2,] -0.540960661 -0.72
## [3,] -0.482147939 -0.46
## [4,] -0.086935218 0.25
## [5,] 0.445629849 1.00
## [6,] 0.335779984 1.03
## [7,] 0.177213639 0.95
## [8,] 0.086953758 0.85
## [9,] 0.011847892 0.72
## [10,] 0.041458101 0.60
## [11,] -0.042534465 0.34
## [12,] -0.111375051 0.09
## [13,] -0.167641953 -0.14
## [14,] -0.014885652 -0.22
## [15,] -0.002299701 -0.40
## [16,] -0.120944178 -0.79
## [17,] 0.005430617 -0.82
##
## $`Elapsed Time`
## Time difference of 0.02900195 secs
```

```r
wnl5(fPK01, dPK02, pNames=c("Ka", "Ka", "V"), IE=c(0.05, 0.1, 30))
wnl5(fPK01, dPK02, pNames=c("k", "Ka", "V"), IE=c(0.05, 0.1, 30))
```

```
## $PE
## Ka Ka V
## k Ka V
## 0.01369891 0.01369939 12.51246734
##
## $WRSS
Expand All @@ -226,6 +291,58 @@ wnl5(fPK01, dPK02, pNames=c("Ka", "Ka", "V"), IE=c(0.05, 0.1, 30))
##
## $SBC
## [1] 46.3926
##
## $`Condition Number`
## [1] 4823.704
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## F H1
## [1,] 0.5861090 1
## [2,] 0.8209606 1
## [3,] 1.0221479 1
## [4,] 1.3369351 1
## [5,] 1.5543700 1
## [6,] 1.6942198 1
## [7,] 1.7727861 1
## [8,] 1.7630460 1
## [9,] 1.7081518 1
## [10,] 1.5585416 1
## [11,] 1.3825342 1
## [12,] 1.2013748 1
## [13,] 1.0276417 1
## [14,] 0.7948854 1
## [15,] 0.6022995 1
## [16,] 0.3309441 1
## [17,] 0.1745693 1
##
## $Residual
## F H1
## [1,] -0.586108992 -1.00
## [2,] -0.540960601 -0.72
## [3,] -0.482147858 -0.46
## [4,] -0.086935098 0.25
## [5,] 0.445630005 1.00
## [6,] 0.335780173 1.03
## [7,] 0.177213856 0.95
## [8,] 0.086954033 0.85
## [9,] 0.011848176 0.72
## [10,] 0.041458394 0.60
## [11,] -0.042534175 0.34
## [12,] -0.111374772 0.09
## [13,] -0.167641692 -0.14
## [14,] -0.014885424 -0.22
## [15,] -0.002299508 -0.40
## [16,] -0.120944050 -0.79
## [17,] 0.005430696 -0.82
##
## $`Elapsed Time`
## Time difference of 0.007999897 secs
```

```r
Expand All @@ -249,25 +366,25 @@ fPK02 = function(THETA) # Prediction function
BA = 1
Dose = 100
colnames(PK02) = c("TIME", "DV", "DV2")
nlr(fPK02, dPK02, pNames=c("Ka", "Ka", "V", "tlag"), IE=c(0.05, 0.1, 30, 20))
nlr(fPK02, dPK02, pNames=c("k", "Ka", "V", "tlag"), IE=c(0.05, 0.1, 30, 20))
```

```
## $Est
## Ka Ka V tlag AddErrVar AddErrSD
## k Ka V tlag AddErrVar AddErrSD
## PE 0.010110127 0.03021029 28.45123 11.495191 0.24179045 0.49172192
## SE 0.006120683 0.02122374 13.06396 3.189944 0.05864321 0.05963046
## RSE 60.540118876 70.25332892 45.91703 27.750249 24.25373334 12.12686667
##
## $Cov
## Ka Ka V tlag
## Ka 3.746276e-05 -1.174219e-04 -7.746652e-02 -3.191263e-03
## k Ka V tlag
## k 3.746276e-05 -1.174219e-04 -7.746652e-02 -3.191263e-03
## Ka -1.174219e-04 4.504470e-04 2.590532e-01 2.241856e-02
## V -7.746652e-02 2.590532e-01 1.706671e+02 8.061794e+00
## tlag -3.191263e-03 2.241856e-02 8.061794e+00 1.017574e+01
## AddErrVar -1.856686e-09 2.585447e-09 4.069133e-06 -4.016032e-06
## AddErrVar
## Ka -1.856686e-09
## k -1.856686e-09
## Ka 2.585447e-09
## V 4.069133e-06
## tlag -4.016032e-06
Expand All @@ -287,20 +404,78 @@ nlr(fPK02, dPK02, pNames=c("Ka", "Ka", "V", "tlag"), IE=c(0.05, 0.1, 30, 20))
## [1] 0.02721811
##
##
## $`Objective Function Value`
## [1] -14.26949
##
## $`-2LL`
## [1] 16.97442
##
## $AIC
## [1] 26.97442
##
## $AICc
## [1] 32.42897
##
## $BIC
## [1] 31.14049
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## F H1
## [1,] -0.1636281 1
## [2,] 0.3468341 1
## [3,] 0.7616992 1
## [4,] 1.3608829 1
## [5,] 1.7271226 1
## [6,] 1.9285323 1
## [7,] 2.0147569 1
## [8,] 1.8956642 1
## [9,] 1.7945198 1
## [10,] 1.5645553 1
## [11,] 1.3319948 1
## [12,] 1.1175826 1
## [13,] 0.9290736 1
## [14,] 0.6968753 1
## [15,] 0.5189470 1
## [16,] 0.2849555 1
## [17,] 0.1556878 1
##
## $Residual
## F H1
## [1,] 0.163628080 -1.00
## [2,] -0.066834101 -0.72
## [3,] -0.221699228 -0.46
## [4,] -0.110882941 0.25
## [5,] 0.272877394 1.00
## [6,] 0.101467695 1.03
## [7,] -0.064756854 0.95
## [8,] -0.045664238 0.85
## [9,] -0.074519831 0.72
## [10,] 0.035444654 0.60
## [11,] 0.008005216 0.34
## [12,] -0.027582587 0.09
## [13,] -0.069073596 -0.14
## [14,] 0.083124750 -0.22
## [15,] 0.081052985 -0.40
## [16,] -0.074955540 -0.79
## [17,] 0.024312152 -0.82
##
## $`Elapsed Time`
## Time difference of 0.03200197 secs
```

```r
wnl5(fPK02, dPK02, pNames=c("Ka", "Ka", "V", "tlag"), IE=c(0.05, 0.1, 30, 20))
wnl5(fPK02, dPK02, pNames=c("k", "Ka", "V", "tlag"), IE=c(0.05, 0.1, 30, 20))
```

```
## $PE
## Ka Ka V tlag
## k Ka V tlag
## 0.01011006 0.03021076 28.45136106 11.49550270
##
## $WRSS
Expand All @@ -325,10 +500,56 @@ wnl5(fPK02, dPK02, pNames=c("Ka", "Ka", "V", "tlag"), IE=c(0.05, 0.1, 30, 20))
##
## $SBC
## [1] 47.14624
##
## $`Condition Number`
## [1] 8669.947
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## F H1
## [1,] -0.1636651 1
## [2,] 0.3468090 1
## [3,] 0.7616832 1
## [4,] 1.3608787 1
## [5,] 1.7271244 1
## [6,] 1.9285368 1
## [7,] 2.0147622 1
## [8,] 1.8956679 1
## [9,] 1.7945227 1
## [10,] 1.5645571 1
## [11,] 1.3319960 1
## [12,] 1.1175836 1
## [13,] 0.9290747 1
## [14,] 0.6968767 1
## [15,] 0.5189487 1
## [16,] 0.2849573 1
## [17,] 0.1556894 1
##
## $Residual
## F H1
## [1,] 0.163665132 -1.00
## [2,] -0.066808992 -0.72
## [3,] -0.221683184 -0.46
## [4,] -0.110878674 0.25
## [5,] 0.272875570 1.00
## [6,] 0.101463167 1.03
## [7,] -0.064762169 0.95
## [8,] -0.045667864 0.85
## [9,] -0.074522695 0.72
## [10,] 0.035442918 0.60
## [11,] 0.008004036 0.34
## [12,] -0.027583634 0.09
## [13,] -0.069074745 -0.14
## [14,] 0.083123304 -0.22
## [15,] 0.081051291 -0.40
## [16,] -0.074957335 -0.79
## [17,] 0.024310629 -0.82
##
## $`Elapsed Time`
## Time difference of 0.00999999 secs
```

## Figures

![](figure1.png)

![](figure2.png)
Loading

0 comments on commit c57b249

Please sign in to comment.