Skip to content

Instantly share code, notes, and snippets.

@agricolamz

agricolamz/code_for_6_task.R

Last active February 18, 2025 09:23
Show Gist options
  • Save agricolamz/762d19a17e3e165efdde25758135c749 to your computer and use it in GitHub Desktop.
Save agricolamz/762d19a17e3e165efdde25758135c749 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
code_for_6_task.R
library(tidyverse)
read_csv("https://raw.githubusercontent.com/agricolamz/2024_HSE_b_da4l/master/data/Coretta_2017_icelandic.csv") |>
filter(speaker == "tt01") ->
vowels
mean_data <- mean(vowels$vowel.dur)
sd_data <- sd(vowels$vowel.dur)
m1 <- function(x) dnorm(x, mean = mean_data, sd = sd_data) *
dnorm(x, mean = 87, sd = 25)
m2 <- function(x) dnorm(x, mean = mean_data, sd = sd_data) *
dnorm(x, mean = 85, sd = 30)
integrate(m1, -Inf, Inf)[[1]]/integrate(m2, -Inf, Inf)[[1]]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment