Appendix
Script for blinded estimation of the standard deviation of the effectiveness of UE-FMA/66 at T1 by the intervention group in R language.# Load the dataset (3 columns: record id, endpoint value, treatment).load(“db.Rdata”).data = db[,c(“record_id”,“endpoint”,“prot”)].# Create a group variable (empty vector).data$group = 0.# Randomly sample one unit from the population.rand = sample(data$record_id,1).# Assign group 1 to this unit.data$group[data$record_id =  = rand] = 1.# Group 1 is assigned to all units that come from the same treatment as the selected patient.data$group = ifelse(data$prot =  = data$prot[data$record_id =  = rand],1,0).# Remove the treatment column from the dataset.data = data[,-3].# Compute the standard deviation of the endpoint (stratified by intervention group).sd = aggregate(data$endpoint, by = list(data$group), FUN = sd).# Print the results.sd.