안녕하세요..^^* 오늘은 R을 활용해서 행렬식을 구해볼건데요. 먼저 여인자전개 공식 없이 하나하나씩 대입해서 한번 해 볼꺼구요.. 그 다음에는 간단하게 한번 해 보도록 하겠습니다. 3곱하기 3행렬입니다. 요렇게 행렬을 표현을 해 봤습니다. 그 다음에 여인자 전개로 하나하나 공식을 넣어보고 2번째로는 det(G)라는 함수를 풀어 썼어요.. 첫번째 방식이 뻘짓이라고 생각하시는 분들 계실텐데.. 뻘짓 제가 생각해도 맞는 것 같네요. 다들 개념은 아시고 문제는 푸셔야 합니다. 개념도 모르고 막 자격증 공부하시는 식으로 외워서 집어넣으면 나중에 되면 머릿속에 들어오는게 하나도 없습니다. 그럼 오늘은 여기서 이만 마치도록 하겠습니다.^^*
안녕하세요..^^* 이번 시간에는 R을 활용을 해서 한번 평행사변형의 넓이를 구해 보도록 하는 시간을 가져 보도록 하겠습니다. 다들 평행사변형 구하는 공식은 아실텐데요. 위에 공식을 적어 놓았어요.. 벡터를 2개 이용해서 할 건데요. 넓이는 밑변 곱하기 높이인데 높이가 X_K*sin세타 이거 모르시면 다시 공부하시고 오세요.. 저는 그럼 지금부터 코딩을 해 볼께여. 벡터 a,b를 열벡터로 저렇게 만들어 줍니다. 그리고 x1과 x2를 다시 만들어 보았어요. 바로 a,b 가셔도 되는데 행렬의 첫째열을 A[,1]이면 첫번째 열의 벡터를 다 나타내라는 거 이런거 한번 해 볼려고 이렇게 했어요. 이제 스퀘어루트 즉 다 더해서 제곱근을 만들어 주는데여. t(x1)%*%x1 하면 다 제곱하는거죠. 이건 행렬의 기본을..
안녕하세요..^^* SPSS로 회귀분석을 실시하였는데 어떻게 해석 해야할지 난감하신 경우가 많으실텐데요.. 오늘은 한번 해석하는 방법을 알려 드리도록 하겠습니다. 저도 배우고 있기 때문에 교과서에 있는 나무 자료를 가져 왔습니다. 과연 부피는 지름과 높이라는 것의 변수로 회귀모형을 만족할까요?? 라고 말하는 것입니다. 즉 y=b0+b1X_i1+b2X_i2를 만족하는지 그리고 표를 만들어주면 어떻게 해석하는지 묻는 것이죠. 일단 분석에 회귀분석 들어가시고 선형을 들어가시면 됩니다. 저는 요렇게 해서 할건데요. 신뢰구간을 체크한 이유는 추정치의 95%신뢰구간을 구해보고 싶어서 입니다. 이것을 해석할려면 SST와 SSE SSR을 알아야 되는데용.. 여길 보시면 SST=SSR+SSE라고 나와 있습니다. 위에 그..
안녕하세요..^^* 혹시 회귀 분석을 하시는데 인터넷에 SPSS 회귀분석 난수생성이라던가 난수 생성을 하려고하는데 잘 안되어서 엄청나게 고민 많으신분 있으실 거라고 믿어요. 저도 진짜 몇시간동안 뱅글뱅글 돌다가 이거 겨우 했는데 다음에 까먹더라도 제 블로그에서 찾기 위해서라도 한번 포스팅을 해 볼까해요.. 어렵지 않은데 자료가 너무 없는게 흠이긴 하네요. 그럼 지금부터 한번 시작을 해 보도록 하겠습니다. 선형 회귀모델이 무엇인지는 아시죠?? ㅎㅎ 그건 다 아신다는 가정하에서 저는 spss를 사용을 해서 한번 해 보도록 하겠습니다. 저만 보는걸 잘 따라오시면 누구나 쉽게 난수를 생성 하실 수 있을 거라고 믿습니다. 자 일단 데이터 X축을 이렇게 만듭니다. 아예 모르시겠으면 저랑 똑같이 따라오세요. 1부터 ..
안녕하세요.^^* 오늘은 R프로그램을 이용해서 교적행렬을 한번 구해 볼까 합니다. 저번 시간에는 기본적인 것들을 알아 봤으니 이제 본격적으로 응용에 들어가 보는 것도 괜찮을 것 같네요. 교적행렬이란 바로 위에 있는 행렬을 뜻합니다. 여기서 X바는 바로 평균을 뜻하는거구요.. 그럼 저희가 평균도 구해야 되고 저 행렬의 원소에 평균을 빼 줘야 합니다. 그 다음에 t(X바)%*%X바까지 해야 되겠죠? 그럼 바로 코딩에 들어가 보도록 할께요. 일단 행렬을 생성해야 되겠네요. 3*3행렬을 1부터 9까지 생성을 했습니다. 그 다음에 nrow(A)라는 것을 씁니다. nrow는 바로 행의 개수를 나타내어 주는 함수입니다. 3이라고 나와 있네요. 그 다음에는 matrix(1,n,1) 3행으로 1을 만들어라는 거겠죠? 1..
안녕하세요..^^* 오늘은 행렬의 덧셈과 뺄셈 행렬식 그리고 대각 성분의 원소나 trace등을 R 프로그램을 이용해서 구해 볼건데요. C언어에 비해서 엄청 쉽습니다. 행렬을 어떻게 생성하시는지 모르는 분들은 앞포스팅에 가셔서 보시길 바랍니다. 그럼 시작하겠습니다. 아주 간단합니다. 여기서 설명하는 것은 두 행렬의 곱은 *가 아니라 %*%이라는 것만 알아두시면 됩니다. X**2는 행렬에서 X%*%X를 의미합니다. 즉 X 제곱이라는 뜻이죠. 저 자리에 3을 넣으면 3제곱이 되겠죠? 아래꺼랑 짤렸네요.. 바로 밑에 있는 것이 X제곱행렬입니다. determinant, transpose, diagonal element, slove는 역행렬을 묻는 것이구요. 마지막 거는 대각행렬의 성분들을 다 더하니까 trace가..
안녕하세요..^^* 요즘 시험기간이라 제 전공공부를 할때가 많습니다. 이렇게 공부하기 귀찮을때면 제가 하는 방법이 있죠. 바로 블로그에 쓰는 겁니다. 이렇게 하면 손으로 쓰는 것 보다 시간도 빠를 뿐만 아니라 남들을 이해시켜야 하기 때문에 저 또한 엄청나게 공부가 됩니다. 오늘은 간단하게 R프로그램을 이용해서 어떻게 행렬을 만드는지 알아 보도록 하겠습니다. R프로그램을 깔고 실행을 시킵니다. 빈 페이지를 열어서 오른쪽에서 코딩하신다음에 복사 붙여넣게 해서 왼쪽 스패이스에 하시면 됩니다. 빈페이지에서 제가 행렬을 몇개 적어 봣습니다. 아주 간단합니다. 행렬을 만드는 방법은 A
저는 통계를 부전공으로 하고 있습니다. 그래서 통계 프로그램을 다루는 기회가 꽤 생기죠. 그래서 저는 일단 통계 패키지를 이용해서 처음으로 프로그램을 짜게 되었습니다. 물론 sas를 이용해서 이죠. 바로 단일품목 판매 프로그램인데요. 그것이 뭔지 지금부터 차근차근 짜 보도록 하겠습니다. 상황은 이렇습니다. 신문팔이 소년이 신문 한부를 150원에 사서 300원에 팔고 못 판 것은 1 부당 100원씩 환불받습니다. 그런데 매일 몇 부를 팔아야 최대 이익을 얻을지 고민합니다. 그래서 100일 동안 매일의 판매 부수를 조사했습니다. 그것을 확률분포표로 나타내었죠. 어떻게 조사했냐 하면 1000일 동안 신문을 팔아 봤더니 판매부수에 따른 확률을 구했습니다. 여기서 판매부수 10 부수가 0.04라는 말은 1000일..