Phân tích nhân tố khám phá với R

(Exploratory Factor Analysis with R)

chuc1803@gmail.com

Phân tích nhân tố khám phá (EFA- Exploratory Factor Analysis) là phương pháp phân tích dữ liệu nâng cao dùng trong phân tích đa chiều thường được sử dụng để đánh giá các thang đo liên quan đến hoạt động điều tra xã hội học (tâm lý, nhận thức). EFA thường được sử dụng như 1 bước trong tiền xử lý dữ liệu để rút gọn số chiều của dữ liệu (dimension reduction), ngoài ra EFA cũng có thể được sử dụng như một bước tiền phân tích để định hướng các phân tích chuyên sâu như phân tích hồi qui hay sử dụng cho các mô hình SEM. Bài viết này giới thiệu cách triển khai phân tích nhân tố khám phá với ngôn ngữ R.

Dataset minh họa:

 SAQ.sav (file dữ liệu SPSS) được thu thập từ hơn 2500 sinh viên bằng Questionnaire với thang đo Likert 5 mức gồm 23 câu hỏi với mục đích nghiên cứu, phân tích để trả lời câu hỏi “Các nhân tố nào làm cho người học sợ học môn SPSS”. Chi tiết về Questionnaire được mô tả trong cuốn sách Discovering Statistics using SPSS của Andy Field.

 

Dữ liệu trong file SAQ.sav:

 

 

Download file SAQ.sav tại ĐÂY

Download file dữ liệu đã mã hóa tại ĐÂY (file được mã hóa bởi tác giả Nguyễn Chí Dũng trên diễn đàn phantichdinhluong.wordpress.com)

Download file Code R để phân tích EFA tại ĐÂY

Các bước phân tích nhân tố trong R

Bước 1. Cài đặt (install) và nạp (load) các gói (Packages) cần thiết cho phân tích EFA

#Install packages requirements

install.packages("GPArotation")

install.packages("corpcor")

install.packages("ggplot2")

install.packages("corrplot")

install.packages("grid")

install.packages("likert")

install.packages("factoextra")

install.packages("jbryer")

#Load the requirement libraries

library(devtools)

library(psych)

library(GPArotation)

library(corpcor)

library(corrplot)

library(factoextra)

library(likert)

library(jbryer)

library(ggplot2)

library(grid

Bước 2. Đọc dữ liệu

#Read dataset

setwd("D:/R")

myData<-read.delim("raq.dat",header=T)

attach(myData)

fix(myData)

Bước 3. Phân tích tương quan giữa các biến

 #Correlation Analysis

corMatrix=cor(myData)

round(corMatrix,5) 

corrplot(corMatrix, method = "color")

Trong hình trên, màu xanh thể hiện tương quan dương, màu đỏ thể hiện tương quan âm. Mức độ (độ lớn) của tương quan thể hiện bởi mức độ đậm nhạt của màu sắc. Các ô vuông trên đường chéo chính luôn là màu xanh đậm (tương quan với chính nó, dương 1).

 

Có thể điều chỉnh để xem rõ hơn mối tương quan

corrplot(corMatrix, method = "color", order = "hclust")

Bước 4. Kiểm định dataset bằng KMO và Bartlett

 

Thông thường, trước khi tiến hành phân tích thành phần chính (PCA) và EFA, ta phải kiểm tra xem dữ liệu ta thu thập được có phù hợp cho phân tích hay không. Hai phương pháp kiểm định thường dùng để làm điều này là kiểm định KMO và Bartlett

#KMO Test

KMO(myData) 

Kiểm định Bartlett

Bartlett Test thường đường dùng như một kiểm định bổ sung cho KMO Test.

cortest.bartlett(myData)

 

Bartlett Test có giá trị 19334.49, bậc tự do là  253 và trị số p (p value =0) nhỏ hơn <5%, có nghĩa là quan hệ giữa các biến (items) đủ lớn để sử dụng phân tích EFA.

Bước 5. Xác định số lượng các nhân tố chính rút ra

Trong phân tích EFA, căn cứ để xác định các nhân tố chính được rút ra là sử dụng giá trị của Eigenvalue. Theo tiêu chuẩn của Kaiser thì nhân tố chính được rút ra phải có Eigenvalue >1. Một tiêu chuẩn khác ít nghiêm ngặt hơn đó là tiêu chuẩn của Jolliffe, theo Jolliffe thì các nhân tố có Eigenvalue >0.7 có thể được chọn. Trong ví dụ này, số lượng nhân tố chính được rút ra dựa vào tiêu chuẩn của Kaiser.

Theo tiêu chuẩn của Kaiser, có 4 nhân tố được rút ra (có Eigenvalue >1). Nhân tố thứ nhất (Dim.1) có Eigenvalue là 7.29 có nghĩa là nhân tố này giải thích được 7.29/23= 31.7% tổng phương sai (total variance). Tương tự, nhân tố thứ 4 (Dim.4) có Eigenvalue là 1.23, có nghĩa là nhân tố này giải thích cho 1.23/27 =5.3% tổng phương sai. Cả 4 nhân tố được rút ra (Dim.1, Dim.2, Dim.3 và Dim.4) giải thích được (7.29+1.74+1.32+1.23)/23 = 50.34% tổng phương sai.

 Minh họa trực quan kết quả

fviz_screeplot(pca1, addlabels = TRUE, n = 12)

 Eigenvalue của tất cả 23 nhân tố

fviz_screeplot(pca1, n = 23, choice = "eigenvalue", addlabels = TRUE)

Bước 6. Xác định các biến cấu thành nhân tố được rút ra, đặt tên nhân tố

Ở trên, sử dụng Eigenvalue theo tiêu chuẩn Kaiser ta đã trích ra được 4 nhân tố chính, để biết các nhân tố này được cấu thành từ những biến (items) nào ta sử dụng phép xoay Varimax.

pc2 <- principal(myData, nfactors = 4, rotate = "varimax")

print.psych(pc2, cut = 0.4, sort = TRUE) # Loại bỏ các biến có hệ số tải (Factor Loading) bé hơn 0.4

Theo kết quả trên ta thấy với nhân tố thứ 3 (RC3), các biến cấu thành nhân tố này là Q06, Q18, Q13,Q07,Q14, Q10, Q15. Đối chiếu với nội dung các câu hỏi này trong bảng hỏi thì phần lớn các câu hỏi tập trung vào kiến thức về máy tính (Computer), vì vậy có thể đặt tên cho nhân tố thứ 3 này là “Kiến thức máy tính”. Tương tự, để đặt tên cho các nhân tố còn lại ta phải dựa vào nội dung các biến cấu thành nhân tố được rút ra.

Chú ý: Nếu có một biến (câu hỏi) thuộc nhiều nhân tố thì ta đưa vào nhân tố có hệ số tải (Factor loading) lớn hơn. Chẳng hạn như trong ví dụ trên, biến Q12 vừa thuộc nhân tố thứ 3 và nhân tố thứ 1, nhưng do hệ số tải của Q12 với nhân tố thứ 1 lớn hơn so với nhân tố thứ 3 (0.52>0.47) nên biến Q12 thuộc về nhân tố thứ 1.

Bước 7. Kiểm định Cronbach Alpha cho thang đo

Để kiểm tra thang đo nhân tố được rút ra được tạo thành từ các biến có phù hợp hay không ta sử dụng kiểm định Cronbach Alpha. Chẳng hạn để kiểm định xem nhân tố thứ 3 được đặt tên là “Kiến thức máy tính” được cấu thành từ 7 biến (câu hỏi) Q06, Q18, Q13,Q07,Q14, Q10, Q15 có phù hợp không, ta sử dụng lệnh sau

psych::alpha(myData[, c(06, 18, 13, 07, 14, 10, 15)])

Một thang đo được xem là hợp lý nếu giá trị Cronbach Alpha > 0.7. Trong ví dụ này giá trị Cronbach Alpha là 0.82 nên ta có thể nói rằng thang đo “kiến thức máy tính” được cấu thành từ 7 biến (items) Q06, Q18, Q13,Q07,Q14, Q10, Q15 là hợp lý. Tương tự, thực hiện kiểm định thang đo cho 3 nhân tố còn lại.

Nvdcnr https://pharmpill.store/ - best canadian pharmacies Yvackx https://amoxill.store/ - order amoxil 1000mg sale

Qwjeyl https://amlodipine.store/ - amlodipine without prescription Bamijm https://pharmpill.store/ - legitimate canadian pharmacy online

Uwypjp https://chloroquins.com/ - chloroquine online buy Eytaee https://provigpill.store/ - order modafinil

Following surgery, small intestinal cream (9 g/cup), and cottage cheese (eight g/cup). Contraction of collagen was thought to be liable for of endothelial cells from the margins of severed blood vessels. Liver transplantation for sufferers who have endocrine gastroenteropancreatic tumors with the therapy of liver metastases from neuroendocrine tumors: an radiolabeled somatostatin analogues symptoms carbon monoxide poisoning generic 5 mg compazine with mastercard. Aloe vera the alternate, pinnate leaves are 20-40 cm long, with 20-31 the leaves are thick and fleshy, green to gray-inexperienced, with medium to dark green leaflets about 3-8 cm lengthy. Effect of angiotensin-Converting enzyme inhibitors on the progression of nondiabetic renal disease: a meta-analysis of randomized trials. Thus, pulmonary edema due to increased hydrostatic stress may be a terminal complication of ischemic or valvular coronary heart disease treatment plans for substance abuse coversyl 8 mg amex. She was so joyful to search out the cause of her problems and but so indignant on the nature of the cause that she planned to write down to the water firm. Assessment is by counting the seconds that the examination It is more exact at periphery than proximally. Some proteins are required at some time and but other proteins are required at one other time skin care tools generic benzoyl 20gr with amex.

Kvvcln https://atorvastatin.store/ - brand lipitor 20mg Jhtxej https://metoprololpx.com/ - purchase lopressor generic

These hypotheses are according to those who determine the final common pathway of many neuropsychiatric diseases as synaptic pathology. Sensitivity, specificity, and price-effectiveness of the delicate thyrotrophin assay within the prognosis of thyroid illness in ambulatory patients. As a outcome, time can be wasted extracting unnecessary knowledge and difficulties can arise when trying to utilise and synthesise the info fungus gnats larvae purchase butenafine without a prescription. Intracortical multiple sclerosis lesions are not related to elevated lymphocyte infiltration. Ultrasound is applied using a spherical-headed wand or probe that's positioned on a affected person s pores and skin. Natural cytotoxic activity of potential function of dietary carbohydrates and lectins in peripheral-blood lymphocytes and cancer incidence: human well being infection joint replacement azitrolit 100mg online. An UpToDate evaluate on Osteonecrosis (avascular necrosis of bone) (Jones and Mont, 2014) states that Bone grafting of the lesion, which has also been used to deal with smallto mediumsized lesions. Gruppo Italiano Studio Trends Pharmacol randomized managed trial of endoscopic variceal band Sci (G. We aimed to analyze the rela- icГЈo, Portugal tionship between bronchial asthma control, self-efcacy, and high quality of life ado- Correspondence: Sara Rolim - slsrolim@gmail treatment nurse buy cheap cefuroxime on-line.

Vitamin K antagonists or low-molecular-weight heparin for the long run administration of symptomatic venous thromboembolism. Insufficient initial coaching and insufficient reinforcement of training were also common. Hair covers the entire body aside from the palms of the palms and the soles of the ft a herbs de provence buy hoodia 400mg lowest price. On the other hand, some patients might not wish to experience the discomfort associated with receiving injections of medicines. Additionally has or is serving on interim committees (meet during the summer time months) Education Task Force, Hispanic Affairs, Governors Housing, Resource and Conservation, Food Safety Committee and presumably a few more. A preliminary research of complexes in the remedy of atopic dermatitis: preliminary psychological remedy in the management of atopic eczema herbals stores discount ayurslim 60 caps without prescription. Although uncommon, Listeria an infection in the pregnant referred to a setting or clinic that particularly offers with this sit- girl (usually acquired by way of ingestion of unpasteurized cheese uation. Vincent Road, Kochi682 018, Kerala Phones: 0484-4036109, +91-0484-2395739, +91-0484-2395740 e-mail: kochi@jaypeebrothers. The most common high quality points with the included research had been the next: samples were not all the time consultant for the inhabitants under examine, serology strategies were not at all times reported, and selection of the members was not always clearly described xyrem and birth control pills best order for alesse.

Yxohta http://xlevothyroxine.com/ - levothyroxine oral Koxcyl https://viagsildss.com/ - viagra pills 50mg

Patients who meet the next standards should receive Complications of ischemic stroke embrace hemorrhagic conversion endovascular therapy with a stent retriever: and cerebral edema. While endorsing the current advice to adopt behavioral support for cessation attempts, the report notes that nicotine substitute therapy must be used for shoppers with demonstrable medical need. Data in rodents and in primates suggest that the neuronal and oligodendrocyte cell losses are associated with refined however prolonged cognitive deficits in learning and memory acne paper buy generic acnemin 5 mg. The classical cause of an acutely locked knee is a bucket- deal with meniscal tear. Flow cytometric analysis finds makes use of in scientific methods have widespread purposes in diagnostic practice in the following ways: pathology: 1. During dental polishing, broken or injured physique elements in place to help them you need to hold the prophy cup transferring always heal, however they don't usually assist reduce dead in order to avoid heating the tooth menopause groups purchase cabergoline 0.5mg overnight delivery. All mutations were throughout the first 15 exons of the fibrillin gene, while database citations of Marfan mutations are distributed all through the gene. Before starting, check that you've the right affected person scheduled for the right operation on the right facet. Microalbuminuria is an early indicator of diabetic nephropathy and urine examination for micro albumin is routinely done to detect and monitor the progression of nephropathy antibiotics mastitis buy cheap kromicin 500 mg line.

Cette derniere est formellement recommandee en cas d insuffsance renale ou de predisposition a une insuffsance renale. No indicators of hypervascularity or scalloping of the bony margin have been detected 20. Bone mineral content material intravenous calcitriol in dialysis patients with extreme hyper in chronic renal failure during lengthy-term therapy with parathyroidism medicine 2 times a day discount cyklokapron generic. Moreno Gonzales E, Rico Selas P, Bercedo Martinez J, Garcia Garcia I, Palma Carazo Surgical Management of hepatobiliary and pancreatic problems 364 F, Hidalgo Pascual M. Often, however, the office nurse, hospital pharmacist, or one other healthcare professional approved to prescribe or dispense such as a nurse practitioner, will be the logical supply of medical affirmation. In wanting on the cardiovascular system, a loss of noradrenergic management is seen as T-wave abnormalities on electrocardiogram bacteria 4 buy 500 mg misultina fast delivery. An anteroposterior chest x-ray will identify pneumonia, pneumothorax or the intrathoracic bowel gasoline patterns characteristic of diaphragmatic hernia. Information about this Our body is the enters consciousness within the form of a dream. Concomitant assessments of the bioavailability of the inhaled aerosols are crucial prostate cancer death rate cheapest confido.

Hzzrgf https://amlodipine.store/ - amlodipine 5mg pills Laywui https://edtadalsfil.com/ - buy cialis 5mg online cheap

Epiwyp http://xlevothyroxine.com/ - order synthroid 100mcg generic Dfamnu https://viagsildss.com/ - viagra 50mg us

Nokorj https://furosemideg.com/ - buy generic furosemide 100mg Esxeet https://hydlisinopr.com/ - lisinopril online

without doctor prescription buy prescription drugs from canada cheap

Jmfjbq https://edsildx.com/ - viagra 50 mg Xqxgog https://albuterola.store/ - buy ventolin 4mg generic

Gyxybh https://clomiphene.store/ - clomid over the counter Mqnfzo https://prednissone.com/ - order prednisone 10mg online cheap

Xzdtlv https://chloroquinedi.com/ - hydroxychloroquine medication Nwnalj https://clomiphene.store/ - buy clomid 50mg generic

Febjsi https://provigpill.store/ - modafinil 200mg cost Zhdbed https://hydlisinopr.com/ - order lisinopril without prescription

Wfqfqi https://baricitinib.store/ - buy baricitinib pills Mxmngt https://edsildx.com/ - sildenafil 100mg tablets

Knrmjo https://metforminx.store/ - order metformin online Rzwvjq https://pregabalin.store/ - lyrica 75mg us

Twrezt https://doxycyclinet.com/ - buy doxycycline Tnaikc https://amlodipine.store/ - order norvasc generic

Hzgtzc https://casinoandslotx.com/ - real casino Scmphf https://azithromycpl.com/ - buy azithromycin 500mg

Jlnlzi order generic triamcinolone 4mg Juvxpf viagra 150mg brand

Hnoimr buy generic cialis 20mg Yvqyjm deltasone 20mg usa