if(!requireNamespace('lattice', quietly = TRUE)){ install.packages('lattice', dependencies=TRUE, repos='http://cran.rstudio.com/') } library(lattice) ## This didn't work due to missing dependencies: ## ## if(!requireNamespace('RVAideMemoire', quietly = TRUE)){ ## install.packages('RVAideMemoire', dependencies=TRUE, repos='http://cran.rstudio.com/') ## } ## require('RVAideMemoire') ## ## Download the source code and install by hand: ## ## https://cran.r-project.org/web/packages/RVAideMemoire/index.html ## ## Source the code you will actuall use: ## if( file.exists('~/Software/RVAideMemoire/R/mood.medtest.R') ){ source('~/Software/RVAideMemoire/R/mood.medtest.R') } source('../../Data/get_visual_fusion_data.R') fusion = get_visual_fusion_data() ##qv = 0.9; fusion[fusion$Time>quantile(fusion$Time, qv),]$Time = fusion[fusion$Time>quantile(fusion$Time, qv),]$Time * 5.0 ## increase the value of the larger samples print(fusion) fusion$Time = jitter(fusion$Time) histogram(~ Time | Treatment, layout=c(1, 2), data=fusion) ## Compare the means with a t-test: ## nv_samples = fusion$Time[ fusion$Treatment =='NV' ] vv_samples = fusion$Time[ fusion$Treatment =='VV' ] res = t.test(nv_samples, vv_samples) print('T-test:') print(res) ## Use a nonparametric method: ## res = mood.medtest(Time ~ Treatment, data=fusion) print('Median Test:') print(res) ## Use the rank-sum test: ## res = wilcox.test(Time ~ Treatment, data=fusion) print('Wilcox Test:') print(res)