Question 1

data(cars)
summary(cars)

Question 2

install.packages("jsonlite")
library(jsonlite)
BTC <- fromJSON("https://min-api.cryptocompare.com/data/v2/histoday?fsym=BTC&tsym=USD&limit=100")
str(BTC)

BTC_data <- BTC$Data$Data  

head(BTC_data)

max_close_price <- max(BTC_data$close, na.rm = TRUE)

max_close_price

Question 3

Identify a topic of interest and give your project a name/title.

Community health in relation to patient satisfaction and quality of life

Phrase 3-5 research questions you would like to explore.

1. How does health service type (rehab, preventive, consultation) relate to patient satisfaction and quality of life scores?

2. What is the relationship between visit frequency and patient satisfaction?

3. What biomechanical characteristics are most strongly associated with high EMG activity levels?

List the data sources that you find that are relevant with your research questions.

Extract one or more relevant datasets associate with your research questions, either import the downloaded dataset(s), extract from APIs, or ethically scrape the web, etc.

library(readr)
community_health_data <- read_csv("~/Downloads/community_health_evaluation_dataset.csv")
Rows: 347 Columns: 12── Column specification ─────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): Gender, Service Type, Visit Frequency, EMG Activity
dbl (8): Participant ID, Age, SES, Step Frequency (steps/min), Stride Length (m), Joint Angle (°)...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
View(community_health_data)

Describe your data extracted, statistically and/or visually.

The Community Health Evaluation Dataset contains synthesized data on 347 participants, capturing demographic details, healthcare usage patterns, and biomechanical metrics to assess community health services and quality. It includes variables such as age, socioeconomic status, visit frequency, stride length, joint angles, and patient satisfaction, aimed at understanding relationships between biomechanical health indicators and quality of life.

colnames(community_health_data) <- c("Participant_ID", "Age", "Gender", "SES", "Service_Type",
                                     "Visit_Frequency", "Step_Frequency", "Stride_Length",
                                     "Joint_Angle", "EMG_Activity", "Patient_Satisfaction",
                                     "Quality_of_Life_Score")

library(ggplot2)
ggplot(community_health_data, aes(x = Visit_Frequency, y = Patient_Satisfaction, fill = Visit_Frequency)) +
  geom_boxplot(outlier.color = "red", outlier.shape = 16, outlier.size = 2) +
  labs(title = "Relationship between Visit Frequency and Patient Satisfaction",
       x = "Visit Frequency",
       y = "Patient Satisfaction") +
  theme_minimal() +
  theme(legend.position = "none")

Key observations:

Weekly visits have the highest median satisfaction score, as represented by the middle line in the box, which is above 7.

Yearly and Monthly visits have lower median satisfaction scores, hovering around 5, with the Monthly group showing slightly more variation in satisfaction.

The spread (interquartile range) of satisfaction scores is largest for the “Weekly” group, indicating that weekly visitors reported a wider range of satisfaction.

Outliers are not shown in this plot, but the whiskers indicate the range of the data within 1.5 times the interquartile range from the quartiles.

library(ggplot2)
ggplot(community_health_data, aes(x = Service_Type, y = Patient_Satisfaction, fill = Service_Type)) +
  geom_boxplot(outlier.color = "red", outlier.shape = 16, outlier.size = 2) + # Outliers are marked in red
  labs(title = "Patient Satisfaction by Health Service Type",
       x = "Health Service Type",
       y = "Patient Satisfaction (1-10)") +
  theme_minimal() +
  theme(legend.position = "none")

Key Observations:

Median Satisfaction Scores: All three health service types (Consultation, Preventive, and Rehab) have similar median patient satisfaction scores, centered around 5. This suggests that the typical satisfaction level is consistent across these service types.

Variation in Satisfaction: The Consultation and Rehab services have wider interquartile ranges (IQRs) compared to Preventive services, indicating that patient satisfaction is more varied in these two categories. Specifically, patients receiving Consultation or Rehab services reported a broader range of satisfaction levels.

Outliers and Extremes: None of the service types show extreme outliers in patient satisfaction, but the Consultation service has the lowest satisfaction score (around 2.5), while the Preventive and Rehab services have a slightly more concentrated distribution around the median. This could suggest that Consultation services may have more room for improvement in terms of patient satisfaction.

Perform necessary data cleaning and manipulation especially if the raw data contains special values or not directly in the format that can answer your research questions.

missing_values <- colSums(is.na(community_health_data))

missing_values <- missing_values[missing_values > 0]

print("Columns with missing values and their counts:")
print(missing_values)

no missing values

List future data preparation work needed if any.

continuing to segment the data to make correlations more visable.

LS0tCnRpdGxlOiAiTWluaSBQcm9qZWN0IgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgojIFF1ZXN0aW9uIDEKYGBge3J9CmRhdGEoY2FycykKc3VtbWFyeShjYXJzKQpgYGAKCiMgUXVlc3Rpb24gMgpgYGB7cn0KaW5zdGFsbC5wYWNrYWdlcygianNvbmxpdGUiKQpsaWJyYXJ5KGpzb25saXRlKQpCVEMgPC0gZnJvbUpTT04oImh0dHBzOi8vbWluLWFwaS5jcnlwdG9jb21wYXJlLmNvbS9kYXRhL3YyL2hpc3RvZGF5P2ZzeW09QlRDJnRzeW09VVNEJmxpbWl0PTEwMCIpCmBgYAoKYGBge3J9CnN0cihCVEMpCgpCVENfZGF0YSA8LSBCVEMkRGF0YSREYXRhICAKCmhlYWQoQlRDX2RhdGEpCgptYXhfY2xvc2VfcHJpY2UgPC0gbWF4KEJUQ19kYXRhJGNsb3NlLCBuYS5ybSA9IFRSVUUpCgptYXhfY2xvc2VfcHJpY2UKYGBgCgojIFF1ZXN0aW9uIDMKCiMjIElkZW50aWZ5IGEgdG9waWMgb2YgaW50ZXJlc3QgYW5kIGdpdmUgeW91ciBwcm9qZWN0IGEgbmFtZS90aXRsZS4KCiMjIyMgQ29tbXVuaXR5IGhlYWx0aCBpbiByZWxhdGlvbiB0byBwYXRpZW50IHNhdGlzZmFjdGlvbiBhbmQgcXVhbGl0eSBvZiBsaWZlCgoKCiMjIFBocmFzZSAzLTUgcmVzZWFyY2ggcXVlc3Rpb25zIHlvdSB3b3VsZCBsaWtlIHRvIGV4cGxvcmUuCgojIyMjIDEuIEhvdyBkb2VzIGhlYWx0aCBzZXJ2aWNlIHR5cGUgKHJlaGFiLCBwcmV2ZW50aXZlLCBjb25zdWx0YXRpb24pIHJlbGF0ZSB0byBwYXRpZW50IHNhdGlzZmFjdGlvbiBhbmQgcXVhbGl0eSBvZiBsaWZlIHNjb3Jlcz8KCiMjIyMgMi4gV2hhdCBpcyB0aGUgcmVsYXRpb25zaGlwIGJldHdlZW4gdmlzaXQgZnJlcXVlbmN5IGFuZCBwYXRpZW50IHNhdGlzZmFjdGlvbj8KCiMjIyMgMy4gV2hhdCBiaW9tZWNoYW5pY2FsIGNoYXJhY3RlcmlzdGljcyBhcmUgbW9zdCBzdHJvbmdseSBhc3NvY2lhdGVkIHdpdGggaGlnaCBFTUcgYWN0aXZpdHkgbGV2ZWxzPwoKCgojIyBMaXN0IHRoZSBkYXRhIHNvdXJjZXMgdGhhdCB5b3UgZmluZCB0aGF0IGFyZSByZWxldmFudCB3aXRoIHlvdXIgcmVzZWFyY2ggcXVlc3Rpb25zLgojIyMjIGh0dHBzOi8vd3d3LmthZ2dsZS5jb20vZGF0YXNldHMveml5YTA3L2NvbW11bml0eS1oZWFsdGgtZXZhbHVhdGlvbi1kYXRhc2V0CgoKCiMjIEV4dHJhY3Qgb25lIG9yIG1vcmUgcmVsZXZhbnQgZGF0YXNldHMgYXNzb2NpYXRlIHdpdGggeW91ciByZXNlYXJjaCBxdWVzdGlvbnMsIGVpdGhlciBpbXBvcnQgdGhlIGRvd25sb2FkZWQgZGF0YXNldChzKSwgZXh0cmFjdCBmcm9tIEFQSXMsIG9yIGV0aGljYWxseSBzY3JhcGUgdGhlIHdlYiwgZXRjLgpgYGB7cn0KbGlicmFyeShyZWFkcikKY29tbXVuaXR5X2hlYWx0aF9kYXRhIDwtIHJlYWRfY3N2KCJ+L0Rvd25sb2Fkcy9jb21tdW5pdHlfaGVhbHRoX2V2YWx1YXRpb25fZGF0YXNldC5jc3YiKQpWaWV3KGNvbW11bml0eV9oZWFsdGhfZGF0YSkKYGBgCgoKCiMjIERlc2NyaWJlIHlvdXIgZGF0YSBleHRyYWN0ZWQsIHN0YXRpc3RpY2FsbHkgYW5kL29yIHZpc3VhbGx5LgoKIyMjIyBUaGUgQ29tbXVuaXR5IEhlYWx0aCBFdmFsdWF0aW9uIERhdGFzZXQgY29udGFpbnMgc3ludGhlc2l6ZWQgZGF0YSBvbiAzNDcgcGFydGljaXBhbnRzLCBjYXB0dXJpbmcgZGVtb2dyYXBoaWMgZGV0YWlscywgaGVhbHRoY2FyZSB1c2FnZSBwYXR0ZXJucywgYW5kIGJpb21lY2hhbmljYWwgbWV0cmljcyB0byBhc3Nlc3MgY29tbXVuaXR5IGhlYWx0aCBzZXJ2aWNlcyBhbmQgcXVhbGl0eS4gSXQgaW5jbHVkZXMgdmFyaWFibGVzIHN1Y2ggYXMgYWdlLCBzb2Npb2Vjb25vbWljIHN0YXR1cywgdmlzaXQgZnJlcXVlbmN5LCBzdHJpZGUgbGVuZ3RoLCBqb2ludCBhbmdsZXMsIGFuZCBwYXRpZW50IHNhdGlzZmFjdGlvbiwgYWltZWQgYXQgdW5kZXJzdGFuZGluZyByZWxhdGlvbnNoaXBzIGJldHdlZW4gYmlvbWVjaGFuaWNhbCBoZWFsdGggaW5kaWNhdG9ycyBhbmQgcXVhbGl0eSBvZiBsaWZlLgoKYGBge3J9CmNvbG5hbWVzKGNvbW11bml0eV9oZWFsdGhfZGF0YSkgPC0gYygiUGFydGljaXBhbnRfSUQiLCAiQWdlIiwgIkdlbmRlciIsICJTRVMiLCAiU2VydmljZV9UeXBlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJWaXNpdF9GcmVxdWVuY3kiLCAiU3RlcF9GcmVxdWVuY3kiLCAiU3RyaWRlX0xlbmd0aCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSm9pbnRfQW5nbGUiLCAiRU1HX0FjdGl2aXR5IiwgIlBhdGllbnRfU2F0aXNmYWN0aW9uIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJRdWFsaXR5X29mX0xpZmVfU2NvcmUiKQoKbGlicmFyeShnZ3Bsb3QyKQpnZ3Bsb3QoY29tbXVuaXR5X2hlYWx0aF9kYXRhLCBhZXMoeCA9IFZpc2l0X0ZyZXF1ZW5jeSwgeSA9IFBhdGllbnRfU2F0aXNmYWN0aW9uLCBmaWxsID0gVmlzaXRfRnJlcXVlbmN5KSkgKwogIGdlb21fYm94cGxvdChvdXRsaWVyLmNvbG9yID0gInJlZCIsIG91dGxpZXIuc2hhcGUgPSAxNiwgb3V0bGllci5zaXplID0gMikgKwogIGxhYnModGl0bGUgPSAiUmVsYXRpb25zaGlwIGJldHdlZW4gVmlzaXQgRnJlcXVlbmN5IGFuZCBQYXRpZW50IFNhdGlzZmFjdGlvbiIsCiAgICAgICB4ID0gIlZpc2l0IEZyZXF1ZW5jeSIsCiAgICAgICB5ID0gIlBhdGllbnQgU2F0aXNmYWN0aW9uIikgKwogIHRoZW1lX21pbmltYWwoKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKQpgYGAKIyMgS2V5IG9ic2VydmF0aW9uczoKIyMjIyBXZWVrbHkgdmlzaXRzIGhhdmUgdGhlIGhpZ2hlc3QgbWVkaWFuIHNhdGlzZmFjdGlvbiBzY29yZSwgYXMgcmVwcmVzZW50ZWQgYnkgdGhlIG1pZGRsZSBsaW5lIGluIHRoZSBib3gsIHdoaWNoIGlzIGFib3ZlIDcuCgojIyMjIFllYXJseSBhbmQgTW9udGhseSB2aXNpdHMgaGF2ZSBsb3dlciBtZWRpYW4gc2F0aXNmYWN0aW9uIHNjb3JlcywgaG92ZXJpbmcgYXJvdW5kIDUsIHdpdGggdGhlIE1vbnRobHkgZ3JvdXAgc2hvd2luZyBzbGlnaHRseSBtb3JlIHZhcmlhdGlvbiBpbiBzYXRpc2ZhY3Rpb24uCgojIyMjIFRoZSBzcHJlYWQgKGludGVycXVhcnRpbGUgcmFuZ2UpIG9mIHNhdGlzZmFjdGlvbiBzY29yZXMgaXMgbGFyZ2VzdCBmb3IgdGhlICJXZWVrbHkiIGdyb3VwLCBpbmRpY2F0aW5nIHRoYXQgd2Vla2x5IHZpc2l0b3JzIHJlcG9ydGVkIGEgd2lkZXIgcmFuZ2Ugb2Ygc2F0aXNmYWN0aW9uLgoKIyMjIyBPdXRsaWVycyBhcmUgbm90IHNob3duIGluIHRoaXMgcGxvdCwgYnV0IHRoZSB3aGlza2VycyBpbmRpY2F0ZSB0aGUgcmFuZ2Ugb2YgdGhlIGRhdGEgd2l0aGluIDEuNSB0aW1lcyB0aGUgaW50ZXJxdWFydGlsZSByYW5nZSBmcm9tIHRoZSBxdWFydGlsZXMuCgpgYGB7cn0KbGlicmFyeShnZ3Bsb3QyKQpnZ3Bsb3QoY29tbXVuaXR5X2hlYWx0aF9kYXRhLCBhZXMoeCA9IFNlcnZpY2VfVHlwZSwgeSA9IFBhdGllbnRfU2F0aXNmYWN0aW9uLCBmaWxsID0gU2VydmljZV9UeXBlKSkgKwogIGdlb21fYm94cGxvdChvdXRsaWVyLmNvbG9yID0gInJlZCIsIG91dGxpZXIuc2hhcGUgPSAxNiwgb3V0bGllci5zaXplID0gMikgKyAjIE91dGxpZXJzIGFyZSBtYXJrZWQgaW4gcmVkCiAgbGFicyh0aXRsZSA9ICJQYXRpZW50IFNhdGlzZmFjdGlvbiBieSBIZWFsdGggU2VydmljZSBUeXBlIiwKICAgICAgIHggPSAiSGVhbHRoIFNlcnZpY2UgVHlwZSIsCiAgICAgICB5ID0gIlBhdGllbnQgU2F0aXNmYWN0aW9uICgxLTEwKSIpICsKICB0aGVtZV9taW5pbWFsKCkgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikKCmBgYAoKIyMgS2V5IE9ic2VydmF0aW9uczoKCiMjIyMgTWVkaWFuIFNhdGlzZmFjdGlvbiBTY29yZXM6IEFsbCB0aHJlZSBoZWFsdGggc2VydmljZSB0eXBlcyAoQ29uc3VsdGF0aW9uLCBQcmV2ZW50aXZlLCBhbmQgUmVoYWIpIGhhdmUgc2ltaWxhciBtZWRpYW4gcGF0aWVudCBzYXRpc2ZhY3Rpb24gc2NvcmVzLCBjZW50ZXJlZCBhcm91bmQgNS4gVGhpcyBzdWdnZXN0cyB0aGF0IHRoZSB0eXBpY2FsIHNhdGlzZmFjdGlvbiBsZXZlbCBpcyBjb25zaXN0ZW50IGFjcm9zcyB0aGVzZSBzZXJ2aWNlIHR5cGVzLgoKIyMjIyBWYXJpYXRpb24gaW4gU2F0aXNmYWN0aW9uOiBUaGUgQ29uc3VsdGF0aW9uIGFuZCBSZWhhYiBzZXJ2aWNlcyBoYXZlIHdpZGVyIGludGVycXVhcnRpbGUgcmFuZ2VzIChJUVJzKSBjb21wYXJlZCB0byBQcmV2ZW50aXZlIHNlcnZpY2VzLCBpbmRpY2F0aW5nIHRoYXQgcGF0aWVudCBzYXRpc2ZhY3Rpb24gaXMgbW9yZSB2YXJpZWQgaW4gdGhlc2UgdHdvIGNhdGVnb3JpZXMuIFNwZWNpZmljYWxseSwgcGF0aWVudHMgcmVjZWl2aW5nIENvbnN1bHRhdGlvbiBvciBSZWhhYiBzZXJ2aWNlcyByZXBvcnRlZCBhIGJyb2FkZXIgcmFuZ2Ugb2Ygc2F0aXNmYWN0aW9uIGxldmVscy4KCiMjIyMgT3V0bGllcnMgYW5kIEV4dHJlbWVzOiBOb25lIG9mIHRoZSBzZXJ2aWNlIHR5cGVzIHNob3cgZXh0cmVtZSBvdXRsaWVycyBpbiBwYXRpZW50IHNhdGlzZmFjdGlvbiwgYnV0IHRoZSBDb25zdWx0YXRpb24gc2VydmljZSBoYXMgdGhlIGxvd2VzdCBzYXRpc2ZhY3Rpb24gc2NvcmUgKGFyb3VuZCAyLjUpLCB3aGlsZSB0aGUgUHJldmVudGl2ZSBhbmQgUmVoYWIgc2VydmljZXMgaGF2ZSBhIHNsaWdodGx5IG1vcmUgY29uY2VudHJhdGVkIGRpc3RyaWJ1dGlvbiBhcm91bmQgdGhlIG1lZGlhbi4gVGhpcyBjb3VsZCBzdWdnZXN0IHRoYXQgQ29uc3VsdGF0aW9uIHNlcnZpY2VzIG1heSBoYXZlIG1vcmUgcm9vbSBmb3IgaW1wcm92ZW1lbnQgaW4gdGVybXMgb2YgcGF0aWVudCBzYXRpc2ZhY3Rpb24uCgoKCiMjIFBlcmZvcm0gbmVjZXNzYXJ5IGRhdGEgY2xlYW5pbmcgYW5kIG1hbmlwdWxhdGlvbiBlc3BlY2lhbGx5IGlmIHRoZSByYXcgZGF0YSBjb250YWlucyBzcGVjaWFsIHZhbHVlcyBvciBub3QgZGlyZWN0bHkgaW4gdGhlIGZvcm1hdCB0aGF0IGNhbiBhbnN3ZXIgeW91ciByZXNlYXJjaCBxdWVzdGlvbnMuCmBgYHtyfQptaXNzaW5nX3ZhbHVlcyA8LSBjb2xTdW1zKGlzLm5hKGNvbW11bml0eV9oZWFsdGhfZGF0YSkpCgptaXNzaW5nX3ZhbHVlcyA8LSBtaXNzaW5nX3ZhbHVlc1ttaXNzaW5nX3ZhbHVlcyA+IDBdCgpwcmludCgiQ29sdW1ucyB3aXRoIG1pc3NpbmcgdmFsdWVzIGFuZCB0aGVpciBjb3VudHM6IikKcHJpbnQobWlzc2luZ192YWx1ZXMpCmBgYAojIyMjIG5vIG1pc3NpbmcgdmFsdWVzCgojIyBMaXN0IGZ1dHVyZSBkYXRhIHByZXBhcmF0aW9uIHdvcmsgbmVlZGVkIGlmIGFueS4KIyMjIyBjb250aW51aW5nIHRvIHNlZ21lbnQgdGhlIGRhdGEgdG8gbWFrZSBjb3JyZWxhdGlvbnMgbW9yZSB2aXNhYmxlLgo=