plot_diff sets to one subject level

[ahsenko]

Hi, I want to use plot_diff function as follows, but not as * subject : factor; set to the value(s): 105798576. How can I compute it for overall differences across subjects?

plot_diff(diff_model, view="timestamp", rm.ranef=F, comp=list(group=c("face", "letter")), cond=list(diff_type=c("i_t_diff")), main="i_t_diff face-letter", ylab = "Difference in pupil dilation (a.u.)", n.grid = length(unique(merge_data$timestamp)), hide.label = TRUE) Summary:

• diff_type : factor; set to the value(s): i_t_diff.
• timestamp : numeric predictor; with 594 values ranging from -750.000000 to 3994.000000.
• subject : factor; set to the value(s): 105798576.

timestamp window(s) of significant difference(s): 1594.000000 - 3098.000000

my model is: diff_model <- bam(diff_values ~ group + diff_type +group:diff_type

• s(timestamp, k = 150, bs = "ad")
• s(timestamp, by=group, bs = "ad", k = 100)
• s(timestamp, by=diff_type, bs = "ad", k = 100)
• s(subject, bs = "re", k = 80) #looks unrequired
• ti(timestamp, group, bs = c("fs", "fs"), k = c(30, 30))
• ti(timestamp, diff_type, group, bs = c("cr", "fs", "fs"), k = c(30, 30, 30))
• ti(timestamp,diff_type, bs = c("fs", "fs"), k = c(30, 30)) #time, group, 3ple
• ti(timestamp, diff_type, by = group, bs = c("cr", "fs"), k = c(30, 30)), data = diff_result, discrete = TRUE, nthreads = 25)

[gaborcsardi]

Hi, this is a read only mirror of CRAN, please see the package authors in the DESCRIPTION file. Look for Maintainer, BugReports and URL. Thanks!