Thanks for your insights into physics topics. I do not doubt, I am maybe sometimes like a little child in these physics things, as I am a statistician (not a physcists).

However, on the other side it may also be interesting for you physicists how a statistician is planning to deal with the piled up amplitude problem by software programming in R. As most of you do not understand R code, I will write pseudocode, instead.

Code: Select all

```
# Assuming the voltage data from a data acqusistion card would a csv data file, which has a voltage figure on every line.
import.csv("mdata_00.csv" as mdata_00)
columns(mdata_00)= c("ln", "volt")
# (Non-noise) signal is greater than 0.1 V
mdata_01= data.frame(mdata_00, signal= (mdata_00$volt>= 0.1))
mdata_01= data.frame(mdata_01, voltmax= roll.max(madata_01$volt, 30))
mdata_01= data.frame(mdata_01, ampbegin= rep(NA, nrow(mdata_01), ampend= rep(NA, nrow(mdata_01))
mdata_01$ampbegin= (madata_01$signal== TRUE & mdata_01$signal.predecessor== FALSE)
mdata_01$ampend= (madata_01$signal== TRUE & mdata_01$signal.succesor== FALSE)
mdata02= select(mdata_01 where mdata_0$sinal==TRUE)
mdata02$ampno=rep(NA, nrow(mdata_02))
j= 1
for i in (nrow(mdata02) {
mdata02$ampno(i)= j
If mdata_01$ampbegin(i)== T {
j += 1
}
}
# Transposing the amplitudes into rows
ampdata_01= transpose(mdata02, by= ampno)
ampdata_01$rowlength = ...
# Assuming a gaussian bell curve for singular amplitudes and testing it for gaussian distribution
ntest_results= apply(test.gaussian(ampdata_01, distrution="normal", alpha= 0.05), by=row)
ntest_results$ampno= ...
# "SQL"-Joining the results back
ampdata_02= merge(ampdata_01, ntest_results, by.x= "ampno", by.y= "ampno")
# Calculating the numerical volt integral of the amplitudes
...
# Creating classes of volt integral
ampdata_02$voltintclass= class(ampdata_02$voltint, step= 0.1)
# Aggregating the sum of volt integral by ntestresult
intvoltdata_01= aggregate(ampdata_02$voltint, by= c("ntestresult", voltintclass")
# Assuming that piled amplitudes contain singular amplitudes in the same proportion as singular amplitudes occur
# Multiplying the singular amplitudes counts so that the integral voltage of multiple amplitudes gets compensated
...
```

As you (hopefully) can imagine now, counting piled up amplitudes correctly by using the statistics software R neither requires a genius nor requires self developing a new mathematical method. Instead simply one of the thousands of R packages is used to apply a widely universal statistical method on the data.

And I thought, as most small particles in physics are best described by a distribution, physics would also always do the same numerical programming (or similar in Matlab, SAS, Pyhton, ...).

But, concluding from your answers, this does not seem to be the case.

[/Your insight to statistical working]