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)


글 찾기

글 갈래