{"id":2280,"date":"2021-02-01T01:01:12","date_gmt":"2021-01-31T16:01:12","guid":{"rendered":"http:\/\/141.164.34.82\/?p=2280"},"modified":"2021-02-01T01:55:39","modified_gmt":"2021-01-31T16:55:39","slug":"shingles-lattice","status":"publish","type":"post","link":"http:\/\/ds.sumeun.org\/?p=2280","title":{"rendered":"Shingles & Lattice"},"content":{"rendered":"
library(lattice)\ndata('BankWages', package='AER')\ndata(mtcars)\n<\/code><\/pre>\n
\uac04\ub2e8\ud558\uac8c \uc870\uac74\ubd80 \uc774\ubcc0\uc218 \uc0b0\uc810\ub3c4 \uadf8\ub9ac\uae30<\/h2>\n
xyplot(x ~ y | g, dat)\n<\/code><\/pre>\n
    \n
  • \ub9cc\uc57d \uc870\uac74 \ubcc0\uc218 g<\/code> \uc5ed\uc2dc \uc5f0\uc18d\ud615\uc774\ub77c\uba74, \uadf8 \uacb0\uacfc\ub294 \ub2e4\uc74c\uacfc \uac19\ub2e4.<\/li>\n<\/ul>\n
    xyplot(qsec ~ hp | mpg, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    \ubaa8\ub4e0 mpg<\/code> \uac12\uc5d0 \ub300\ud574 qsec ~ hp<\/code>\uc758 \uc0b0\uc810\ub3c4\ub97c \uadf8\ub9ac\uace0 \uc788\uc5b4\uc11c \uadf8\ub2e4\uc9c0 \uc720\uc6a9\ud558\uc9c0 \uc54a\ub2e4.<\/p>\n
    table(mtcars$mpg)\n<\/code><\/pre>\n
    ## \n## 10.4 13.3 14.3 14.7   15 15.2 15.5 15.8 16.4 17.3 17.8 18.1 18.7 19.2 19.7   21 21.4 21.5 22.8 24.4   26 27.3 30.4 \n##    2    1    1    1    1    2    1    1    1    1    1    1    1    2    1    2    2    1    2    1    1    1    2 \n## 32.4 33.9 \n##    1    1\n<\/pre>\n
    \ub3d9\uc218\ub85c \uad6c\uac04 \ub098\ub204\uae30<\/h2>\n
    \ubcc0\uc218 mpg<\/code>\ub97c \uad6c\uac04\ubcc4\ub85c \ub098\ub220\uc11c qsec\uc640 hp\uc758 \uc0b0\uc810\ub3c4\ub97c \uadf8\ub9ac\uace0\uc790 \ud55c\ub2e4\uba74
    \n\ub2e4\uc74c\uacfc \uac19\uc774 \ud560 \uc218 \uc788\ub2e4.<\/p>\n
    # equal.count(mtcars$mpg,  # \uad6c\uac04 \ubcc0\uc218\n#             number=3,    # \uad6c\uac04 \uac2f\uc218\n#             overlap=0)   # \uad6c\uac04\uc774 \uacb9\uce58\ub294 \uc815\ub3c4\nmpgequal <- equal.count(mtcars$mpg, number=3, overlap=0)\nxyplot(qsec ~ hp | mpgequal, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    \uc138 \uad6c\uac04\uc774 \uc124\uc815\ub418\uc5c8\uc73c\uba74, \ud55c \uad6c\uac04 \uc548\uc5d0\ub294 \uac70\uc758 \ub3d9\uc77c\ud55c \uac2f\uc218\uc758 \uc810\uc774 \ud45c\uc2dc\ub41c\ub2e4.
    \n\uc2e4\uc81c \ud45c\uc2dc\ub418\ub294 \uc810\uc758 \uac2f\uc218\ub294 mpgequal \ubcc0\uc218\ub97c \uc0b4\ud3b4\ubd04\uc73c\ub85c\uc368 \ud655\uc778 \uac00\ub2a5\ud558\ub2e4.<\/p>\n
    mpgequal\n<\/code><\/pre>\n
    ## \n## Data:\n##  [1] 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3 15.2 10.4 10.4 14.7 32.4 30.4 33.9 21.5 15.5\n## [23] 15.2 13.3 19.2 27.3 26.0 30.4 15.8 19.7 15.0 21.4\n## \n## Intervals:\n##     min   max count\n## 1 10.35 16.45    11\n## 2 17.25 21.45    11\n## 3 21.35 33.95    12\n## \n## Overlap between adjacent intervals:\n## [1] 0 2\n<\/pre>\n
    mpg<\/code> = 10.35-16.45 \uad6c\uac04\uc5d0\ub294 11\uac1c\uc758 \uc810\uc774 \ub098\ud0c0\ub0a8\uc744 \ud655\uc778\ud560 \uc218 \uc788\ub2e4. \uad6c\uac04\uc774 \uacb9\uce58\ub294 \uacbd\uc6b0\ub3c4 \uc788\ub294\ub370 \uad6c\uac042(17.25-21.45)\uc640 \uad6c\uac043(21.35-33.95)\ub294 \uacb9\uce58\ub294 \uc810\uc774 \uc788\uace0, \uac2f\uc218\ub294 2\uac1c\uc774\ub2e4(Overlap between adjacent intervals).<\/p>\n
    \uad6c\uac04\uc758 \uac2f\uc218\ub97c 5\uac1c\ub85c \ub298\ub9ac\uace0, overlap=<\/code> \uc5ed\uc2dc 0.2\ub85c \ub298\ub824\ubcf4\uba74 \ub2e4\uc74c\uacfc \uac19\uc740 \uacb0\uacfc\ub97c \uc5bb\ub294\ub2e4.<\/p>\n
    mpgequal <- equal.count(mtcars$mpg, number=5, overlap=0.2)\nxyplot(qsec ~ hp | mpgequal, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    mpgequal\n<\/code><\/pre>\n
    ## \n## Data:\n##  [1] 21.0 21.0 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 17.8 16.4 17.3 15.2 10.4 10.4 14.7 32.4 30.4 33.9 21.5 15.5\n## [23] 15.2 13.3 19.2 27.3 26.0 30.4 15.8 19.7 15.0 21.4\n## \n## Intervals:\n##     min   max count\n## 1 10.35 15.25     8\n## 2 15.15 18.15     8\n## 3 17.75 21.05     8\n## 4 20.95 24.45     8\n## 5 22.75 33.95     9\n## \n## Overlap between adjacent intervals:\n## [1] 2 2 2 3\n<\/pre>\n
    \ub4f1\uad6c\uac04 \uad6c\ud68d<\/h2>\n
    \uad6c\uac04\uc758 \uae38\uc774\ub97c \uc77c\uc815\ud558\uac8c \ud558\ub294 \uac8c \ub354 \ub098\uc740 \uacbd\uc6b0\ub3c4 \uc788\ub2e4. \uc774\ub7f4 \ub54c \uc190\uc27d\uac8c \uc4f8 \uc218 \uc788\ub294 \ubc29\ubc95\uc740 \ub2e4\uc74c\uacfc \uac19\ub2e4.<\/p>\n
    require(dplyr)\nmtcars <- mtcars %>% mutate(mpgCut = cut(mpg, breaks=5, include.lowest=TRUE))\nxyplot(qsec ~ hp | mpgCut, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    cut<\/code> \ud568\uc218\ub97c \uc4f0\uba74 \uc190\uc27d\uac8c \ub4f1\uad6c\uac04\uc73c\ub85c \uad6c\ud68d\uc774 \uac00\ub2a5\ud558\ub2e4. \uc774\ub54c include.lowest=TRUE<\/code>\ub97c \uc78a\uc9c0 \ub9d0\uc790. (include.lowest=FALSE<\/code>\ub97c \ud558\uba74 \uad6c\uac04\uc774 (b1', b2], (b2,b3], …, (bn-1, bn] \uc774 \ub418\uace0,  include.lowest=TRUE<\/code>\ub97c \ud558\uba74 \uad6c\uac04\uc774 [b1, b2], (b2, b3], … (bn-1, bn]\uc774 \ub41c\ub2e4. \uc774\ub54c b1'\ub294 b1\ubcf4\ub2e4 \uc870\uae08 \uc791\uc544\uc9c4\ub2e4.)<\/p>\n
    \ub9cc\uc57d \uad6c\uac04 \uc0ac\uc774\uc5d0 \uc870\uae08\uc758 overlap<\/code>(\uacb9\uce68)\uc774 \ud544\uc694\ud55c \uacbd\uc6b0\ub294 \uc5b4\ub5bb\uac8c \ud574\uc57c \ud560\uae4c? \uc774\ubbf8 \ub9c8\ub828\ub418\uc5b4 \uc788\ub294 \ud568\uc218\ub85c\ub294 \ub2e4\uc18c \uae4c\ub2e4\ub85c\uc6b4 \ubd80\ubd84\uc774 \uc788\uc5b4\uc11c \ub2e4\uc74c\uc758 equal.interval()<\/code> \ud568\uc218\ub97c \ub9cc\ub4e4\uc5c8\ub2e4.<\/p>\n
    equal.interval = function(v, breaks=10, na.rm=FALSE, overlap=0) {\n  stopifnot(inherits(breaks, 'numeric') & \n              inherits(v, 'numeric'))\n  if (length(breaks) == 1) {\n    int = (max(v, na.rm=na.rm)-min(v, na.rm=na.rm))\/breaks\n    br1 = min(v, na.rm=na.rm) + seq(0,breaks-1)*int - overlap\/2*int\n    br2 = min(v, na.rm=na.rm) + seq(1,breaks)*int + overlap\/2*int\n    return(shingle(v, cbind(br1,br2)))\n  } else {\n    br1 = breaks[-length(breaks)]\n    br2 = breaks[-1]\n    int = br2 - br1\n    br2 = br2 + overlap\/2*int\n    br1 = br1 -overlap\/2*int\n    #print(int)\n    #print(br1)\n    #print(br2)\n    return(shingle(v, cbind(br1, br2)))\n  }\n}\n<\/code><\/pre>\n
    \uad6c\uac04\uc744 \uadf8\ub824\ubcf4\uba74 \ub2e4\uc74c\uacfc \uac19\ub2e4.<\/p>\n
    ## Testing\nx <- rnorm(1000)\neqi1 = equal.interval(x, breaks=3, overlap=0.1)\neqi2 = equal.interval(x, breaks=c(-1,0,1), overlap=0.1)\n\nplot(eqi1)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    plot(eqi2)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    \uc774\uc81c \uc2e4\uc81c \ub370\uc774\ud130\uc5d0 \uc801\uc6a9\ud574 \ubcf4\uc790. \ub3d9\uc77c\ud55c \uc870\uac74\uc5d0 overlap=<\/code>\uc744 \ub298\ub824\ub098\uac00\ub2c8 mpg<\/code>\uc5d0 \uc5b4\ub290 \uc815\ub3c4 \uc77c\uc815\ud558\uac8c \uace0\uc815\ud588\uc744 \ub54c qsec<\/code>\uc640 hp<\/code>\uc758 \uad00\uacc4 \uc880 \ub354 \ub69c\ub837\ud558\uac8c \ubcf4\uc774\ub294 \ub4ef \ud558\ub2e4(?).<\/p>\n
    mpgInt = equal.interval(mtcars$mpg, breaks = 6, overlap=0.1)\nxyplot(qsec ~ hp | mpgInt, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    mpgInt = equal.interval(mtcars$mpg, breaks = 6, overlap=0.2)\nxyplot(qsec ~ hp | mpgInt, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n
    mpgInt = equal.interval(mtcars$mpg, breaks = 6, overlap=0.3)\nxyplot(qsec ~ hp | mpgInt, mtcars)\n<\/code><\/pre>\n
    \"plot<\/p>\n","protected":false},"excerpt":{"rendered":"
    library(lattice) data('BankWages', package='AER') data(mtcars) \uac04\ub2e8\ud558\uac8c \uc870\uac74\ubd80 \uc774\ubcc0\uc218 \uc0b0\uc810\ub3c4 \uadf8\ub9ac\uae30 xyplot(x ~ y | g, dat) \ub9cc\uc57d \uc870\uac74 \ubcc0\uc218 g \uc5ed\uc2dc \uc5f0\uc18d\ud615\uc774\ub77c\uba74, \uadf8 \uacb0\uacfc\ub294 \ub2e4\uc74c\uacfc \uac19\ub2e4. xyplot(qsec ~ hp | mpg, mtcars) \ubaa8\ub4e0 mpg \uac12\uc5d0 \ub300\ud574 qsec ~ hp\uc758 \uc0b0\uc810\ub3c4\ub97c \uadf8\ub9ac\uace0 \uc788\uc5b4\uc11c \uadf8\ub2e4\uc9c0 \uc720\uc6a9\ud558\uc9c0 \uc54a\ub2e4. table(mtcars$mpg) ## ## 10.4 13.3 14.3 14.7 15 15.2 15.5 15.8 […]<\/p>\n","protected":false},"author":1,"featured_media":2283,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":[],"categories":[28,165],"tags":[169,477],"jetpack_featured_media_url":"http:\/\/ds.sumeun.org\/wp-content\/uploads\/2021\/02\/Shingles.png","_links":{"self":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts\/2280"}],"collection":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=2280"}],"version-history":[{"count":6,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts\/2280\/revisions"}],"predecessor-version":[{"id":2287,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts\/2280\/revisions\/2287"}],"wp:featuredmedia":[{"embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/media\/2283"}],"wp:attachment":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=2280"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=2280"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=2280"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}