Coding/R
-
이 다음이 스크립트랑 함수파트인데... 함수 정의하고 스크립트 실행하는게 링크가 없다... 뭐 어쩌라는겨... 아 참고로 이번에도 라이브러리 하나 깔아야됩니다 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 ..
R 배워보기-8.4. Other interesting graphs이 다음이 스크립트랑 함수파트인데... 함수 정의하고 스크립트 실행하는게 링크가 없다... 뭐 어쩌라는겨... 아 참고로 이번에도 라이브러리 하나 깔아야됩니다 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 ..
2022.08.23 -
참고로 이번꺼는 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..
R 배워보기-8.3. Basic graphs with standard graphics참고로 이번꺼는 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..
2022.08.23 -
이거는 솔직히 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건 이상의 경고들을 발견..
R 배워보기-8.2. Miscellaneous이거는 솔직히 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건 이상의 경고들을 발견..
2022.08.23 -
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로 standard curve 그리기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를 재도 무게가 어느 정도인지 모..
2022.08.22 -
실습용 데이터는 어지간하면 가상으로 만드는 편이지만, R에는 내장데이터가 겁나 풍부하다. 무슨 패키지 깔면 데이터 드리는 수준... 오늘은 그래서 본인 컴퓨터에 있는 R 내장 데이터 목록을 싹 털었다. 덤으로 ggplot편에 나온 데이터 출처 가르쳐드림. 들어가기 전에 소환하고 싶은 내장 데이터가 있다면 > dat=data(BJsales) 걍 이렇게 부르면 된다. > data(baseball) 경고메시지(들): In data(baseball) : 데이터셋 ‘baseball’을 찾을 수 없습니다 # 라이브러리가 필요한 건 그냥 부르면 에러뜬다 > library(plyr) > data(baseball) # 라이브러리를 부르고 부르자 라이브러리가 있어야 하는 건 라이브러리 부르고 불러야된다.. Q. 그 데이터..
R의 내장 데이터 (부제: 공공데이터 어떻게 받아요?)실습용 데이터는 어지간하면 가상으로 만드는 편이지만, R에는 내장데이터가 겁나 풍부하다. 무슨 패키지 깔면 데이터 드리는 수준... 오늘은 그래서 본인 컴퓨터에 있는 R 내장 데이터 목록을 싹 털었다. 덤으로 ggplot편에 나온 데이터 출처 가르쳐드림. 들어가기 전에 소환하고 싶은 내장 데이터가 있다면 > dat=data(BJsales) 걍 이렇게 부르면 된다. > data(baseball) 경고메시지(들): In data(baseball) : 데이터셋 ‘baseball’을 찾을 수 없습니다 # 라이브러리가 필요한 건 그냥 부르면 에러뜬다 > library(plyr) > data(baseball) # 라이브러리를 부르고 부르자 라이브러리가 있어야 하는 건 라이브러리 부르고 불러야된다.. Q. 그 데이터..
2022.08.22 -
그래프 제목(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")) 물론..
R 배워보기-8.1. ggplot2로 그래프 그리기 (하)그래프 제목(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")) 물론..
2022.08.22 -
참고로 말씀드리는거지만... 분량 ㄹㅇ 역대급임... 노션으로 거의 팔만대장경 나온 듯. 데이터 관련된 얘기는 다른 글에서 다루겠습니다. 들어가기 전에 install.packages("ggplot2") 혹시나... ggplot2가 껄려있지 않다... 깔고 오세요... > library(ggplot2) 어디가요 깔았으면 불러여지. 본인은 저기다가 아예 디렉토리까지 고정으로 박고 시작했다. 막대그래프(geom_bar()) 나 저 geom 자꾸 점으로 읽어... 클났음... 아무튼! 막대그래프를 그릴 때 쓸 공공데이터는 제주도의 야채 생산 현황에 대한 공공데이터이다. 연산 채소구분대분류 채소구분소분류 면적.ha. 생산량.톤. 조수입.백만원. 1 20-21 노지채소 월동무 5056 359575 106434 2 ..
R 배워보기-8.1. ggplot2로 그래프 그리기 (상)참고로 말씀드리는거지만... 분량 ㄹㅇ 역대급임... 노션으로 거의 팔만대장경 나온 듯. 데이터 관련된 얘기는 다른 글에서 다루겠습니다. 들어가기 전에 install.packages("ggplot2") 혹시나... ggplot2가 껄려있지 않다... 깔고 오세요... > library(ggplot2) 어디가요 깔았으면 불러여지. 본인은 저기다가 아예 디렉토리까지 고정으로 박고 시작했다. 막대그래프(geom_bar()) 나 저 geom 자꾸 점으로 읽어... 클났음... 아무튼! 막대그래프를 그릴 때 쓸 공공데이터는 제주도의 야채 생산 현황에 대한 공공데이터이다. 연산 채소구분대분류 채소구분소분류 면적.ha. 생산량.톤. 조수입.백만원. 1 20-21 노지채소 월동무 5056 359575 106434 2 ..
2022.08.22 -
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..
R 배워보기-7. Statistical analysis (하)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..
2022.08.21 -
들어가기 전에 1) R의 꽃, R의 주목적이다보니 분량이 어마무시하게 길다. 네이버에서 글 분량갖고 자르진 않겠지만 그래서 상, 하편으로 나눠서 작성할 예정이다. (cookbook 소챕터 3:3으로 나눔) 2) 그리고 이새기들 자꾸 안 알려주고 뜬금없이 라이브러리 소환하세요 하는데 오늘도 깔아야 할 패키지가 있다. 하편에 나오는 것도 있고 상편에 나오는 것도 있는데 그냥 다 깔고 가자. install.packages("car") install.packages("ggplot2") install.packages("exact2x2") install.packages("irr") 3) 처음 보는 분석들이 많아서 나도 이게 뭔지 모른다. 따라서 이론적인 설명은 패스. 회귀분석과 상관계수 참고로 여기서 만드는 데이터..
R 배워보기-7. Statistical analysis (상)들어가기 전에 1) R의 꽃, R의 주목적이다보니 분량이 어마무시하게 길다. 네이버에서 글 분량갖고 자르진 않겠지만 그래서 상, 하편으로 나눠서 작성할 예정이다. (cookbook 소챕터 3:3으로 나눔) 2) 그리고 이새기들 자꾸 안 알려주고 뜬금없이 라이브러리 소환하세요 하는데 오늘도 깔아야 할 패키지가 있다. 하편에 나오는 것도 있고 상편에 나오는 것도 있는데 그냥 다 깔고 가자. install.packages("car") install.packages("ggplot2") install.packages("exact2x2") install.packages("irr") 3) 처음 보는 분석들이 많아서 나도 이게 뭔지 모른다. 따라서 이론적인 설명은 패스. 회귀분석과 상관계수 참고로 여기서 만드는 데이터..
2022.08.21 -
내일 예고: 통계분석 들어가기 때문에 골치아파질 예정 교수님 죄송합니다 여러번 외칠 예정 이동평균 계산하기 이동평균: 전체 데이터 집합의 여러 하위 집합에 대한 일련의 평균을 만들어 데이터 요소를 분석하는 계산(솔직히 뭐 하는건지는 모르겠음) 난 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() 이게 모눈을 킨다고 다 이쁜 그래프가 아니그등요... 아..
R 배워보기- 6.5. Manipulating data-Sequential data내일 예고: 통계분석 들어가기 때문에 골치아파질 예정 교수님 죄송합니다 여러번 외칠 예정 이동평균 계산하기 이동평균: 전체 데이터 집합의 여러 하위 집합에 대한 일련의 평균을 만들어 데이터 요소를 분석하는 계산(솔직히 뭐 하는건지는 모르겠음) 난 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() 이게 모눈을 킨다고 다 이쁜 그래프가 아니그등요... 아..
2022.08.20 -
들어가기 전에 아니 새기들아 깔아야 하는 라이브러리가 있으면 미리 좀 알려달라고!!! (깊은 분노) 아니 어느 레시피에서 재료설명도 없이 주저리 주저리 레시피 쓰다가 존내 당연하다는 듯 여러분 다들 집에 맨드레이크 있으시죠? 맨드레이크를 채썰어주세요. 하면서 레시피를 쓰냐!!! 집에 왜 그런게 있죠 아니 외가에서 무 받아온게 사람 모양이더라고 아무튼... 좀 개빡치긴 했지만... 라이브러리 깔고 가세요... install.packages("tidyr") install.packages("reshape2") install.packages("doBy") 테이블 가로세로 바꾸기 테이블은 보통 가로로 길거나 세로로 길거나 둘 중 하나이다. 캡처는 못했지만, 전전직장에서 일하면서 SQL로 정리해뒀던 샘플 표는 가로..
R 배워보기- 6.4. Manipulating data-Restructing data들어가기 전에 아니 새기들아 깔아야 하는 라이브러리가 있으면 미리 좀 알려달라고!!! (깊은 분노) 아니 어느 레시피에서 재료설명도 없이 주저리 주저리 레시피 쓰다가 존내 당연하다는 듯 여러분 다들 집에 맨드레이크 있으시죠? 맨드레이크를 채썰어주세요. 하면서 레시피를 쓰냐!!! 집에 왜 그런게 있죠 아니 외가에서 무 받아온게 사람 모양이더라고 아무튼... 좀 개빡치긴 했지만... 라이브러리 깔고 가세요... install.packages("tidyr") install.packages("reshape2") install.packages("doBy") 테이블 가로세로 바꾸기 테이블은 보통 가로로 길거나 세로로 길거나 둘 중 하나이다. 캡처는 못했지만, 전전직장에서 일하면서 SQL로 정리해뒀던 샘플 표는 가로..
2022.08.20 -
들어가기 전에 작은 시범조교를 하나(아니고 넷) 준비했음.. 다운 ㄱㄱ 각 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 배워보기- 6.3. Manipulating data-Data Frames들어가기 전에 작은 시범조교를 하나(아니고 넷) 준비했음.. 다운 ㄱㄱ 각 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..
2022.08.20 -
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..
R 배워보기- 6.2. Manipulating data-FactorsR의 데이터에는 벡터와 팩터가 있다. 그리고 숫자벡터-문자벡터-팩터간에 변환이 가능하다. 어쨌든 가능함. 팩터란 무엇인가 뮤츠씨가 좋아하는거 그건 팩트고 아무튼 벡터와 달리 팩터를 단식으로 뽑게 되면 한 가지 요소가 더 나오게 된다. 그것이 바로 '레벨'이다. > 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..
2022.08.20 -
이거 쿡복 보니까 시리즈가 개 많고... 분량이 그냥 종류별로 있습니다... 농담같지만 실화임. 그래서 세부적으로 나갈거예요... 근데 데이터프레임에 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..
R 배워보기- 6.1. Manipulating data-General이거 쿡복 보니까 시리즈가 개 많고... 분량이 그냥 종류별로 있습니다... 농담같지만 실화임. 그래서 세부적으로 나갈거예요... 근데 데이터프레임에 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..
2022.08.20 -
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 ..
R 배워보기-5. 데이터 불러오고 쓰기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 ..
2022.08.20 -
하는김에 어? 이거 되나? 해서 찾아봄. 미분 > f D(f,"x") a * (5 * x^4) + b * (2 * x) 미분은 이런 식으로 한다. 근데 계산은 안 해준다. 믿고보는 알파신이 짱이다 > f D(f,"x") exp(x) y=e^x를 미분하면 이게 정상이다. 으아아 미분했는데 왜 그대로야 정상입니다 손님 편미분 마려워요? > f D(f,'x') 2 * x * y^2 + 2 * y # x로 미분해주세요 > D(f,'y') x^2 * (2 * y) + 2 * x + 1 # y로 미분해주세요 저 D에서 기호를 바꿔주면 알아서 편미분이 된단다. 위 코드는 각각 x로 한 번, y로 한 번 미분한 결과. (그러니까 어쨌든 한 번 미분했다) > D(D(f,'x'),'y') 2 * x * (2 * y) + ..
번외편-R로 미적분 하기하는김에 어? 이거 되나? 해서 찾아봄. 미분 > f D(f,"x") a * (5 * x^4) + b * (2 * x) 미분은 이런 식으로 한다. 근데 계산은 안 해준다. 믿고보는 알파신이 짱이다 > f D(f,"x") exp(x) y=e^x를 미분하면 이게 정상이다. 으아아 미분했는데 왜 그대로야 정상입니다 손님 편미분 마려워요? > f D(f,'x') 2 * x * y^2 + 2 * y # x로 미분해주세요 > D(f,'y') x^2 * (2 * y) + 2 * x + 1 # y로 미분해주세요 저 D에서 기호를 바꿔주면 알아서 편미분이 된단다. 위 코드는 각각 x로 한 번, y로 한 번 미분한 결과. (그러니까 어쨌든 한 번 미분했다) > D(D(f,'x'),'y') 2 * x * (2 * y) + ..
2022.08.20 -
공식 번외편도 따로 나갑니다. as.formula() y~x y ~ x 이놈은 문자고 as.formula("y~x") y ~ x 이놈은 공식이다. 뭔 차이인지는 모르겠다. > measurevar="y" > groupvars=c("x1","x2","x3") paste(measurevar,paste(groupvars,collapse="+"),sep="=") [1] "y=x1+x2+x3" 이놈도 문자인데 as.formula(paste(measurevar, paste(groupvars, collapse=" + "), sep=" ~ ")) y ~ x1 + x2 + x3 이놈은 공식이다. 역시 뭔 차이인지 모른다. 공식의 구조 t=y~x1+x2 이렇게 공식을 만들어서 확인해보면(아, 참고로 f=as.formula(..
R 배워보기-4. 공식공식 번외편도 따로 나갑니다. as.formula() y~x y ~ x 이놈은 문자고 as.formula("y~x") y ~ x 이놈은 공식이다. 뭔 차이인지는 모르겠다. > measurevar="y" > groupvars=c("x1","x2","x3") paste(measurevar,paste(groupvars,collapse="+"),sep="=") [1] "y=x1+x2+x3" 이놈도 문자인데 as.formula(paste(measurevar, paste(groupvars, collapse=" + "), sep=" ~ ")) y ~ x1 + x2 + x3 이놈은 공식이다. 역시 뭔 차이인지 모른다. 공식의 구조 t=y~x1+x2 이렇게 공식을 만들어서 확인해보면(아, 참고로 f=as.formula(..
2022.08.20 -
쿡북 분량은 일단 짧다. 다음것도 짧다. grep() v=c("피카츄","피츄","라이츄","에몽가","따라큐") w=c("Alticuno","moltres","zapdos","lugia","Ho-oh") 참고로 오늘의 시범조교다. grep("따라큐",v) [1] 5 grep은 문자를 찾아주는 함수. grep("alticuno",w) integer(0) ...인데 대소문자를 가린다. grep("alticuno",w,ignore.case=TRUE) [1] 1 그래서 ignore.case=TRUE를 줘야 대소문자 상관 없이 찾아준다. ^와 $ 각각 문자의 시작과 끝을 지정해서 찾아주는 것. grep("^피",v) [1] 1 2 # 시작 문자 지정하기 grep("츄$",v) [1] 1 2 3 # 끝 문자 지정..
R 배워보기-3. 문자열쿡북 분량은 일단 짧다. 다음것도 짧다. grep() v=c("피카츄","피츄","라이츄","에몽가","따라큐") w=c("Alticuno","moltres","zapdos","lugia","Ho-oh") 참고로 오늘의 시범조교다. grep("따라큐",v) [1] 5 grep은 문자를 찾아주는 함수. grep("alticuno",w) integer(0) ...인데 대소문자를 가린다. grep("alticuno",w,ignore.case=TRUE) [1] 1 그래서 ignore.case=TRUE를 줘야 대소문자 상관 없이 찾아준다. ^와 $ 각각 문자의 시작과 끝을 지정해서 찾아주는 것. grep("^피",v) [1] 1 2 # 시작 문자 지정하기 grep("츄$",v) [1] 1 2 3 # 끝 문자 지정..
2022.08.20 -
파이썬에서는 똑같은 걸 해보려면 모듈을 불러야 하는데 얘는 모듈 X까! 걍 해! 같은 느낌... 이 다음편 문자열인데 분량 개짧습니다(스포일러) 난수 만들기 runif(1) [1] 0.7232427 runif()를 쓰면 0부터 1까지 중 아무 숫자나 하나 출력한다. runif(4) [1] 0.8477728 0.4359127 0.4291748 0.4625472 저 괄호 안에 숫자는 범위가 아니고 개수 지정하는거다. runif(5,min=0,max=100) [1] 87.75978 48.65714 87.73802 18.87537 75.90590 범위는 이런 식으로 지정한다. floor(runif(6,min=0,max=7)) [1] 1 4 6 4 4 4 소수점 떼뿌라!!! floor는 밑에 설명할건데 그거 보면..
R 배워보기-2. Numbers파이썬에서는 똑같은 걸 해보려면 모듈을 불러야 하는데 얘는 모듈 X까! 걍 해! 같은 느낌... 이 다음편 문자열인데 분량 개짧습니다(스포일러) 난수 만들기 runif(1) [1] 0.7232427 runif()를 쓰면 0부터 1까지 중 아무 숫자나 하나 출력한다. runif(4) [1] 0.8477728 0.4359127 0.4291748 0.4625472 저 괄호 안에 숫자는 범위가 아니고 개수 지정하는거다. runif(5,min=0,max=100) [1] 87.75978 48.65714 87.73802 18.87537 75.90590 범위는 이런 식으로 지정한다. floor(runif(6,min=0,max=7)) [1] 1 4 6 4 4 4 소수점 떼뿌라!!! floor는 밑에 설명할건데 그거 보면..
2022.08.20 -
이거 미디움에도 올려야하나 좀 고민인게 일단 쿡북 분량이 생각보다 좀 되고 미디움에는 코드블럭이 없음... 네? 설치요? 구글가서 R 설치하는법 찾아보세요. R studio도 같이 깔아야됨. Python과의 차이점 1) Python과 달리 R은 인덱스 번호가 1부터 시작이다 2) 음수 인덱싱이 Python과 달리 빼고 출력하라는 의미이다(파이썬은 맨 뒤에꺼 달라는 얘기) 3) Python pandas는 기본적으로 결측값을 빼고 계산하지만 R은 결측값을 넣은 상태에서 계산한다. 물론 둘 다 반대로 설정하는 옵션이 존재. 4) R은 결측값을 서브셋 이용해서 대체할 수는 있지만 dropna()같은 기능은 없는 듯 하다. Data Indexing v = c(1,2,4,8,16,32) 일단 벡터는 이런 식이다. ..
R 배워보기-1. Basics이거 미디움에도 올려야하나 좀 고민인게 일단 쿡북 분량이 생각보다 좀 되고 미디움에는 코드블럭이 없음... 네? 설치요? 구글가서 R 설치하는법 찾아보세요. R studio도 같이 깔아야됨. Python과의 차이점 1) Python과 달리 R은 인덱스 번호가 1부터 시작이다 2) 음수 인덱싱이 Python과 달리 빼고 출력하라는 의미이다(파이썬은 맨 뒤에꺼 달라는 얘기) 3) Python pandas는 기본적으로 결측값을 빼고 계산하지만 R은 결측값을 넣은 상태에서 계산한다. 물론 둘 다 반대로 설정하는 옵션이 존재. 4) R은 결측값을 서브셋 이용해서 대체할 수는 있지만 dropna()같은 기능은 없는 듯 하다. Data Indexing v = c(1,2,4,8,16,32) 일단 벡터는 이런 식이다. ..
2022.08.20