Sample size calculations

Code:

m = c(1.5, 2, 3, 4, 5, 6, 7, 8, 9, 10)
r = 0:9/10
SDbase = 43.8
effect = 45.8 - 28

alpha = 0.05
power = 0.8

n = function(r, m) {
    SDdiff = sqrt(2) * SDbase * sqrt(1 - r)
    round(2 * (SDdiff/effect)^2 * (qnorm(1 - alpha/2) + qnorm(power))^2 * (1 + 
        1/m)^2)
}
result = outer(r, m, n)

Results (m in columns:, r in rows)

## 
## ----- ----- ----- ----- ----- ----- ----- ----- ----- ----- -----
##       1.5   2     3     4     5     6     7     8     9     10   
## 0     528   428   338   297   274   259   248   241   235   230  
## 0.1   475   385   304   267   246   233   223   217   211   207  
## 0.2   422   342   270   238   219   207   199   192   188   184  
## 0.3   370   299   237   208   192   181   174   168   164   161  
## 0.4   317   257   203   178   164   155   149   144   141   138  
## 0.5   264   214   169   149   137   129   124   120   117   115  
## 0.6   211   171   135   119   109   103   99    96    94    92   
## 0.7   158   128   101   89    82    78    74    72    70    69   
## 0.8   106   86    68    59    55    52    50    48    47    46   
## 0.9   53    43    34    30    27    26    25    24    23    23   
## ----- ----- ----- ----- ----- ----- ----- ----- ----- ----- -----