Blog

  • 평균치 비교

    이번 숙제에 R을 사용해 보기 위해서 계속 공부중이다. Odds ratio 를 구하는 함수는 없는 것 같지만 계속 찾아볼 예정이고 평균치 비교는 비교적 쉽게 찾을 수 있었다. 이건 비통계전공자의 전투족에 해당하는 내용이다. 🙂

    기본적으로 함수는 t.test 이다. 간단하다. 그래서 어렵다. 아파트나 자동차와 마찬가지로 옵션을 넣어주어야 한다.
    사용하는 기본 옵션은 다음과 같다.
    var.equal = T(RUE) or F(ALSE); 두 집단이 등분산이면 TRUE를 아니면 FALSE를 선택한다.

    등분산 여부는 사전에 var.test를 사용해서 확인한다.
    paired = T(RUE) or F(ALSE); 두 집단이 짝을 이루면 TRUE, 아니면 FALSE 이다.
    conf.level = 0.95 ; 입력을 안하면 기본적으로 0.95이다.

    예시)
    t.test(subset(gumjin.raw, sex==”1″)[,3], subset(gumjin.raw, sex==”2″)[,3],
            paired=FALSE, var.equal=TRUE, conf.level=0.95)

    설명)
      기본적으로 사용한 Data는 gumjin.raw 라는 항목으로 저장되어 있다.
      subset(gumjin.raw, sex==”1″)[,3]  : gumjin.raw 항목에서 sex 항목에 “1”이라고 선택되어진 자료들에서 3번째 열에 해당하는 자료를 불러오는 것을 말한다.

    결과)
            Two Sample t-test

    data:  subset(gumjin.raw, sex == “1”)[, 3] and subset(gumjin.raw, sex == “2”)[, 3]
    t = -1.2357, df = 188, p-value = 0.2181
    alternative hypothesis: true difference in means is not equal to 0
    95 percent confidence interval:
     -15.897745   3.651791
    sample estimates:
    mean of x mean of y
     180.5364  186.6593

    이렇게 결과가 나온다. 🙂

  • Excel 불러오기

    http://remantu.com/r/tutorial/load-data

    엑셀 파일인 XLS 파일을 R에서 불러오는 방법을 역시나 인터넷에서 찾았다. CSV 등으로 변환해서 불러와도 되지만 그 방법보다 좋은 점은 특정 Sheet를 지정해수 불러올 수 있다는 점이다.

    library(RODBC)
    xls = odbcConnectExcel(“파일명.xls”)
    Data 명 = sqlFetch(xls, “엑셀 Sheet 명”)

    이게 홈페이지에 있던 방법이고, 이를 응용해서 한 줄을 줄일 수도 있다.

    Data 명 = sqlFetch(odbcConnectExcel(“파일명.xls”), “엑셀 Sheet 명”)

    이 방법의 단점은 Windows 에서 실행한 R에서만 사용할 수 있다는 것이다. Mac 에서 실행한 R에서는 사용할 수가 없다. 기존의 방법대로 CVS 파일 등으로 변환을 해야 한다. ㅡㅡ;;

  • 눈빠지게..

    눈빠지게..

      
      환자 병력은 기억이 나지를 않는다. 지금 현미경 사진에 찍힌 이 부분 때문에 Toxoplasma 로 진단이 되었던 증례였던 것으로는 기억이 난다. 신경과 전공의 선생님이 이 사진을 다시 확대해서 고화질로 찍어달라고 해서 쉽게 찾을 수 있는 줄 알았다. 하지만, 찾아도 찾아도 도무지 보이지를 않았다. 400배율로 봤지만 잘 안보였다. 그렇게 하루가 지나고 이틀이 지났다. 삼일째 되던날 결국 마음을 붙잡고 400배율로 시야가 겹치도록 계속 뒤지다 보니 30분쯤 지났을 때 결국에 찾을 수 있었다. 방심 했으면 또 놓쳤을테지만 말이다.
      40배 대물렌즈로도 작게 나와서 100배 대물렌즈를 사용해서 (그러니깐 1000배) 찍었다. 좀 옛날 모델이어서 그런지 상이 조금 흐리게 나왔다. 분명히, 이걸로 증례 보고한다고 했던 것 같은데, 아직도 연락이 없는게… ‘ㅇ’
     
      다시 봐도 참 잘 찍혔다. 또 찍으라고 하면 울어버릴테다. 🙂

  • KSP 1573

    KSP 1573


      어영부영하고 있으니깐 벌써 금요일이 다가온다. 제출자 발표는 간결해야하니깐 짧게 중요한 것만 언급하는 재주를 필요로 하는데 얼마나 또 고민을 해야할지.. ㅠㅠ