lines = '3.84 3.54 3.96 3.58 3.50 3.34 3.76 3.74 3.56 3.58 3.70 3.64 3.54 3.72 3.58 4.04 3.42 3.86 3.94 4.04' con = textConnection( lines ) avg = scan( con, quiet=TRUE ) close(con) lines = '0.6 0.4 1.1 0.3 0.9 0.6 0.5 0.9 0.5 1.0 0.9 0.1 0.6 0.6 1.1 0.6 0.9 1.3 0.5 1.4' con = textConnection( lines ) range = scan( con, quiet=TRUE ) close(con) DF = data.frame( average=avg, range=range ) sample_1 = c( 4.0, 4.2, 3.6, 3.8, 3.6 ) sample_2 = c( 3.3, 3.7, 3.4, 3.6, 3.7 ) sample_20 = c( 4.2, 4.9, 3.6, 3.5, 4.0 ) lines = '3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9' con = textConnection( lines ) x_hist_mid = scan( con, quiet=TRUE ) close(con) lines = '5 3 4 5 6 10 13 11 9 10 11 2 5 4 1 0 0 0 0 1' con = textConnection( lines ) x_hist_cnt = scan( con, quiet=TRUE ) close(con) # From the histogram data lets construct the raw data: # raw_data = c() for( ii in 1:length(x_hist_cnt) ){ raw_data = c( raw_data, rep( x_hist_mid[ii], x_hist_cnt[ii] ) ) } x_breaks = c( 3, c( x_hist_mid[2:length(x_hist_mid)] + x_hist_mid[1:(length(x_hist_mid)-1)] )/2, 5.1 )