ich habe einen Datensatz (agg_data) mit geschätzten Parametern eines Linear Ballistic Accumulator Models (LBA).
cond1 cond2 s vs b t0 A
...
Aus diesen Parametern würde ich gerne Verteilungen generieren. Das Package mit welchem ich arbeite (glba) bietet hierfür folgende Funktion an:
Code: Alles auswählen
rlba(n, b, A, vs, s, t0, st0 = 0, truncdrifts = 0)
Hieraus ergibt sich ein neuer Datensatz mit zwei Variablen ("rt", "resp").
Was ich jetzt gerne machen würde, ist, die rlba-Funktion für jede Reihe meines agg_data-Datensatzes anwenden und die sich daraus ergebenden Datensätze in eigenen files speichern.
Was ich bisher versucht habe, ist folgendes:
Code: Alles auswählen
for (row in 1:nrow(agg_data))
{tmp<-rlba(1000, b =(agg_data[,3]+agg_data[,4]), A = agg_data[,4], vs = c(agg_data[,5], 1agg_data[,5]), s=agg_data[,6], t0 =agg_data[,7], st0 = 0, truncdrifts = TRUE)
myfile <- file.path("myfilepath", paste0(cond1, "_", cond2, ".dat"))
write.table(tmp, file = myfile, sep = "\t", row.names = FALSE, col.names = TRUE,
quote = FALSE, append = FALSE)
rm(tmp, myfile)
}
Besten Dank schon mal!
Viele Grüße
Frida