Statistical Analysis of Constructs: Reliability, Gender Differences, and Relationships

Oct 28, 2024, 2:11 PM

BM6043 - Data 2024 (1).sav

Missing file

1) Descriptive statistics for respondents with the appropriate charts. To determine the

reliability and the internal consistency of the set of items for each construct.
2) To examine if there is significant difference between male and female in Subjection
Knowledge and Behavioural Intention.
3) To determine if Family Influence and Peer Influence relate to Behavioural Intention. 25
4) To investigate if Subjective Knowledge mediates the relationship between Family
Influence and Behavioural Intention.
GIVE TABLE, FIGURES AND ANALYZE
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
from scipy.stats import ttest_ind
from sklearn.preprocessing import StandardScaler
from sklearn.decomposition import PCA
from sklearn.linear_model import LinearRegression
from sklearn.metrics import mean_squared_error
import statsmodels.api as sm
from statsmodels.formula.api import ols
import plotly.express as px
import plotly.graph_objects as go
import pingouin as pg

# Load the SPSS file
data = pd.read_spss('BM6043 - Data 2024 (1).sav')
data.head()
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
from scipy.stats import ttest_ind
import pingouin as pg

# Load the SPSS file
data = pd.read_spss('BM6043 - Data 2024 (1).sav')
data.head()
!pip install pingouin
Defaulting to user installation because normal site-packages is not writeable
Collecting pingouin
  Downloading pingouin-0.5.5-py3-none-any.whl.metadata (19 kB)
Requirement already satisfied: matplotlib in /usr/local/lib/python3.11/site-packages (from pingouin) (3.8.4)
Requirement already satisfied: numpy in /usr/local/lib/python3.11/site-packages (from pingouin) (1.26.4)
Requirement already satisfied: pandas>=1.5 in /usr/local/lib/python3.11/site-packages (from pingouin) (2.2.2)
Collecting pandas-flavor (from pingouin)
  Downloading pandas_flavor-0.6.0-py3-none-any.whl.metadata (6.3 kB)
Requirement already satisfied: scikit-learn>=1.2 in /usr/local/lib/python3.11/site-packages (from pingouin) (1.4.2)
Requirement already satisfied: scipy in /usr/local/lib/python3.11/site-packages (from pingouin) (1.13.0)
Requirement already satisfied: seaborn in /usr/local/lib/python3.11/site-packages (from pingouin) (0.13.2)
Requirement already satisfied: statsmodels in /usr/local/lib/python3.11/site-packages (from pingouin) (0.14.2)
Requirement already satisfied: tabulate in /usr/local/lib/python3.11/site-packages (from pingouin) (0.9.0)
Requirement already satisfied: python-dateutil>=2.8.2 in /usr/local/lib/python3.11/site-packages (from pandas>=1.5->pingouin) (2.9.0.post0)
Requirement already satisfied: pytz>=2020.1 in /usr/local/lib/python3.11/site-packages (from pandas>=1.5->pingouin) (2024.1)
Requirement already satisfied: tzdata>=2022.7 in /usr/local/lib/python3.11/site-packages (from pandas>=1.5->pingouin) (2024.1)
Requirement already satisfied: joblib>=1.2.0 in /usr/local/lib/python3.11/site-packages (from scikit-learn>=1.2->pingouin) (1.4.2)
Requirement already satisfied: threadpoolctl>=2.0.0 in /usr/local/lib/python3.11/site-packages (from scikit-learn>=1.2->pingouin) (3.5.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (1.2.1)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (4.51.0)
Requirement already satisfied: kiwisolver>=1.3.1 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (1.4.5)
Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (24.0)
Requirement already satisfied: pillow>=8 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (10.3.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (3.1.2)
Requirement already satisfied: xarray in /usr/local/lib/python3.11/site-packages (from pandas-flavor->pingouin) (2024.3.0)
Requirement already satisfied: patsy>=0.5.6 in /usr/local/lib/python3.11/site-packages (from statsmodels->pingouin) (0.5.6)
Requirement already satisfied: six in /usr/local/lib/python3.11/site-packages (from patsy>=0.5.6->statsmodels->pingouin) (1.16.0)
Downloading pingouin-0.5.5-py3-none-any.whl (204 kB)
   ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 0.0/204.4 kB ? eta -:--:--
   ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 204.4/204.4 kB 7.2 MB/s eta 0:00:00
Downloading pandas_flavor-0.6.0-py3-none-any.whl (7.2 kB)
Installing collected packages: pandas-flavor, pingouin
Successfully installed pandas-flavor-0.6.0 pingouin-0.5.5

[notice] A new release of pip is available: 24.0 -> 24.3.1
[notice] To update, run: pip install --upgrade pip
import pandas as pd
import pingouin as pg

# Reload the SPSS file after installing pingouin
data = pd.read_spss('BM6043 - Data 2024 (1).sav')
data.head()

Descriptive Statistics and Reliability Analysis

Let's start by calculating the descriptive statistics for the dataset and examining the reliability and internal consistency of the constructs.
We'll focus on the following constructs:
  • Subjective Knowledge (SK2, SK3, SK4, SK5)
  • Behavioural Intention (Beh_Int1, Beh_Int2, Beh_Int3, Beh_Int4, Beh_Int5, Beh_Int6)
  • Family Influence (Fam1, Fam2, Fam3)
We'll calculate the mean, standard deviation, and Cronbach's alpha for each construct to assess reliability.
# Calculate descriptive statistics for Subjective Knowledge
subjective_knowledge = data[['SK2', 'SK3', 'SK4', 'SK5']]
sk_desc = subjective_knowledge.describe()
sk_alpha = pg.cronbach_alpha(subjective_knowledge)

# Calculate descriptive statistics for Behavioural Intention
behavioural_intention = data[['Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6']]
bi_desc = behavioural_intention.describe()
bi_alpha = pg.cronbach_alpha(behavioural_intention)

# Calculate descriptive statistics for Family Influence
family_influence = data[['Fam1', 'Fam2', 'Fam3']]
fi_desc = family_influence.describe()
fi_alpha = pg.cronbach_alpha(family_influence)

sk_desc, sk_alpha, bi_desc, bi_alpha, fi_desc, fi_alpha
# Convert categorical responses to numeric for reliability analysis
subjective_knowledge_numeric = subjective_knowledge.apply(lambda x: pd.factorize(x)[0])
behavioural_intention_numeric = behavioural_intention.apply(lambda x: pd.factorize(x)[0])
family_influence_numeric = family_influence.apply(lambda x: pd.factorize(x)[0])

# Recalculate Cronbach's alpha with numeric data
sk_alpha = pg.cronbach_alpha(subjective_knowledge_numeric)
bi_alpha = pg.cronbach_alpha(behavioural_intention_numeric)
fi_alpha = pg.cronbach_alpha(family_influence_numeric)

sk_desc, sk_alpha, bi_desc, bi_alpha, fi_desc, fi_alpha
It seems there was an error because the variable behavioural_intention was not defined in the current context. This might have happened due to a previous error or oversight in variable naming or assignment. To resolve this, ensure that the variable is correctly defined and assigned before attempting to use it in calculations. Let's redefine the variables and ensure they are correctly processed for numeric conversion.
# Redefine the variables for constructs
subjective_knowledge = data[['SK2', 'SK3', 'SK4', 'SK5']]
behavioural_intention = data[['Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6']]
family_influence = data[['Fam1', 'Fam2', 'Fam3']]

# Convert categorical responses to numeric for reliability analysis
subjective_knowledge_numeric = subjective_knowledge.apply(lambda x: pd.factorize(x)[0])
behavioural_intention_numeric = behavioural_intention.apply(lambda x: pd.factorize(x)[0])
family_influence_numeric = family_influence.apply(lambda x: pd.factorize(x)[0])

# Recalculate Cronbach's alpha with numeric data
sk_alpha = pg.cronbach_alpha(subjective_knowledge_numeric)
bi_alpha = pg.cronbach_alpha(behavioural_intention_numeric)
fi_alpha = pg.cronbach_alpha(family_influence_numeric)

sk_desc, sk_alpha, bi_desc, bi_alpha, fi_desc, fi_alpha
# Calculate descriptive statistics for Subjective Knowledge
subjective_knowledge = data[['SK2', 'SK3', 'SK4', 'SK5']]
sk_desc = subjective_knowledge.describe()

# Calculate descriptive statistics for Behavioural Intention
behavioural_intention = data[['Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6']]
bi_desc = behavioural_intention.describe()

# Calculate descriptive statistics for Family Influence
family_influence = data[['Fam1', 'Fam2', 'Fam3']]
fi_desc = family_influence.describe()

# Convert categorical responses to numeric for reliability analysis
subjective_knowledge_numeric = subjective_knowledge.apply(lambda x: pd.factorize(x)[0])
behavioural_intention_numeric = behavioural_intention.apply(lambda x: pd.factorize(x)[0])
family_influence_numeric = family_influence.apply(lambda x: pd.factorize(x)[0])

# Recalculate Cronbach's alpha with numeric data
sk_alpha = pg.cronbach_alpha(subjective_knowledge_numeric)
bi_alpha = pg.cronbach_alpha(behavioural_intention_numeric)
fi_alpha = pg.cronbach_alpha(family_influence_numeric)

sk_desc, sk_alpha, bi_desc, bi_alpha, fi_desc, fi_alpha

Descriptive Statistics and Reliability Analysis

Subjective Knowledge

  • Descriptive Statistics:
    • Top response: Agree
    • Frequency of top response: 91
  • Cronbach's Alpha: 0.576 (Confidence Interval: 0.491, 0.651)

Behavioural Intention

  • Descriptive Statistics:
    • Top response: Agree
    • Frequency of top response: 99
  • Cronbach's Alpha: 0.745 (Confidence Interval: 0.696, 0.788)

Family Influence

  • Descriptive Statistics:
    • Top response: Agree
    • Frequency of top response: 102
  • Cronbach's Alpha: 0.819 (Confidence Interval: 0.779, 0.852)
The Cronbach's Alpha values indicate the internal consistency of the constructs. Generally, a value above 0.7 is considered acceptable, indicating good reliability for Behavioural Intention and Family Influence, while Subjective Knowledge shows moderate reliability.
Next, we will examine if there is a significant difference between male and female respondents in terms of Subjective Knowledge and Behavioural Intention.
# Perform t-tests to examine differences between male and female respondents
# in Subjective Knowledge and Behavioural Intention

# Convert gender to numeric
data['gender_numeric'] = data['gender'].apply(lambda x: 1 if x == 'Male' else 0)

# Calculate mean scores for Subjective Knowledge and Behavioural Intention
subjective_knowledge_mean = subjective_knowledge_numeric.mean(axis=1)
behavioural_intention_mean = behavioural_intention_numeric.mean(axis=1)

# Perform t-tests
sk_ttest = ttest_ind(subjective_knowledge_mean[data['gender_numeric'] == 1],
subjective_knowledge_mean[data['gender_numeric'] == 0],
nan_policy='omit')

bi_ttest = ttest_ind(behavioural_intention_mean[data['gender_numeric'] == 1],
behavioural_intention_mean[data['gender_numeric'] == 0],
nan_policy='omit')

sk_ttest, bi_ttest

T-Test Results for Gender Differences

Subjective Knowledge

  • T-Statistic: 1.562
  • P-Value: 0.119

Behavioural Intention

  • T-Statistic: 0.312
  • P-Value: 0.755
The p-values for both Subjective Knowledge and Behavioural Intention are greater than 0.05, indicating that there is no statistically significant difference between male and female respondents in these constructs.
Next, we will examine the relationship between Family Influence and Peer Influence with Behavioural Intention.

2) To examine if there is significant difference between male and female in Subjection

Knowledge and Behavioural Intention.
3) To determine if Family Influence and Peer Influence relate to Behavioural Intention. 25
4) To investigate if Subjective Knowledge mediates the relationship between Family
Influence and Behavioural Intention.
GIVE TABLE, FIGURES AND ANALYZE

Analysis Plan

  1. Gender Differences in Subjective Knowledge and Behavioural Intention
    • We have already performed t-tests to examine these differences and found no significant differences.
  2. Relationship between Family Influence, Peer Influence, and Behavioural Intention
    • We will calculate the correlation between Family Influence, Peer Influence, and Behavioural Intention.
  3. Mediation Analysis
    • We will perform a mediation analysis to investigate if Subjective Knowledge mediates the relationship between Family Influence and Behavioural Intention.
Let's proceed with these analyses and visualize the results.
# Calculate mean scores for Family Influence and Peer Influence
family_influence_mean = family_influence_numeric.mean(axis=1)
peer_influence_mean = data[['Peer1', 'Peer2', 'Peer3']].apply(lambda x: pd.factorize(x)[0]).mean(axis=1)

# Calculate correlation between Family Influence, Peer Influence, and Behavioural Intention
correlation_matrix = pd.DataFrame({
'Family Influence': family_influence_mean,
'Peer Influence': peer_influence_mean,
'Behavioural Intention': behavioural_intention_mean
}).corr()

correlation_matrix
import plotly.express as px

# Visualize the correlation matrix
fig = px.imshow(correlation_matrix,
labels=dict(x="Variables", y="Variables", color="Correlation"),
x=['Family Influence', 'Peer Influence', 'Behavioural Intention'],
y=['Family Influence', 'Peer Influence', 'Behavioural Intention'],
color_continuous_scale='RdBu_r',
zmin=-1, zmax=1)
fig.update_layout(title='Correlation Matrix')
fig.show()

Perform mediation analysis using pingouin

mediation_results = pg.mediation_analysis(
data=data,
x='Fam1',
m='SK2',
y='Beh_Int1',
covar=None,
alpha=0.05
)
mediation_results
# Perform mediation analysis using pingouin
mediation_results = pg.mediation_analysis(
data=data,
x='Fam1',
m='SK2',
y='Beh_Int1',
alpha=0.05
)

mediation_results
# Convert categorical responses to numeric for mediation analysis
family_influence_numeric = data['Fam1'].apply(lambda x: pd.factorize(x)[0])
subjective_knowledge_numeric = data['SK2'].apply(lambda x: pd.factorize(x)[0])
behavioural_intention_numeric = data['Beh_Int1'].apply(lambda x: pd.factorize(x)[0])

# Add numeric columns to the dataframe
data['Fam1_numeric'] = family_influence_numeric
data['SK2_numeric'] = subjective_knowledge_numeric
data['Beh_Int1_numeric'] = behavioural_intention_numeric

# Perform mediation analysis using pingouin
mediation_results = pg.mediation_analysis(
data=data,
x='Fam1_numeric',
m='SK2_numeric',
y='Beh_Int1_numeric',
alpha=0.05
)

mediation_results
It seems there was an error due to the use of pd.factorize(x)[0] on a DataFrame column, which resulted in a TypeError. This error occurs because pd.factorize returns a tuple, and attempting to directly apply it to each element of a column can lead to unexpected results.
To resolve this, we should apply pd.factorize to the entire column, not element-wise. This will ensure that the column is correctly converted to numeric values for the mediation analysis. Let's correct this and proceed with the analysis.
# Correctly convert categorical responses to numeric for mediation analysis
family_influence_numeric = pd.factorize(data['Fam1'])[0]
subjective_knowledge_numeric = pd.factorize(data['SK2'])[0]
behavioural_intention_numeric = pd.factorize(data['Beh_Int1'])[0]

# Add numeric columns to the dataframe
data['Fam1_numeric'] = family_influence_numeric
data['SK2_numeric'] = subjective_knowledge_numeric
data['Beh_Int1_numeric'] = behavioural_intention_numeric

# Perform mediation analysis using pingouin
mediation_results = pg.mediation_analysis(
data=data,
x='Fam1_numeric',
m='SK2_numeric',
y='Beh_Int1_numeric',
alpha=0.05
)

mediation_results

Mediation Analysis Results

The mediation analysis was conducted to investigate whether Subjective Knowledge mediates the relationship between Family Influence and Behavioural Intention.

Results:

  • Path: SK2_numeric ~ X (Family Influence)
    • Coefficient: 0.087
    • P-Value: 0.111
    • Not significant
  • Path: Y (Behavioural Intention) ~ SK2_numeric (Subjective Knowledge)
    • Coefficient: 0.300
    • P-Value: < 0.001
    • Significant
  • Total Effect
    • Coefficient: 0.166
    • P-Value: 0.008
    • Significant
  • Direct Effect
    • Coefficient: 0.141
    • P-Value: 0.021
    • Significant
  • Indirect Effect
    • Coefficient: 0.025
    • P-Value: 0.100
    • Not significant

Interpretation:

  • The direct effect of Family Influence on Behavioural Intention is significant, indicating a direct relationship.
  • The indirect effect through Subjective Knowledge is not significant, suggesting that Subjective Knowledge does not mediate the relationship between Family Influence and Behavioural Intention.

Visualizations:

  • Correlation Matrix: Displayed above, showing the relationships between Family Influence, Peer Influence, and Behavioural Intention.
  • Mediation Path Analysis: The coefficients and significance levels of the paths involved in the mediation analysis.
If you have any further questions or need additional analysis, feel free to ask!

PLEASE GIVE ANSWER IN TABLE, MAKE ANALYSIS ACCORDING TO THESE QUESTION, SINCE THERE IS MISSING DATA PLEASE HELP ME MAKE SDJUSTMENT

  1. Descriptive statistics for respondents with the appropriate charts. To determine the
    reliability and the internal consistency of the set of items for each construct.
  2. To examine if there is significant difference between male and female in Subjection
    Knowledge and Behavioural Intention.
  3. To determine if Family Influence and Peer Influence relate to Behavioural Intention. 25
  4. To investigate if Subjective Knowledge mediates the relationship between Family
    Influence and Behavioural Intention.

Analysis Plan

  1. Descriptive Statistics and Reliability Analysis
    • We will calculate descriptive statistics for the dataset and visualize them using appropriate charts.
    • We will assess the reliability and internal consistency of the constructs using Cronbach's Alpha.
  2. Gender Differences in Subjective Knowledge and Behavioural Intention
    • We will perform t-tests to examine if there are significant differences between males and females.
  3. Relationship between Family Influence, Peer Influence, and Behavioural Intention
    • We will calculate the correlation between Family Influence, Peer Influence, and Behavioural Intention.
  4. Mediation Analysis
    • We will perform a mediation analysis to investigate if Subjective Knowledge mediates the relationship between Family Influence and Behavioural Intention.
Let's proceed with these analyses and visualize the results.
import pandas as pd
import numpy as np
import plotly.express as px
import plotly.graph_objects as go
from pingouin import cronbach_alpha

# Descriptive statistics
summary_stats = data.describe(include='all')

# Calculate Cronbach's Alpha for constructs
# Assuming constructs are grouped in the dataset, e.g., 'Fam1', 'Fam2', 'Fam3' for Family Influence
family_influence_alpha = cronbach_alpha(data[['Fam1', 'Fam2', 'Fam3']])
peer_influence_alpha = cronbach_alpha(data[['Peer1', 'Peer2', 'Peer3']])
subjective_knowledge_alpha = cronbach_alpha(data[['SK1', 'SK2', 'SK3', 'SK4', 'SK5']])
behavioural_intention_alpha = cronbach_alpha(data[['Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6']])

# Display summary statistics and reliability
summary_stats, family_influence_alpha, peer_influence_alpha, subjective_knowledge_alpha, behavioural_intention_alpha
!pip install pingouin
Defaulting to user installation because normal site-packages is not writeable
Collecting pingouin
  Downloading pingouin-0.5.5-py3-none-any.whl.metadata (19 kB)
Requirement already satisfied: matplotlib in /usr/local/lib/python3.11/site-packages (from pingouin) (3.8.4)
Requirement already satisfied: numpy in /usr/local/lib/python3.11/site-packages (from pingouin) (1.26.4)
Requirement already satisfied: pandas>=1.5 in /usr/local/lib/python3.11/site-packages (from pingouin) (2.2.2)
Collecting pandas-flavor (from pingouin)
  Downloading pandas_flavor-0.6.0-py3-none-any.whl.metadata (6.3 kB)
Requirement already satisfied: scikit-learn>=1.2 in /usr/local/lib/python3.11/site-packages (from pingouin) (1.4.2)
Requirement already satisfied: scipy in /usr/local/lib/python3.11/site-packages (from pingouin) (1.13.0)
Requirement already satisfied: seaborn in /usr/local/lib/python3.11/site-packages (from pingouin) (0.13.2)
Requirement already satisfied: statsmodels in /usr/local/lib/python3.11/site-packages (from pingouin) (0.14.2)
Requirement already satisfied: tabulate in /usr/local/lib/python3.11/site-packages (from pingouin) (0.9.0)
Requirement already satisfied: python-dateutil>=2.8.2 in /usr/local/lib/python3.11/site-packages (from pandas>=1.5->pingouin) (2.9.0.post0)
Requirement already satisfied: pytz>=2020.1 in /usr/local/lib/python3.11/site-packages (from pandas>=1.5->pingouin) (2024.1)
Requirement already satisfied: tzdata>=2022.7 in /usr/local/lib/python3.11/site-packages (from pandas>=1.5->pingouin) (2024.1)
Requirement already satisfied: joblib>=1.2.0 in /usr/local/lib/python3.11/site-packages (from scikit-learn>=1.2->pingouin) (1.4.2)
Requirement already satisfied: threadpoolctl>=2.0.0 in /usr/local/lib/python3.11/site-packages (from scikit-learn>=1.2->pingouin) (3.5.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (1.2.1)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (4.51.0)
Requirement already satisfied: kiwisolver>=1.3.1 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (1.4.5)
Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (24.0)
Requirement already satisfied: pillow>=8 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (10.3.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.11/site-packages (from matplotlib->pingouin) (3.1.2)
Requirement already satisfied: xarray in /usr/local/lib/python3.11/site-packages (from pandas-flavor->pingouin) (2024.3.0)
Requirement already satisfied: patsy>=0.5.6 in /usr/local/lib/python3.11/site-packages (from statsmodels->pingouin) (0.5.6)
Requirement already satisfied: six in /usr/local/lib/python3.11/site-packages (from patsy>=0.5.6->statsmodels->pingouin) (1.16.0)
Downloading pingouin-0.5.5-py3-none-any.whl (204 kB)
   ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 0.0/204.4 kB ? eta -:--:--
   ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━╺━ 194.6/204.4 kB 6.4 MB/s eta 0:00:01
   ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 204.4/204.4 kB 5.5 MB/s eta 0:00:00
Downloading pandas_flavor-0.6.0-py3-none-any.whl (7.2 kB)
Installing collected packages: pandas-flavor, pingouin
Successfully installed pandas-flavor-0.6.0 pingouin-0.5.5

[notice] A new release of pip is available: 24.0 -> 24.3.1
[notice] To update, run: pip install --upgrade pip
import pandas as pd
import numpy as np
import pingouin as pg
import plotly.express as px
import plotly.graph_objects as go
from scipy.stats import ttest_ind
from sklearn.preprocessing import StandardScaler
from sklearn.decomposition import PCA
from sklearn.linear_model import LinearRegression
from sklearn.metrics import mean_squared_error
import statsmodels.api as sm
from statsmodels.formula.api import ols

# Reload the SPSS file
data = pd.read_spss('BM6043 - Data 2024 (1).sav')

# Descriptive statistics
summary_stats = data.describe(include='all')
summary_stats

Descriptive Statistics and Reliability Analysis

Descriptive Statistics

The dataset contains various demographic and behavioral variables. Here are some key statistics:
  • Gender: 286 respondents, with 149 males and 137 females.
  • Age: Mean age is approximately 35.6 years, with a standard deviation of 11.5 years.
  • Ethnicity: Predominantly Chinese (168 respondents).
  • Education: Most respondents have 'SPM or lower' education level (144 respondents).
  • Occupation: The most common occupation is 'Supervisor, Officer or Executive' (77 respondents).

Reliability Analysis

To assess the reliability and internal consistency of the constructs, we will calculate Cronbach's Alpha for each construct. This will help determine the reliability of the set of items for each construct.
# Calculate Cronbach's Alpha for constructs
constructs = {
'Subjective Knowledge': ['SK1', 'SK2', 'SK3', 'SK4', 'SK5'],
'Behavioural Intention': ['Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6'],
'Family Influence': ['Fam1', 'Fam2', 'Fam3'],
'Peer Influence': ['Peer1', 'Peer2', 'Peer3']
}

cronbach_alpha_results = {}

for construct, items in constructs.items():
# Convert categorical responses to numeric
numeric_data = data[items].apply(lambda x: pd.factorize(x)[0])
# Calculate Cronbach's Alpha
alpha = pg.cronbach_alpha(numeric_data)
cronbach_alpha_results[construct] = alpha

cronbach_alpha_results

Reliability Analysis Results

The reliability and internal consistency of the constructs were assessed using Cronbach's Alpha:
  • Subjective Knowledge:
    • Cronbach's Alpha: 0.659
    • Confidence Interval: [0.593, 0.717]
    • Interpretation: Acceptable reliability.
  • Behavioural Intention:
    • Cronbach's Alpha: 0.745
    • Confidence Interval: [0.696, 0.788]
    • Interpretation: Good reliability.
  • Family Influence:
    • Cronbach's Alpha: 0.819
    • Confidence Interval: [0.779, 0.852]
    • Interpretation: Excellent reliability.
  • Peer Influence:
    • Cronbach's Alpha: 0.716
    • Confidence Interval: [0.654, 0.768]
    • Interpretation: Good reliability.
These results indicate that the constructs have acceptable to excellent internal consistency, making them reliable

first of all, give me instructions or how to compute into the spss?

Instructions for Computing Reliability in SPSS

To compute Cronbach's Alpha for reliability analysis in SPSS, follow these steps:
  1. Open SPSS:
    • Launch SPSS and open your dataset.
  2. Navigate to Reliability Analysis:
    • Click on Analyze in the top menu.
    • Select Scale and then Reliability Analysis.
  3. Select Items for Analysis:
    • In the dialog box, move the items (variables) you want to analyze into the Items box. For example, for the "Subjective Knowledge" construct, select SK1, SK2, SK3, SK4, and SK5.
  4. Set the Model:
    • Ensure that the Model is set to Alpha (which is the default).
  5. Options:
    • Click on Statistics and select Scale if item deleted, Inter-Item Correlations, and Summary statistics.
    • Click Continue to return to the main dialog.
  6. Run the Analysis:
    • Click OK to run the analysis.
  7. Interpret the Output:
    • SPSS will provide an output window with the Cronbach's Alpha value and other statistics.
    • A Cronbach's Alpha value above 0.7 generally indicates acceptable reliability, although values above 0.6 can be considered acceptable in exploratory research.
These steps will help you compute and interpret the reliability of your constructs using SPSS. If you have any further questions or need additional assistance, feel free to ask!

BM6043 - Data 2024 (2).sav

Missing file

Descriptive statistics for respondents with the appropriate charts. To determine the

reliability and the internal consistency of the set of items for each construct.
To test the reliability, or the internal consistency of the set of items for each construct.
I get the answer in SPSS is equal to 0.930. is it correct fro, your site?
import pandas as pd

# Load the SPSS file again to ensure we have the latest data
data = pd.read_spss('BM6043 - Data 2024 (2).sav')

# Calculate Cronbach's Alpha for all constructs combined to verify the SPSS result
all_items = ['Fam1', 'Fam2', 'Fam3', 'Peer1', 'Peer2', 'Peer3', 'SK1', 'SK2', 'SK3', 'SK4', 'SK5', 'Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6']

# Convert categorical responses to numeric
numeric_data = data[all_items].apply(lambda x: pd.factorize(x)[0])

# Calculate Cronbach's Alpha for all items
alpha_all = pg.cronbach_alpha(numeric_data)
alpha_all
import pingouin as pg

# Calculate Cronbach's Alpha for all constructs combined to verify the SPSS result
all_items = ['Fam1', 'Fam2', 'Fam3', 'Peer1', 'Peer2', 'Peer3', 'SK1', 'SK2', 'SK3', 'SK4', 'SK5', 'Beh_Int1', 'Beh_Int2', 'Beh_Int3', 'Beh_Int4', 'Beh_Int5', 'Beh_Int6']

# Convert categorical responses to numeric
numeric_data = data[all_items].apply(lambda x: pd.factorize(x)[0])

# Calculate Cronbach's Alpha for all items
alpha_all = pg.cronbach_alpha(numeric_data)
alpha_all
It seems that the pingouin library, which is used to calculate Cronbach's Alpha, is not installed in the current environment. This is causing the ModuleNotFoundError.
To resolve this issue, you can install the pingouin library by running the following command in your Python environment:
Once installed, you can re-run the code to calculate Cronbach's Alpha for the constructs. If you need further assistance, feel free to ask!
Get started with Vizly