BookmarkSubscribeRSS Feed
mariia8b
Calcite | Level 5

Hello! 

 

SAS newbie here, so I hope someone can help me!

 

A little context: 

 

I did an experiment where I attached nine acceloremeters to a blood turner.
The acceloremeters detects movement and collects data every 60 second as 'Counts per minute' (CPM). The blood turner turned for 10 minutes. That means I have 10 measurements for each accelerometer (9 x 10 measurements in total). The rotation was constant and all the nine accelerometers were attached at the same time. 

 

The aim of the experiment was to examine whether the accelerometers measures the same movement - maybe some kind of reliability test(?). Since the accelerometers were exposed to the same rotation my hypothesis is that there is no difference between the measurements. However, I am unsure about what is the right approach to this including which statistical test I should do? 

 

To make a long story short: I want to examine if the different accelerometers are measuring the same when the exposure, the rotation, is the same.

 

I hope someone can help. I'm really desperate 🙂 

 

Best regards
Maria

 

My data: 

Timeaccelerometer_IDCPM
00:01123871,13
00:02124357,09
00:03123883,17
00:04123843,27
00:05123810,54
00:06123824,35
00:07123827,39
00:08123819,94
00:09123827,7
00:10123834,48
00:01223482,77
00:02223323,05
00:03223514,34
00:04223507,99
00:05223520,11
00:06223495,93
00:07223482,15
00:08223501,74
00:09223487,59
00:10223309,01
00:01324367,35
00:02325030,93
00:03323924,86
00:04323824,59
00:05323805,19
00:06323789,9
00:07323800,91
00:08323782,73
00:09323798,33
00:10323833,54
00:01423788,85
00:02423621,23
00:03423759,49
00:04423810,91
00:05423821,34
00:06423800,15
00:07423802,32
00:08423820,26
00:09423794,8
00:10423808,73
00:01523523,54
00:02523563
00:03523471,06
00:04523504,21
00:05523503,65
00:06523521,73
00:07523506,33
00:08523491,45
00:09523524,15
00:10523523,19
00:01623501,65
00:02623612,43
00:03623531,49
00:04623567,37
00:05623592,66
00:06623587,56
00:07623569,51
00:08623577,08
00:09623569,36
00:10623565,01
00:01723668,53
00:02724571,37
00:03723833,24
00:04723755,22
00:05723765,99
00:06723740,07
00:07723724,57
00:08723760,98
00:09723727,77
00:10723740,25
00:01823789,84
00:02824796,97
00:03823912,49
00:04823818,4
00:05823820,1
00:06823822,9
00:07823804,95
00:08823807,48
00:09823830,52
00:10823806,69
00:01924538,05
00:02926018,59
00:03924001,34
00:04923740,09
00:05923735,53
00:06923679,07
00:07923726,67
00:08923695,23
00:09923701,06
00:10923543,47
2 REPLIES 2
StatDave
SAS Super FREQ

Since it's always good to start by plotting the data, the picture from the following makes them look pretty different  I assume that the CPM values are decimal values even though you say they are counts, suggesting integers. 

proc sgplot;
series y=cpm x=time/group=id;
run;

There appears to be several IDs that spike upwards at time=2 and several that don't. If you want to compare the centers of the values, it's probably best to use a nonparametric test that doesn't assume a distribution for the values, such as this.

proc npar1way wilcoxon;
class id;
var cpm;
run;

Results suggest the centers of the ID distributions are not the same which is also suggested in the plot even when ignoring time 2.

SteveDenham
Jade | Level 19

Do you really want to test that the accelerometers are the "same", or is it sufficient to test for "at least one accelerometer is clearly different from the others". The reason is that the first is an equivalence test where two one-sided tests (TOST), or some generalization of that, would be used (ANOM ?), and you have to state what the boundaries on the equivalence interval are. For the second, an ANOVA would suffice (although you may want to fit a generalized linear model, depending on the residuals and a histogram plot). The issue with the equivalence test approach is that there are 36 pairwise comparisons, all of which would have to pass the equivalence criteria in order to say that the instruments were the "same". If one of the accelerometers is accepted as a standard, the number of comparisons drops to 8, and that may be tractable.

 

SteveDenham

SAS Innovate 2025: Call for Content

Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!

Submit your idea!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 2 replies
  • 466 views
  • 0 likes
  • 3 in conversation