이 다음이 스크립트랑 함수파트인데... 함수 정의하고 스크립트 실행하는게 링크가 없다... 뭐 어쩌라는겨... 아 참고로 이번에도 라이브러리 하나 깔아야됩니다 install.packages("ellipse") ㅇㅋ ㄱㄱ Correlation matrix 사실 여기는 이거 하나밖에 없음...ㅋㅋㅋㅋㅋㅋ 근데 내가 이걸 어디서 본 것 같은디... (가물가물) > set.seed(955) > vvar wvar xvar yvar zvar data head(data) vvar wvar xvar yvar zvar 1 -4.252354 5.1219288 16.02193 -15.156368 -4.086904 2 1.702318 -1.3234340 15.83817 -24.063902 3.468423 3 4.323054 ..

참고로 이번꺼는 ggplot 안 데려와도 된다. 근데 라이브러리 깔긴 해야됨... install.packages("sm") install.packages("car") 네 두개 깔고 오세여. 히스토그램과 밀도 곡선 > set.seed(1) > rating=rnorm(200) > head(rating) [1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078 -0.8204684 > rating2=rnorm(200,mean=.8) > head(rating2) [1] 1.2094018 2.4888733 2.3865884 0.4690922 -1.4852355 3.2976616 다들 이쯤되면 알잖음? 히스토그램은 역사와 전통의 난수생성... > cond=factor(rep..

이거는 솔직히 8.1에 비하면 분량은 짧아요... 근데 ggplot은 불러야됨 그래프를 저⭐장 그래프 두갠가 그리긴 했음... 제주도 야채 서브셋으로... > pdf("plots.pdf") > plot(plot) 50건 이상의 경고들을 발견되었습니다 (이들 중 처음 50건을 확인하기 위해서는 warnings()를 이용하시길 바랍니다). > dev.off() pdf("파일명.pdf")만 쓰면 빈 pdf파일이 나오고 밑에 저장할 그래프를 하나씩 쓰면 페이지당 하나씩 저장된다. > pdf("plots.pdf") > plot(plot) 50건 이상의 경고들을 발견되었습니다 (이들 중 처음 50건을 확인하기 위해서는 warnings()를 이용하시길 바랍니다). > plot(plot2) 50건 이상의 경고들을 발견..

https://blog.naver.com/pokemonms/222606583751 Bradford assay 이게 뭐 하는거냐면 단백질 농도 보는 실험이다. 1. 뭐야 이거 어케해요 이게 Bradford assay용 시약이다.... blog.naver.com Bradford assay는 단백질의 무게를 확인하기 위해 진행하는 실험이다. Bradford assay용 시약을 섞고 OD595를 재면 되는데, 그러기 위해서 Standard curve가 필요하다. 일정한 무게의 단백질(BSA; Bovine serum albumin)을 용해한 다음 Bradford assay용 시약을 섞고 OD595를 측정하고, 이런 식으로 standard curve를 그린다. 이게 없으면 OD595를 재도 무게가 어느 정도인지 모..

실습용 데이터는 어지간하면 가상으로 만드는 편이지만, R에는 내장데이터가 겁나 풍부하다. 무슨 패키지 깔면 데이터 드리는 수준... 오늘은 그래서 본인 컴퓨터에 있는 R 내장 데이터 목록을 싹 털었다. 덤으로 ggplot편에 나온 데이터 출처 가르쳐드림. 들어가기 전에 소환하고 싶은 내장 데이터가 있다면 > dat=data(BJsales) 걍 이렇게 부르면 된다. > data(baseball) 경고메시지(들): In data(baseball) : 데이터셋 ‘baseball’을 찾을 수 없습니다 # 라이브러리가 필요한 건 그냥 부르면 에러뜬다 > library(plyr) > data(baseball) # 라이브러리를 부르고 부르자 라이브러리가 있어야 하는 건 라이브러리 부르고 불러야된다.. Q. 그 데이터..

그래프 제목(ggtitle()) 김후추씨의 조언대로 그래프를 만든 신입 데이터분석가. 그런데 문제가 하나 있다. "그래프에 제목을 넣고 싶은데... 어떻게 해야 할까요? " > ggplot(data=data_siseol,aes(x=채소구분소분류,y=생산량.톤.,fill=채소구분소분류))+geom_bar(stat="identity")+ggtitle("제주도 시설재배 야채 생산량") 이렇게요. > ggplot(data=data_siseol,aes(x=채소구분소분류,y=생산량.톤.,fill=채소구분소분류))+geom_bar(stat="identity")+ggtitle("제주도 시설재배 야채 생산량")+theme(plot.title=element_text(lineheight=1.5,face="bold")) 물론..

참고로 말씀드리는거지만... 분량 ㄹㅇ 역대급임... 노션으로 거의 팔만대장경 나온 듯. 데이터 관련된 얘기는 다른 글에서 다루겠습니다. 들어가기 전에 install.packages("ggplot2") 혹시나... ggplot2가 껄려있지 않다... 깔고 오세요... > library(ggplot2) 어디가요 깔았으면 불러여지. 본인은 저기다가 아예 디렉토리까지 고정으로 박고 시작했다. 막대그래프(geom_bar()) 나 저 geom 자꾸 점으로 읽어... 클났음... 아무튼! 막대그래프를 그릴 때 쓸 공공데이터는 제주도의 야채 생산 현황에 대한 공공데이터이다. 연산 채소구분대분류 채소구분소분류 면적.ha. 생산량.톤. 조수입.백만원. 1 20-21 노지채소 월동무 5056 359575 106434 2 ..

3부작으로 나눌 걸 그랬나... 라이브러리 다 깔았지들? Logistic regression(로지스틱 회귀분석) > dat=subset(mtcars,select=c(mpg,am,vs)) > dat mpg am vs Mazda RX4 21.0 1 0 Mazda RX4 Wag 21.0 1 0 Datsun 710 22.8 1 1 Hornet 4 Drive 21.4 0 1 Hornet Sportabout 18.7 0 0 Valiant 18.1 0 1 Duster 360 14.3 0 0 Merc 240D 24.4 0 1 Merc 230 22.8 0 1 Merc 280 19.2 0 1 Merc 280C 17.8 0 1 Merc 450SE 16.4 0 0 Merc 450SL 17.3 0 0 Merc 450SLC 15..

들어가기 전에 1) R의 꽃, R의 주목적이다보니 분량이 어마무시하게 길다. 네이버에서 글 분량갖고 자르진 않겠지만 그래서 상, 하편으로 나눠서 작성할 예정이다. (cookbook 소챕터 3:3으로 나눔) 2) 그리고 이새기들 자꾸 안 알려주고 뜬금없이 라이브러리 소환하세요 하는데 오늘도 깔아야 할 패키지가 있다. 하편에 나오는 것도 있고 상편에 나오는 것도 있는데 그냥 다 깔고 가자. install.packages("car") install.packages("ggplot2") install.packages("exact2x2") install.packages("irr") 3) 처음 보는 분석들이 많아서 나도 이게 뭔지 모른다. 따라서 이론적인 설명은 패스. 회귀분석과 상관계수 참고로 여기서 만드는 데이터..

내일 예고: 통계분석 들어가기 때문에 골치아파질 예정 교수님 죄송합니다 여러번 외칠 예정 이동평균 계산하기 이동평균: 전체 데이터 집합의 여러 하위 집합에 대한 일련의 평균을 만들어 데이터 요소를 분석하는 계산(솔직히 뭐 하는건지는 모르겠음) 난 sequential data라길래 파이썬처럼 시퀀스형 데이터가 있나 했더니 연속형 데이터 말하는건가봄. 전구간에서 미분 가능한가요 NA 들어가면 짤없을 예정 > set.seed(1) > x=1:300 > y=sin(x)+rnorm(300,sd=1) > y[295:300]=NA > plot(x, y, type="l", col=grey(.5)) 일단 뒤에 여백의 미를 줄 예정이다. (마른세수) > grid() 이게 모눈을 킨다고 다 이쁜 그래프가 아니그등요... 아..

들어가기 전에 아니 새기들아 깔아야 하는 라이브러리가 있으면 미리 좀 알려달라고!!! (깊은 분노) 아니 어느 레시피에서 재료설명도 없이 주저리 주저리 레시피 쓰다가 존내 당연하다는 듯 여러분 다들 집에 맨드레이크 있으시죠? 맨드레이크를 채썰어주세요. 하면서 레시피를 쓰냐!!! 집에 왜 그런게 있죠 아니 외가에서 무 받아온게 사람 모양이더라고 아무튼... 좀 개빡치긴 했지만... 라이브러리 깔고 가세요... install.packages("tidyr") install.packages("reshape2") install.packages("doBy") 테이블 가로세로 바꾸기 테이블은 보통 가로로 길거나 세로로 길거나 둘 중 하나이다. 캡처는 못했지만, 전전직장에서 일하면서 SQL로 정리해뒀던 샘플 표는 가로..

들어가기 전에 작은 시범조교를 하나(아니고 넷) 준비했음.. 다운 ㄱㄱ 각 csv파일의 내용물을 R로 불러오면 > df=read.csv('/home/koreanraichu/example.csv',sep=";") > df ID Interesred.in Class 1 kimlab0213 Python Basic 2 ahn_0526 Python Medium 3 peponi01 R Basic 4 kuda_koma R Expert 5 comma_life Java Basic 6 wheresjohn Java Medium 7 hanguk_joa Python Expert 8 sigma_00 R Basic 9 kokoatalk Java Basic (example, 구분자 세미콜론) > df2=read.csv('/home/k..

R의 데이터에는 벡터와 팩터가 있다. 그리고 숫자벡터-문자벡터-팩터간에 변환이 가능하다. 어쨌든 가능함. 팩터란 무엇인가 뮤츠씨가 좋아하는거 그건 팩트고 아무튼 벡터와 달리 팩터를 단식으로 뽑게 되면 한 가지 요소가 더 나오게 된다. 그것이 바로 '레벨'이다. > v=factor(c("A","B","C","D","E","F")) > v [1] A B C D E F Levels: A B C D E F > w=factor(c("35S Promoter","pHellsgate","pStargate","pWatergate","pHellsgate")) > w [1] 35S Promoter pHellsgate pStargate pWatergate pHellsgate Levels: 35S Promoter pHellsg..

이거 쿡복 보니까 시리즈가 개 많고... 분량이 그냥 종류별로 있습니다... 농담같지만 실화임. 그래서 세부적으로 나갈거예요... 근데 데이터프레임에 csv 불러오는 건 생각 좀 해봐야겠음. 분량이 정말 살인적입니다. 농담 아님. 아 참고로 데이터프레임 정리하기에 라이브러리가 하나 필요합니다. plyr이라고... 그거 없이 하는 방법도 있긴 한데 나중에 뭉텅이로 처리하려면 plyr 필요해요. install.packages("plyr") 설치 ㄱㄱ. sort() 벡터는 sort()로 정렬한다. 그죠 여기 벡터가 나온다는 건 데이터프레임도 있다 이겁니다. (스포일러) > v=sample(1:25) > v [1] 11 2 12 18 23 21 22 14 3 19 13 9 1 16 5 20 6 10 25 8 4..

read.csv() data data=read.csv("http://www.cookbook-r.com/Data_input_and_output/Loading_data_from_a_file/datafile.csv") > data First Last Sex Number 1 Currer Bell F 2 2 Dr. Seuss M 49 3 Student 21 웹에서 가져오는 것도 된다. read.table() data data=read.table("http://www.cookbook-r.com/Data_input_and_output/Loading_data_from_a_file/datafile.csv",strip.white=TRUE) > data V1 V2 1 First ,"Last","Sex","Number" 2 ..