﻿# Experiment 4: 2022/23 Replication of Field-Dependency and Cybersickness
# Study run by J. Josupeit, L. Greim, and S. Sanchez Rivas


# this file contains adapted code for 4 special cases regarding the R-script "VR_RFT_Preprocessing_RFT.R".
This affected the participants 54, 10, 72, and 1 to 5, which are listed in R as Y <- 1, 40, 79 and 76, 77, 78, 13, 14.
Adapt the code in "VR_RFT_Preprocessing_RFT.R" for the each special case as follows:


# for participant 54 (Y ← 1) replace the lines 140-189 with:

Unity_Missing <- Unity_New[Unity_New$Trialnumber == '1',]
Unity_Missing['Z_Rod'][Unity_Missing['Z_Rod'] == 22.1000] <- 22
Unity_New <- rbind(Unity_New, Unity_Missing)
Unity_New <- Unity_New %>%
  mutate(
    Constel = case_when(
      (Z_Rod == '-22' & Z_Frame == '-33') ~ "-22_-33",
      (Z_Rod == '-22' & Z_Frame == '33') ~ "-22_33",
      (Z_Rod == '-22' & Z_Frame == '0') ~ "-22_0",
      (Z_Rod == '22' & Z_Frame == '-33') ~ "22_-33",
      (Z_Rod == '22' & Z_Frame == '33') ~ "22_33",
      (Z_Rod == '22' & Z_Frame == '0') ~ "22_0",
      (Z_Rod == '-11' & Z_Frame == '-33') ~ "-11_-33",
      (Z_Rod == '-11' & Z_Frame == '33') ~ "-11_33",
      (Z_Rod == '-11' & Z_Frame == '0') ~ "-11_0",
      (Z_Rod == '11' & Z_Frame == '-33') ~ "11_-33",
      (Z_Rod == '11' & Z_Frame == '33') ~ "11_33",
      (Z_Rod == '11' & Z_Frame == '0') ~ "11_0"
    )
  )

Unity_New_noNA <- na.omit(Unity_New)
Unity_New_noNA <- Unity_New_noNA[, c('Constel', 'Trialnumber')]
Unity_New_NA <- Unity_New[rowSums(is.na(Unity_New)) > 0,]
Unity_New_NA <-
  merge(Unity_New_noNA, Unity_New_NA, by = 'Trialnumber', all = TRUE)
Unity_New_NA <- Unity_New_NA [, c(-9)]
Unity_New_NA <- rename(Unity_New_NA, c('Constel' = 'Constel.x'))
Unity_New_NA <- Unity_New_NA [c(-2),]
Unity_New_NA <- Unity_New_NA [c(-1),]
Unity_New_NA <- Unity_New_NA %>%
  mutate(Condi = case_when(
    (
      Constel == "22_-33" |
        Constel == "-22_33" |
        Constel == "11_-33" | Constel == "-11_33"
    ) ~ "Incompatible",
    (
      Constel == "-22_-33" |
        Constel == "22_33" |
        Constel == "-11_-33" | Constel == "11_33"
    ) ~ "Compatible",
    (
      Constel == "22_0" |
        Constel == "-22_0" |
        Constel == "11_0" | Constel == "-11_0"
    ) ~ "Neutral"
  ))


# for participants 1-5 (Y <- 76, 77, 78, 13, 14) replace line 63-71 with:

#first 3 entries from Unity (not filtered)
Unity_Constellation <-
  ddply(Unity, 'Trialnumber', function(group) {
        return(
          data.frame(Z_Rod = head(group[['Z_Rod']], n = 3)) %>%
            data.frame(Z_Frame = head(group[['Z_Frame']], n = 3)) %>%
            data.frame(Time = head(group[['TimestampNew']], n = 3))
        )
  })

#last entry of the first 3
Unity_Constellation <-
  ddply(Unity_Constellation, 'Trialnumber', function(group) {
        return(
          data.frame(Z_Rod = tail(group[['Z_Rod']], n = 1)) %>%
            data.frame(Z_Frame = tail(group[['Z_Frame']], n = 1)) %>%
            data.frame(Time = tail(group[['Time']], n = 1))
    )
  })


# for participant 10 (Y ← 40) replace line 160-189 with:

Unity_Missing <- Unity_New[Unity_New$Trialnumber == '40',]
Unity_New_noNA <- na.omit(Unity_New)
Unity_New_noNA <- Unity_New_noNA[, c('Constel', 'Trialnumber')]
Unity_New_NA <- Unity_New[rowSums(is.na(Unity_New)) > 0,]
last_row <- tail(Unity_Missing, n = 1)
Unity_New_NA <- rbind(Unity_New_NA, last_row)
Unity_New_NA <-
  merge(Unity_New_noNA, Unity_New_NA, by = 'Trialnumber', all = TRUE)
Unity_New_NA <- Unity_New_NA [, c(-9)]
Unity_New_NA <- rename(Unity_New_NA, c('Constel' = 'Constel.x'))
Unity_New_NA <- Unity_New_NA %>%
  mutate(Condi = case_when((
    Constel == "22_-33" |
      Constel == "-22_33" |
      Constel == "11_-33" | Constel == "-11_33"
  ) ~ "Incompatible",
  (
    Constel == "-22_-33" |
      Constel == "22_33" |
      Constel == "-11_-33" | Constel == "11_33"
  ) ~ "Compatible",
  (
    Constel == "22_0" |
      Constel == "-22_0" |
      Constel == "11_0" | Constel == "-11_0"
  ) ~ "Neutral"
  ))
Unity_New_NA <- Unity_New_NA[-c(39), ]


# for participant 72 (Y ← 79) replace line 160-189 with:

Unity_Missing <- Unity_New[Unity_New$Trialnumber == '13',]
Unity_New_noNA <- na.omit(Unity_New)
Unity_New_noNA <- Unity_New_noNA[, c('Constel', 'Trialnumber')]
Unity_New_NA <- Unity_New[rowSums(is.na(Unity_New)) > 0,]
Unity_New_NA <- rbind(Unity_New_NA, Unity_Missing)
Unity_New_NA <-
  merge(Unity_New_noNA, Unity_New_NA, by = 'Trialnumber', all = TRUE)
Unity_New_NA <- Unity_New_NA [, c(-9)]
Unity_New_NA <- rename(Unity_New_NA, c('Constel' = 'Constel.x'))
Unity_New_NA <- Unity_New_NA %>%
  mutate(Condi = case_when((
    Constel == "22_-33" |
      Constel == "-22_33" |
      Constel == "11_-33" | Constel == "-11_33"
  ) ~ "Incompatible",
  (
    Constel == "-22_-33" |
      Constel == "22_33" |
      Constel == "-11_-33" | Constel == "11_33"
  ) ~ "Compatible",
  (
    Constel == "22_0" |
      Constel == "-22_0" |
      Constel == "11_0" | Constel == "-11_0"
  ) ~ "Neutral"
  ))

 






