레이블이 통계 프로그램 R인 게시물을 표시합니다. 모든 게시물 표시
레이블이 통계 프로그램 R인 게시물을 표시합니다. 모든 게시물 표시

2010년 6월 3일 목요일

자우루스 C1000에 통계 프로그램 R 설치하기

옛날부터 인기 있던 통계 프로그램은 SAS와 SPSS 등이지요. 요즘은 R이라는 오픈 소스 프로그램이 인기라고 합니다. 통계학 전공자 사이에서는 이미 R로 '평정'이 되었다나요. 간단한 소개는 아래 링크를 참고하세요:

http://wagnerian.textcube.com/385

그런데 R은 데스크톱 컴퓨터나 노트북 컴퓨터에서 쓰라고 만든 프로그램입니다. PDA용으로 따로 나온 것은 없는 듯해요. 휴대용 기기에서도 쓸 수 있으면 좋을 텐데 말이죠. UMPC로 분류되는 기기가 그나마 대안이 될 수 있겠고, 그 가운데 OQO가 제가 알기로 가장 작습니다. 그러나 OQO 참 비싸죠. ^^;

그런데 리눅스를 쓸 줄 알면 다른 방법이 또 있습니다. 몇몇 휴대용 기기에 우분투, 데비안 등을 설치할 수도 있더란 말이지요. 그 가운데 가장 싸고 작고 가벼운 기기는 자우루스(Zaurus) 시리즈인 듯합니다. 제가 중고로 산 자우루스 C1000에 데비안 리눅스를 설치한 얘기가 바로 아래 글에 있습니다.

http://wagnerian.textcube.com/647

이제 R을 설치해 봤습니다. 아래 링크에 자세한 설명이 나옵니다. 영어로요. ^^;

http://bibs.snu.ac.kr/R/bin/linux/debian/README.html

짧게 설명하자면, 먼저 /etc/apt/sources.list 파일을 열어서 아래와 같이 한 줄 써줍니다.

deb http://bibs.snu.ac.kr/R/bin/linux/debian lenny-cran/

인터넷 연결한 다음 터미널로 명령어를 써줍니다.

sudo apt-get update

뭔가 줄줄 받아옵니다. NOKEY 인가 어쩌고 하는 오류 메시지가 나오는데, 이상한 프로그램 함부로 설치하지 말라는 얘기이므로 무시해도 되는 듯합니다.

또 명령어를 써줍니다.

sudo apt-get install r-base

뭘 잔뜩 받아와서 설치합니다. 무선랜이라면 매우 오랜 시간이 걸립니다. -_-;;
중간에 뭘 물어보기도 하는데, 그냥 엔터, 또는 'y'라고 해주면 됩니다.

▶ R을 띄워보자

설치가 끝나면 R을 실행시켜 봅시다.

R

실행이 안 된다고요? 'r'을 소문자로 입력했죠? 대문자로 입력하지 않으면 실행되지 않습니다.

(사실은 제가 소문자로 써놓고 왜 실행이 안 될까, 한참 고민했습니다. -_-;;)

위 사진은 .csv 파일을 불러와서 간단한 상자그림(Boxplot)을 그리는 과정을 재구성해 본 것입니다. 자료는 아래 사이트에서 긁어온 주요국가 투표율 현황입니다. 한국 투표율(2008년)이 가장 낮군요.

자료 출처는 다음과 같습니다:

☞ 자작나무숲, "의무투표제가 필요하다" http://www.betulo.co.kr/1589 (2010년 6월 3일 읽음.)

이 자료에 어제 투표율을 더해서 그래프를 그려 봤습니다. 오른쪽 그림은 동유럽과 남미 등 우리보다 못하다고 믿고 싶은(…) 나라를 빼고 한국도 뺀, 그야말로 선진국들만 모은 것입니다. 선진국이 어째 투표율이 더 높군요. T-검정 결과 실제로 선진국 투표율이 전체 투표율보다 95% 신뢰구간에서 유의미하게 높다고 나왔습니다.

우리나라 투표율은 2008년 46.1%, 2010년 54.5%입니다.

▶ 주의사항

리눅스(유닉스)는 텍스트 문서에 쓰이는 줄바꿈 문자(Carriage return)가 윈도우즈와 다르므로 윈도우즈로 작업한 .csv 파일을 읽어오려면 리눅스 표준에 맞게 바꿔줘야 합니다. 자세한 내용은 아래 링크를 참고하세요:

http://goo.gl/dbct
http://en.wikipedia.org/wiki/Carriage_return

또 .csv 파일에 한글이 있으면 UTF8 형식으로 저장해야 합니다.

▶ USB 키보드 연결하기

이번에는 USB 키보드를 연결해 봅시다. USB 호스트(Host) 케이블이 있어야 합니다. 쉽게 볼 수 없는 케이블인데, 중고로 사면 같이 팔기도 하고 아니면 〈아이리버〉에서도 팝니다. 자세한 내용은 구글 검색해 보세요.

키보드를 연결하기에 앞서 터미널을 열고 이렇게 써줍니다.

modprobe ohci-hcd

이제 자우루스 ― USB 호스트케이블 ― USB 키보드를 연결해 줍니다. 글자가 써지나요? ^^;

참고로, USB 메모리를 읽으려면 해줘야 할 게 더 있습니다.

modprobe usb-storage
mkdir /mnt/usb
mount /dev/sda1 /mnt/usb

그런데 이렇게 하면 한글 파일 및 디렉터리가 깨져 보입니다. /etc/fstab을 고쳐주세요.

/dev/sda1 /mnt/usb auto noauto,umask=002,iocharset=utf8,codepage=949 0 0

2009년 4월 10일 금요일

오픈 소스 통계 프로그램 R

내가 배운 게 SAS인데 몇 년 손 놓고 있다가 (어둠의 경로로) 구해다 깔기 귀찮아서 스크립트만 만들어서 다른 사람한테 돌려달라 부탁했다. 그런데 R이라는 오픈소스 프로그램도 있다네. 이미 결과는 냈지만 다음에는 이걸 함 도전해 봐야겠다. (나중에 붙임. 써보니 대박이다. SAS는 이제 쓸 일 없을 듯.)

http://www.r-project.org/
http://lib.stat.cmu.edu/R/CRAN/index.html
R Graph Gallery: http://addictedtor.free.fr/graphiques/
R 사용자 모임 http://ruser.cyworld.com/

http://user.cs.tu-berlin.de/~ulfi/cgi-bin/r-online/r-online.cgi
R Online

http://www.statmethods.net/
sample code를 필요한 토픽별로 잘 분류해서 정리해놓은 사이트

http://www.personality-project.org/R/
예제까지 있는 사이트

http://www.remantu.com/r
한글 사이트

http://tolstoy.newcastle.edu.au/R/
Mailing list archive. 뻘질문 엄금.

도움말 사용법:
?lm
??optimization
?'topic of interest'

Tinn-R이라는 윈도우즈용 에디터도 있단다.

아래는 Shapiro-Wilk 검정 스크립트 예

#######################
era <- c(rep('B', 10), rep('A', 13) )
x <- c(39, 26, 31, 54, 42, 46, 33, 36, 48, 46, 40, 27, 12, 16, 29, 20, 30, 12, 24, 34, 20, 35, 27)
shapiro.test(x)
shapiro.test(x[era=='A'])
shapiro.test(x[era=='B'])
#######################


> shapiro.test(x)

Shapiro-Wilk normality test

data: x
W = 0.9809, p-value = 0.9209


> shapiro.test(x[era=='A'])

Shapiro-Wilk normality test

data: x[era == "A"]
W = 0.9634, p-value = 0.8051


> shapiro.test(x[era=='B'])

Shapiro-Wilk normality test

data: x[era == "B"]
W = 0.9798, p-value = 0.964

==========
붙임. 간단한 그래프 명령 예:
hist(x)
plot(density(x))
boxplot(x)


글 찾기

글 갈래