Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 16 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
16
Dung lượng
337,49 KB
Nội dung
TIỂU LUẬN HỌC PHẦN : PP PHÂN TÍCH DỮ LIỆU NGÀNH DỆT MAY Xử lý phân tích liệu ảnh hưởng thông số dán tới độ bền liên kết dán dựng Yêu cầu 1: Mô tả biến với đặc trưng thống kê biến: số trung bình, trung vị, độ lệch chuẩn, phương sai, độ biến thiên, khoảng tin cậy 95% - Đọc liệu dandung.csv gọi object dandung > dandung=read.csv("D:/dan dung.csv",header=T) > attach(dandung) > names(dandung) [1] "nhietdodan" "aplucdan" "thoigiandan" [4] "dobenxelienketdan" "dogian" Mô tả đặc trưng thống kê biến: "nhietdodan" > summary(nhietdodan) >sd=sd(nhietdodan) >var=var(nhietdodan) >cv=sd(nhietdodan)*100/mean(nhietdodan) >c(sd,var,cv) >mean(nhietdodan)-sd(nhietdodan) >mean(nhietdodan)+sd(nhietdodan) Ta bảng đặc trưng thống kê biến: "nhietdodan" ST T Số thống kê Giá trị Số trung bình (Mean) Trung vị (Median) Độ lệch chuẩn (Sd) Phương sai (Var) Độ biến thiên (Cv) Khoảng tin cậy 95% 160 160 14.433757 208.333333 9.021098 145.5662; 174.4338 Mô tả đặc trưng thống kê biến: "aplucdan" > summary(aplucdan) >sd=sd(aplucdan) >var=var(aplucdan) >cv=sd(aplucdan)*100/mean(aplucdan) >c(sd,var,cv) >mean(aplucdan)-sd(aplucdan) >mean(aplucdan)+sd(aplucdan) Ta bảng đặc trưng thống kê biến: : "aplucdan" ST T Số thống kê Giá trị Số trung bình (Mean) Trung vị (Median) Độ lệch chuẩn (Sd) Phương sai (Var) Độ biến thiên (Cv) Khoảng tin cậy 95% 2,5 2,5 0.7216878 0.5208333 28.8675135 1.778312; 3.221688 Mô tả đặc trưng thống kê biến: " thoigiandan" >summary(thoigiandan) >sd=sd(thoigiandan) >var=var(thoigiandan) >cv=sd(thoigiandan)*100/mean(thoigiandan) >c(sd,var,cv) >mean(thoigiandan)-sd(thoigiandan) >mean(thoigiandan)+sd(thoigiandan) Ta bảng đặc trưng thống kê biến: : " thoigiandan " ST T Số thống kê Giá trị Số trung bình (Mean) Trung vị (Median) Độ lệch chuẩn (Sd) Phương sai (Var) Độ biến thiên (Cv) Khoảng tin cậy 95% 16 16 4.330127 18.750000 27.063294 11.66987; 20.33013 Mô tả đặc trưng thống kê biến: " dobenxelienketdan" >summary(dobenxelienketdan) >sd=sd(dobenxelienketdan) >var=var(dobenxelienketdan) >cv=sd(dobenxelienketdan)*100/mean(dobenxelienketdan) >c(sd,var,cv) >mean(dobenxelienketdan)-sd(dobenxelienketdan) >mean(dobenxelienketdan)+sd(dobenxelienketdan) Ta bảng đặc trưng thống kê biến: : " dobenxelienketdan " ST T Số thống kê Giá trị Số trung bình (Mean) Trung vị (Median) Độ lệch chuẩn (Sd) Phương sai (Var) Độ biến thiên (Cv) Khoảng tin cậy 95% 81,36 81,00 7.968061 63.490000 9.793586 73.39194; 89.32806 Mô tả đặc trưng thống kê biến: " dogian" >summary(dogian) >sd=sd(dogian) >var=var(dogian) >cv=sd(dogian)*100/mean(dogian) >c(sd,var,cv) >mean(dogian)-sd(dogian) >mean(dogian)+sd(dogian) Ta bảng đặc trưng thống kê biến: : " dogian " ST Số thống kê Giá trị T Số trung bình (Mean) Trung vị (Median) Độ lệch chuẩn (Sd) Phương sai (Var) Độ biến thiên (Cv) Khoảng tin cậy 95% 12,36 12,00 4.283009 18.344167 34.652177 8.076991; 16.64301 Yêu cầu 2: Vẽ biểu đồ phân bố biến - Ta sử dụng lệnh hist R để vẽ biểu đồ phân bố biến cố, để dễ dàng so sánh biến cách khách quan nhất, ta vẽ đồ thị cho tất biến cách chia ảnh thành cửa sổ với dòng cột: >op=par(mfrow=c(2,3)) - biểu đồ phân phối biến : "nhietdodan" >hist(nhietdodan, prob=T) >lines(density(nhietdodan)) - biểu đồ phân phối biến : "aplucdan" >hist(aplucdan, prob=T) >lines(density(aplucdan)) - biểu đồ phân phối biến : " thoigiandan " >hist(thoigiandan, prob=T) >lines(density(thoigiandan)) - biểu đồ phân phối biến : " dobenxelienketdan " >hist(dobenxelienketdan, prob=T) >lines(density(dobenxelienketdan)) - biểu đồ phân phối biến : " dogian " >hist(dogian, prob=T) >lines(density(dogian)) 150 160 170 180 0.02 0.00 2.0 2.5 3.0 3.5 Histogram of dogian 0.08 0.06 0.02 0.04 Density 0.03 0.02 0.01 0.00 0.00 10 15 dogian 10 12 14 16 18 20 22 thoigiandan 0.10 Histogram of dobenxelienketdan 0.04 aplucdan dobenxelienketdan 0.04 Density 1.5 nhietdodan 65 70 75 80 85 90 95 0.06 0.08 0.8 0.0 0.00 0.2 0.4 Density 0.6 0.03 0.02 0.01 Density 140 Density Histogram of thoigiandan 0.10 Histogram of aplucdan 0.04 Histogram of nhietdodan 20 Yêu cầu 3: Vẽ biểu đồ tương quan đôi biến liên tục - Vẽ biểu đồ tương quan đôi biến: "nhietdodan, "aplucdan", "thoigiandan" , "dobenxelienketdan","dogian" >dat=cbind(nhietdodan,aplucdan,thoigiandan,dobenxelienketdan,dogia n) >pairs.panels(dat) Ta biểu đồ: Nhận xét : Đường màu đỏ thể xu hướng biến đổi cặp biến Đường dóng ngang thể hệ số tương quan hai biến Đồ thị biểu thị mối liên quan cặp biến đôi Yêu cầu 4: Xác định mơ hình tuyến tính đa biến tối ưu thể ảnh hưởng các thông số dán dựng tới độ bền mối liên kết dán Mơ hình tuyến tính đa biến tối ưu thể ảnh hưởng thông số dán dựng tới độ bền mối liên kết dán mơ hình tham số giải thích nhiều liệu Ta có phương trình hồi quy: > m1=lm(dobenxelienketdan~nhietdodan) >m2=lm(dobenxelienketdan~aplucdan) >m3=lm(dobenxelienketdan ~thoigiandan) >m4=lm(dobenxelienketdan~ dogian) >m5=lm(dobenxelienketdan~nhietdodan+aplucdan+thoigiandan+dogi an) >summary(m1) p-value>2.2 e-16 , mơ hình khơng có ý nghĩa thống kê dobenxelienketdan=100,56-0.12x nhietdodan nhiệt độ tăng lên độ C thị độ bền liên kết dán giảm 0,12 N >summary(m2) p-value>2.2 e-16 , mơ hình khơng có ý nghĩa thống kê phương trình: dobenxelienketdan=108,16-10.72xaplucdan áp lực dán tăng lên 1N độ bền liên kết giảm 10,72N >summary(m3) p-value>2.2 e-16 , mô hình khơng có ý nghĩa thống kê phương trình: dobenxelienketdan=109,94-1,78xthoigiandan thời gian dán tăng 1s độ bền liên kết dán giảm 1,78 N >summary(m4) p-valuesummary(m5) Mơ hình nhietdodan R2 Sai số chuẩn aplucdan thoigianda dogian n 7.945 1.948 1.948 0.049600 4.428000 N/A 0.6088 1.100000 0.04725 0.9427 0.9427 0.9944 0.9985 Phương trình hồi quy tuyến tính đa biến: Dobenxelienketdan=86.770000-0.049600xnhietdodan4.428000xaplucdan+ 1.100000xdogian p-valueinstall.packages("BMA") >install.packages("ggplot2") > library(BMA) > library(ggplot2) > newdata =data.frame(nhietdodan,aplucdan,thoigiandan,dogian, dobenxelienketdan) > newdata=na.omit(newdata) > yvar = newdata[,1] > xvars = newdata[,-1] > reg = lm(dobenxelienketdan ~ , data=newdata) > step(reg, direction="both") > bma= bicreg(xvars, yvar, strict=FALSE, OR=20) > summary(bma) BMA trình bày kết mơ hình đánh giá tối ưu Cột thứ liệt kê biến số độc lập Cột thứ hai trình bày xác suất biến độc lập có ảnh hưởng đến độ bền liên kết mối dán Ta có "dogian" “dobenxelienketdan “có xác xuất ảnh hưởng đến độ bền mối liên kết dán 100% cịn hai biến "aplucdan" "thoigiandan" 50%.Như mơ hình với biến "dogian", “dobenxelienketdan “ mơ hình tối ưu Cột EV cột SD trình bày trị số trung bình độ lệch chuẩn hệ số cho biến độc lập Cột ước tính hệ số ảnh hưởng mơ hình Mơ hình có xác suất xuất 50% Cột ước tính hệ số ảnh hưởng mơ hình Mơ hình có xác suất xuất 50% Thể kết dạng biểu đồ > imageplot.bma(bma) Models selected by BMA aplucdan thoigiandan dogian dobenxelienketdan Model # Như mơ hình biến "dogian" “dobenxelienketdan “ mơ hình tối ưu thể ảnh hưởng các thông số dán dựng tới độ bền mối liên kết dán ... độ dao động độ bền mối liên kết dán phương trình: dobenxelienketdan=86,77-0,05xnhietdogian4,42aplucdan+1,1dogian Mơ hình giải thích 99,85 dao động độ bền liên kết mối dán Mỗi 1N tăng độ bền liên. .. thể hệ số tương quan hai biến Đồ thị biểu thị mối liên quan cặp biến đôi u cầu 4: Xác định mơ hình tuyến tính đa biến tối ưu thể ảnh hưởng các thông số dán dựng tới độ bền mối liên kết dán Mơ... 99,44 % độ dao động độ bền mối liên kết dán phương trình: dobenxelienketdan=58,43+1,85dogian độ giãn tăng lên 1cm độ bền liên kết dán tăng 1,85N >summary(m5) Mơ hình nhietdodan R2 Sai số chuẩn