I. EXAMPLES OF USE

1. 	4+4						# result: 8
2. 	(4+4)/2^3				# result: 1
3. 	abs((4+4)/2^3 - 30)		# result: 29
4. 	mean (c(1,2,3))			# result: 2
5. 	sd (c(1,2,3))			# result: 1
6. 	sample (c(1,2,3),1)		# result: any one random number from the set, i.e. 1 or 2 or 3
7. 	qchisq(.95, df=1)		# result: 3.841459

II. REAL DATA: Use of definite article in BNC64 – relative frequencies per 10k

// calculate mean of all 64 values

mean(c(294.62, 238.77, 285.4, 366.05, 257.86, 347.29, 354.55, 248.91, 325.27, 390.48, 304.5, 302.17, 312.99, 424.1, 324.47, 295.63, 300.41, 311.62, 342.86, 343.47, 301.82, 318.22, 377.53, 331.43, 158.2, 360.65, 322.83, 407.28, 142.19, 333.42, 338.24, 380.7, 277.01, 189.27, 239.81, 276.36, 189.73, 305.19, 233.33, 290.14, 208.62, 246.82, 298.14, 291.6, 270.09, 81.87, 307.3, 198.85, 297.49, 276.29, 379.49, 264.2, 273.59, 244.29, 384.37, 275.91, 403.3, 355.04, 306.37, 150.32, 343.2, 275.76, 323.22, 322.95))

// calculate SD (standard deviation)  of all 64 values

sd(c(294.62, 238.77, 285.4, 366.05, 257.86, 347.29, 354.55, 248.91, 325.27, 390.48, 304.5, 302.17, 312.99, 424.1, 324.47, 295.63, 300.41, 311.62, 342.86, 343.47, 301.82, 318.22, 377.53, 331.43, 158.2, 360.65, 322.83, 407.28, 142.19, 333.42, 338.24, 380.7, 277.01, 189.27, 239.81, 276.36, 189.73, 305.19, 233.33, 290.14, 208.62, 246.82, 298.14, 291.6, 270.09, 81.87, 307.3, 198.85, 297.49, 276.29, 379.49, 264.2, 273.59, 244.29, 384.37, 275.91, 403.3, 355.04, 306.37, 150.32, 343.2, 275.76, 323.22, 322.95))

// calculate mean of random 10 values

mean(sample(c(294.62, 238.77, 285.4, 366.05, 257.86, 347.29, 354.55, 248.91, 325.27, 390.48, 304.5, 302.17, 312.99, 424.1, 324.47, 295.63, 300.41, 311.62, 342.86, 343.47, 301.82, 318.22, 377.53, 331.43, 158.2, 360.65, 322.83, 407.28, 142.19, 333.42, 338.24, 380.7, 277.01, 189.27, 239.81, 276.36, 189.73, 305.19, 233.33, 290.14, 208.62, 246.82, 298.14, 291.6, 270.09, 81.87, 307.3, 198.85, 297.49, 276.29, 379.49, 264.2, 273.59, 244.29, 384.37, 275.91, 403.3, 355.04, 306.37, 150.32, 343.2, 275.76, 323.22, 322.95), 10))

// calculate SD of random 10 values
sd(sample(c(294.62, 238.77, 285.4, 366.05, 257.86, 347.29, 354.55, 248.91, 325.27, 390.48, 304.5, 302.17, 312.99, 424.1, 324.47, 295.63, 300.41, 311.62, 342.86, 343.47, 301.82, 318.22, 377.53, 331.43, 158.2, 360.65, 322.83, 407.28, 142.19, 333.42, 338.24, 380.7, 277.01, 189.27, 239.81, 276.36, 189.73, 305.19, 233.33, 290.14, 208.62, 246.82, 298.14, 291.6, 270.09, 81.87, 307.3, 198.85, 297.49, 276.29, 379.49, 264.2, 273.59, 244.29, 384.37, 275.91, 403.3, 355.04, 306.37, 150.32, 343.2, 275.76, 323.22, 322.95), 10))