library(psych)
library(retimes)
library(gamlss.dist)
source('ICM_Simulations_for_Dissertation_abs.R')
ICM_Fit_Abs <- function(
Observed_Data1 = c(0, -38, -54, 94, 83, 39),
Observed_Data2 = c(0, -39, -34, 75, 82, 5),
Num_Sims = 10000
){
Simulated_Data <- list()
FitData <- list()
hrecord <- vector()
cor1 <- vector()
cor2 <- vector()
chi1 <- vector()
pval1 <- vector()
chi2 <- vector()
pval2 <- vector()
L1Record <- vector()
comprecord <- vector()
RMS1 <- vector()
RMS2 <- vector()
for(i in 1:Num_Sims){
competition = runif(1, 0.5, 0.58)
L1Strength = runif(1, 1, 10)
h = runif(1, 0.001, 10) #Language Activation parameter on switch trials = h*3
x = ICMSimulation_Abs(Simulations = 1, Participants = 1, Num_SubBlocks = 8, Trial_Comparisons = c(1,2,3,4,5,6), Comp = competition, h1 = h, h2 = h, L1_Strength = L1Strength, L2_Strength = L1Strength, NoiseMu = .0001, NoiseTau = .0001, NoiseSigma = .001, p1 = 0.35, c =.01)
Simulated_Data1 <- x$Simulation.Results[seq(1, 11, 2),]$Mean.RT
Simulated_Data1 <- Simulated_Data1 - Simulated_Data1[1]
Simulated_Data2 <- x$Simulation.Results[seq(2, 12, 2),]$Mean.RT
Simulated_Data2 <- Simulated_Data2 - Simulated_Data2[1]
Simulated_Data3 <- data.frame(Simulated_Data1, Simulated_Data2, Simulation = i)
Simulated_Data[[i]] <- Simulated_Data3
chidata1 <- data.frame(Simulated_Data1 +50000 , Observed_Data1 +50000) #Make Sure chi square isn't negative
chidata2 <- data.frame(Simulated_Data2 +50000 , Observed_Data2 + 50000)
cor1[i] <- cor(Simulated_Data1, Observed_Data1)
cor2[i] <- cor(Simulated_Data2, Observed_Data2)
hrecord[i] = h
L1Record[i] = L1Strength
comprecord[i] = competition
chi1[i] <- suppressWarnings(as.numeric(chisq.test(chidata1)[1]))
pval1[i] <- suppressWarnings(as.numeric(chisq.test(chidata1)[3]))
chi2[i] <- suppressWarnings(as.numeric(chisq.test(chidata2)[1]))
pval2[i] <- suppressWarnings(as.numeric(chisq.test(chidata2)[3]))
RMS1[i] <- sqrt(mean((Simulated_Data1 - Observed_Data1)^2))
RMS2[i] <- sqrt(mean((Simulated_Data2 - Observed_Data2)^2))
}
df1 <- data.frame(cor1, cor2, hrecord, L1Record, comprecord, chi1, pval1, chi2, pval2, RMS1, RMS2)
df2 <- do.call("rbind", Simulated_Data)
resultsList <- list(fits = df1, SimData = df2)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.