pwrLaw {astrochron} | R Documentation |

## Generate power law (1/f) noise surrogates

### Description

Generate power law (1/f) noise surrogates, following the algorithm of Timmer and Konig (1995).

### Usage

```
pwrLaw(npts=1024,dt=1,mean=0,sdev=1,beta=2,fcut=0,nsim=1,genplot=T,verbose=T)
```

### Arguments

`npts` |
number of data points for 1/f surrogate time series |

`dt` |
sampling interval |

`mean` |
mean value for 1/f surrogate series |

`sdev` |
standard deviation for 1/f surrogate series |

`beta` |
power law coefficient. Positive number will yield a negative slope. |

`fcut` |
frequency cutoff: below this frequency a plateau will be modeled. Set to zero (default) for no plateau. |

`nsim` |
Number of surrogate series to generate |

`genplot` |
generate summary plots (T or F) |

`verbose` |
verbose output (T or F) |

### Details

These simulations use the random number generator of Matsumoto and Nishimura (1998). Power law noise series are generated following the algorithm of Timmer and Konig (1995).

### References

M. Matsumoto, and T. Nishimura, (1998),
*Mersenne Twister: A 623-dimensionally equidistributed uniform pseudo-random number generator*, ACM Transactions on Modeling and Computer Simulation, 8, 3-30.

J. Timmer and K. Konig (1995),
*On Generating Power Law Noise*, Astronomy and Astrophysics: v. 300, p. 707-710.

*astrochron*version 1.3 Index]