디시인사이드 갤러리

갤러리 이슈박스, 최근방문 갤러리

갤러리 본문 영역

재획하면서 공부하기 #1

ㅇㅇ갤로그로 이동합니다. 2024.08.30 20:09:04
조회 55 추천 0 댓글 0
														

regression

classification


X : INPUT

x1, x2, x3 ... : feature

Y : OUTPUT


Y = F(X)


ex) Y = w1x1 + w2x2 + e

w1, w2 : weight, parameter

e : error, loss function

Σei^2 : cost function

cost function을 최소화 시키는 w1, w2 찾는 것이 F를 찾는것


모델 결정 -> 파라미터 추정


모델에는 다중선형회귀모델, 로지스틱회귀모델, 의사결정나무모델, 인공신경망모델, ...


각 모델별로 파라미터 추정 알고리즘이 있다.


다중선형회귀모델의 경우 least square estimation 사용,

로지스틱회귀모델의 경우 conjugate gradient algorithm 사용,

인공신경망모델의 경우 back-propagation algorithm 사용




선형회귀모델

Y = w0 + w1x1 + ... + wnxn

가정) e : error는 정규분포를 따른다. 

파라미터 추정) least square estimation

Σei^2 을 최소화 하는 w0, w1, ...,wn 을 찾는것



코딩 직접 해보기 데이터갯수에 따라 plot 하기

plot된 데이터를 바탕으로 직선 구하기, 직선 플롯하기, cost func 구하기


파라미터에 대한 구간추정

점 추정량 값 플마 상수 * 점추정량에 대한 표준편차

여기서 상수는 t_a/2,n-2 : 유의수준 1-a 하에서 자유도가 n-2인 t분포의 값


결정 계수 R^2 = SSR/SST = 1 - SSE/SST

SSR/SST = 1 : 확정적인 관계

SSR/SST = 0 : X변수가 Y변수를 설명하는데 전혀 도움이 되지 않는다.

유의하지 않은 변수가 추가되어도 항상 증가함

-> Adjusted R^2 사용



## y=2x+50+error 를 따르는 n개의 데이터셋 구성, error~N(0,1)
## x는 100이하의 자연수
## 실험목표: n=4, n=8, n=16 인 케이스에 대해 데이터셋 plot 해보기.
## 각 케이스에 대해 데이터셋을 바탕으로 선형회귀모델 y=w0+w1x 구해보기.
## 각 케이스에 대해 cost function  값 계산하기


import random
import matplotlib.pyplot as plt
n = 4
X_input_4 = []
Y = []
for i in range(n):
    X_input_4.append(random.randint(1, 100))
    error = random.gauss(0, 1)
    Y.append(2*X_input_4[i]+50+error)
print(X_input_4)
print(Y)
plt.scatter(X_input_4, Y)
plt.xlabel('X_input_4')
plt.ylabel('Y')
plt.title('Scatter plot of X_input_4 and Y')
plt.show()


import numpy as np
from sklearn.linear_model import LinearRegression

# Convert X_input_4 to a NumPy array and reshape it
X = np.array(X_input_4).reshape(-1, 1)

# Create a linear regression model
model = LinearRegression()

# Fit the model to the data
model.fit(X, Y)

# Get the coefficients
w0 = model.intercept_
w1 = model.coef_[0]

# Print the equation of the line
print("y = {:.2f} + {:.2f}x".format(w0, w1))

# Generate points for the line
X_line = np.linspace(min(X_input_4), max(X_input_4), 100)
Y_line = w0 + w1 * X_line

# Plot the data points and the best fit line
plt.scatter(X_input_4, Y)
plt.plot(X_line, Y_line, color='red')
plt.xlabel('X_input_4')
plt.ylabel('Y')
plt.title('Linear Regression')
plt.show()


추천 비추천

0

고정닉 0

2

댓글 영역

전체 댓글 0
등록순정렬 기준선택
본문 보기

하단 갤러리 리스트 영역

왼쪽 컨텐츠 영역

갤러리 리스트 영역

갤러리 리스트
번호 제목 글쓴이 작성일 조회 추천
설문 자식 사교육에 돈 엄청 쓸 것 같은 스타는? 운영자 24/09/09 - -
AD <소드 오브 콘발라리아> 첫 대규모 업데이트 실시! 운영자 24/09/13 - -
8410259 뉴비 오늘 매할일 다함 이제 집감요 ㅋㅋ 메갤러(121.144) 16:50 20 0
8410258 롤대남티비화면에비추면방통위에서경고들어올듯…. 나옹이갤로그로 이동합니다. 16:50 22 0
8410255 근우애미창년새끼 칼로 찔러죽이고싶네 [2] ㅇㅇ(175.114) 16:49 45 3
8410254 지든노 [4] 망고초밥갤로그로 이동합니다. 16:48 26 0
8410253 지구온난화로 난리칠거면 미국,중국,인도한테해야댐 [2] 구르몽갤로그로 이동합니다. 16:48 21 0
8410251 똥마렵다 메갤러(121.144) 16:48 9 0
8410250 예슬이 롤갤에서 노는중인데 좀 데려가라 [1] ㅇㅇ(106.102) 16:48 23 1
8410249 ㄴ 살다살다 이런 경우없는 새끼는 처음보면 개추 ㅇㅇ(118.235) 16:48 17 0
8410248 더스크 14분에 70퍼 까는데 [5] 메갤러(58.79) 16:47 52 0
8410247 한번만 도와주세요 메갤러(58.235) 16:47 13 0
8410245 한남들 얼마나 빻앗으면 여자들만 치어풀보여주냐 [2] ㅇㅇ(61.39) 16:47 51 0
8410244 이 씹새들 리린 기록말살형 들어갔냐? ㅇㅇ(182.221) 16:47 11 0
8410242 백인도 한녀랑 결혼하면 이렇게 되는구나 시발 ㅋㅋㅋㅋ [2] ㅇㅇ(211.36) 16:46 48 0
8410241 감귤66작 돈없어서 못한거아님??? 메갤러(118.235) 16:46 26 0
8410240 엔임나연버 이 미친년 메갤러(211.235) 16:46 53 2
8410239 구마유시는 왜 진 직스 메타에도 페바뎊보다 후달림? [1] ㅇㅇ갤로그로 이동합니다. 16:46 20 0
8410238 9월중순인데도 33도라니 [12] 쯔위안갤로그로 이동합니다. 16:46 53 0
8410235 [공 짜] 배달의민족 10,000원짜리 쓸사람 써 내일까지임! [1] ㅇㅇ(223.39) 16:45 22 0
8410233 앱솔랩스 코인 교환 질문 메갤러(39.117) 16:45 18 0
8410232 리선족 서버이름 뭘로 바뀔지가 ㅈㄴ 기대되네 [2] ㅇㅇ(118.235) 16:45 28 0
8410231 부엉이를 좋아하는 여고딩.....JPG ㅇㅇ(210.98) 16:45 42 0
8410230 뭐야 vip 부스터 시세왜이ㅏ래 늦게팔면 비싼거아님? [3] 메갤러(121.144) 16:44 46 0
8410229 요즘도 친척들 만나고 그럼???????????????????? ㅇㅇ갤로그로 이동합니다. 16:44 19 0
8410228 지금 마약이벤 없는거지? ㅇㅇ(122.34) 16:44 14 0
8410227 비왓는데 어케ㅜ더 더워지냐 [1] 얏바리갤로그로 이동합니다. 16:44 16 0
8410226 리갤 쟤네는 언제까지 저럴 생각인거지? ㅇㅇ갤로그로 이동합니다. 16:44 24 0
8410225 저번주는좀 시원햇던거같은데 ㅅㅂ 정서기갤로그로 이동합니다. 16:44 12 0
8410224 출간도서가 환장하는거 ㅇㅇ갤로그로 이동합니다. 16:44 16 0
8410222 원더베리vs플레티넘애플 [3] 10분갤로그로 이동합니다. 16:43 31 0
8410221 영화관 걸어서 10분거리라 좋다 얏바리갤로그로 이동합니다. 16:43 21 0
8410220 고브 몬컬 확률업 체감 좆되노 ㅇㅇ(39.117) 16:42 16 0
8410219 가을인데시발날씨좀풀려라 10도정도로 [4] 정서기갤로그로 이동합니다. 16:42 33 0
8410218 구마유시 직스는 ㄹㅇ신명호네 [2] 윌비백갤로그로 이동합니다. 16:42 47 0
8410217 메이플하고나서 전자제품사는게 꺼려짐 ㅇㅇ(118.235) 16:42 37 0
8410216 오랫만에 다시 해보려는데 자본 얼마나 있어야됨? ㅇㅇ(182.217) 16:42 16 0
8410215 페이커는 증명해낼거임 메갤러(221.164) 16:42 15 0
8410212 사냥 쉬고싶어도 자리가 없으니까 쉴수가업네 ㅋㅋ 메갤러(118.38) 16:42 12 0
8410211 ㄴ 한근우면 탈갤함 메갤러(106.101) 16:41 11 0
8410210 다큰 어른이 란도셀매고 유혹을 하며 닝얀갤로그로 이동합니다. 16:41 30 0
8410209 스카 메결할 여캐본캐인사람 있음? 메갤러(118.235) 16:41 12 0
8410208 딱 1년전에 ㅇㅇ(118.235) 16:41 11 0
8410207 요즘 1000억 스타트하면 어디까지함? [5] 산타갤로그로 이동합니다. 16:40 44 0
8410206 ㄴ 이새끼 한근우급 병신 메갤러(106.101) 16:40 10 0
8410205 형들 새매물찾음 22억 이륙 고고? [1] 메갤러(121.144) 16:40 38 0
8410204 이거 70억인데 사면 병신호구임? [6] 메갤러(121.160) 16:40 60 2
8410203 스햐등장 스톤멍키ㅇㄷ(119.64) 16:40 14 0
8410202 템 일괄로 좀 못사나 메갤러(106.101) 16:40 11 0
8410201 썬콜오리진빙결스택영향받음? 메갤러(222.104) 16:40 9 0
8410200 날씨진짜병신같네 ㅇㅇ(115.138) 16:39 17 0
8410199 재벌집 막내 가슴 ㅇㅇ(117.111) 16:39 50 0
갤러리 내부 검색
제목+내용게시물 정렬 옵션

오른쪽 컨텐츠 영역

실시간 베스트

1/8

뉴스

디시미디어

디시이슈

1/2