회귀분석... 아직도 난 이게 뭔지 모르겠어...회귀분서억관찰된 연속형 변수들에 대해 두 변수 사이의 모형을 구한 뒤 적합도를 측정해 내는 분석 방법이라고 하는데, 그래서 이게 뭔데요? 하나(단순선형-y = β₀ + β₁x + ε) 혹은 그 이상(다중선형-단순선형에서 x가 오지게 많아짐)의 독립변수들이 종속변수에 미치는 영향을 추정할 수 있는 통계기법 이다. 일차식으로 나오는것도 있고, 로지스틱이라고 해서 0 또는 1만 있는 거 할 때도 있는데 뒤에껀 무슨 변환식 보면 갑자기 자연로그가 튀어나오더니 시그모이드 등장함... 하다보면 이항 다항 튀어나오고 난리납니다. 다중공선성회귀분석을 망칠수도 있는 요인이다. 이게 뭐냐고? 여러분들 멘델의 유전법칙 아십니까? 우열의 법칙이랑 독립의 법칙이라고 있는데, 독립..
이걸 왜 알아야 하냐면 결과 주고 옳은 것은? 옳지 않은 것은? 해석하라고 나옵니다. 연관, 회귀(주로 다항), 또 뭐 있었는데 아무튼... 그래서 분석도 해볼 수 있으면 직접 찍먹정도는 해보는게 좋다는거임. 연관분석이 뭔데요?내가 겜덕이니까 겜 위주로 비유를 좀 해보겠음. 요즘 스위치 2 나왔잖아요? 연관분석에는 크게 장바구니 분석, 서열분석 두 개가 있는데, 장바구니 분석은 스위치 2를 사는 고객들이 뭘 같이 사는가(액정필름, 파우치, 냥발(스틱 커버) 이런거)이고 서열분석은 스위치 2를 구매한 고객들이 이후에 뭘 사는가를 분석하는거다. 게임기로는 잘 안 와닿는 분들을 위해 한가지 더 예를 들어보자면, 내가 지금쓰고 있는 폰이 아이폰 15 프맥인데 얘를 살 때 액정필름이랑 케이스를 같이 샀었다. (아..
일단 이항분포가 뭐냐… 특정 확률(p)을 가진 베르누이 시행을 n번 독립적으로 반복했을 때, 성공하는 횟수(X)에 대한 이산 확률 분포라고 한다. Pass or Fail 뭐 이런건데, 여기서 중요한 건 결과가 두개라고 확률까지 반반이라는 건 아니라는 얘기. 그것은 하등 근거없는 편견이다. 근데 확률이 너무 낮아서 이항분포 때릴 수 있나 모르겠음.확률이 얼마길래?우리가 여기서 해볼 건 1/4096(no빛부 야생), 1/1365(빛부), 1/512(빛부+국제교배)이다. R에서 이항분포 때려보기> x_val=0:10> y_val1=dbinom(x_val,100,1/4096)일단 위는 x가 0부터 10까지라는 얘기이고, 아래가 이항분포 그 뭐시기다. 저걸 쉽게 풀어주자면 1) 1/4096의 확률을 가지고 있는 ..
이 다음이 스크립트랑 함수파트인데... 함수 정의하고 스크립트 실행하는게 링크가 없다... 뭐 어쩌라는겨... 아 참고로 이번에도 라이브러리 하나 깔아야됩니다 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) 처음 보는 분석들이 많아서 나도 이게 뭔지 모른다. 따라서 이론적인 설명은 패스. 회귀분석과 상관계수 참고로 여기서 만드는 데이터..