cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
Browse apps to extend the software in the new JMP Marketplace
Choose Language Hide Translation Bar
manphoolfageria
Level III

Paired t test

Paired t test.pngHi I want to compare between two treatments: 1) No fertilizer and 2) Fertilizer. I wanted to know whether there is an increase in yield due to the fertilizer application. Two parallel blocks were selected, in the first block no fertilizer was applied whereas in the second block the fertilizer was applied. Fifteen samples were harvested from each block. I would like to use the paired t test and there are different methods in JMP. I thought to use the distribution and then test means. My question is that is adding the column “n” for “frequency” is correct? Please explain why or why not. If I use n as frequency in the data analysis the results changes as the degree of freedom increase from 14 to 29. Paired t test.png

 

12 REPLIES 12
cwillden
Super User (Alumni)

Re: Paired t test

No, do not add a frequency column.  Each pair is 1 data point in a paired t-test.  A paired t-test is mathematically equivalent to a 1-sample t-test on the paired differences.

-- Cameron Willden
manphoolfageria
Level III

Re: Paired t test

Thanks so much for your quick reply. I understand that the paired t test can be performed by using the following four methods: I am wondering whether all four methods are equally good. If not which one is the best and why? I got identical results.

 

1. Analyze > distribution- mean test

Picture1.png

2. Analyze-specialized modeling-matched pairs

Picture2.png

3. Fit Y by X > compare means> Each pair, students t test- detailed comparison report 

Picture3.png

4. Fit Model platform to do a two-way ANOVA

Picture4.png

cwillden
Super User (Alumni)

Re: Paired t test

In #4 you should have Rep as a random effect in order to get an equivalent result to the other methods.  There is actually another way, which is to plot your Diff column in the Distribution platform and do a 1-sample t-test.   Since they should all give the same answer, you could say they are equally good.  The Matched Pairs platform is specialized for this problem, so you'll get the most specialized/relevant output in your report there.  However, that platform requires the data to be in a different shape than the usual 'long' format.  If I'm want to avoid splitting my table to get it in the right shape, I might just use Fit Y by X with the pair ID column as the block, or Fit Model with a random effect on the pair ID column.

-- Cameron Willden
manphoolfageria
Level III

Re: Paired t test

Thanks so much again,

Sorry, could you please expain the following:

1. How to use rep as a random effect in #4? Could you please explain it by including a table.   

2. Could you include another table explaining Fit Y by X with the pair ID column as the block.

cwillden
Super User (Alumni)

Re: Paired t test

Here's an example using Gosset's Corn data from the JMP's Sample Data Library.  I added a pair ID column and then stacked the 2 yield columns (Kiln and Regular) to get to a long data format.  Make sure the pair ID column has a nominal modeling type.

 

From Fit Model:

1. Add the Yield column as the Y, Response variable.

2. Add the Treatment and Pair columns into the Construct Model Effects box.

3. Select (highlight) 'Pair' in the Construct Model Effects box and make it random by clicking the Attributes red-arrow > Random.

 

MatchedPairsInFitModel.PNG

 

 

You'll see the p-value for Yield has the same p-value as the Matched Pairs platform for the same data in the original wider format.MatchedPairsPvals.PNG

 

 

 

-- Cameron Willden
cwillden
Super User (Alumni)

Re: Paired t test

Here's how to get the same p-value in Fit Y by X with the same stacked data table for Gosset's Corn:FitYbyXMatchedPairs.PNG

 

 

FitYbyXResults.PNG

-- Cameron Willden

Re: Paired t test

Please clarify the basis for pairing the samples. How were two samples paired together?

manphoolfageria
Level III

Re: Paired t test

Field MapField Map

Hi Mark,

We wanted to compare between two seed types in potatoes (planting Cut seed vs planting whole seed). Cut seed is used as control which is a common practice. This trial was conducted at the grower’s field. The field was divided into two strips blocks, the first block (in yellow) was planted with control (cut seed) and the second block (light blue) was planted with the whole seed. These strip blocks are really long 400-foot. In total, twenty four samples (12 from each strip block) were harvested.  Each sample was harvested from a 10-foot long row plot and converted the yield in CWT per acre. The field variation was huge. The trial was conducted at grower’s field. So, we could not randomize the treatments. I have analyzed the data using the paired t with distribution and test mean, see below:

 

My question:

I understand that we should use the paired t test.

One of the great statistician suggested that we should create an extra column frequency and use it for the analysis. I have analyzed the data with and without using the frequency.

I wanted to know which method is incorrect and why?

Which method is correct?

I would appreciate your time.

Best regards

 

 

  • I have created difference and frequency columns and both

Jeff_Perkinson_1-1581778225562.png

 

 

See the results below: They are significant P >0.0232

Jeff_Perkinson_2-1581778225562.png

  • I have created difference and frequency column but this time I did not use the frequency column

Jeff_Perkinson_3-1581778225562.png

Results see below: Non significant P >0.1208

Jeff_Perkinson_4-1581778225563.png

Re: Paired t test

You description of blocks is yellow and blue, so the blocking (pairs) does not make sense,

 

If the cut and whole seed plantings occurred next to each other, but the 12 strips are separated over a variable fied, the blocking (pairs) make sense. Is the field homogeneous within a pair? That assumption is important.

 

Yu do not need any more than the two response columns, CONTROL and WHOLE SEED. Use the Matched Pairs platform in the Analyze menu. Enter both of the data columns in the Y role in the order above. Do not enter any other data columns in this analysis.

 

It seems that you expect a higher yield for the whole seed treatment (alternative hypothesis), so an upper-tailed paired t-test is called for.