RenewalProcess[rdist]
represents a renewal process with interarrival times distributed according to rdist.
RenewalProcess
RenewalProcess[rdist]
represents a renewal process with interarrival times distributed according to rdist.
Details
- RenewalProcess is a continuous- or discrete-time and discrete-state process.
- RenewalProcess is a discrete-state and continuous-time or discrete-time process depending on rdist.
- The state
is the number of events in the interval 0 to
and
. - The distribution rdist can be any continuous or discrete distribution with positive domain.
- RenewalProcess can be used with such functions as Mean, PDF, Probability, and RandomFunction.
Examples
open all close allBasic Examples (2)
Scope (6)
Continuous Interarrival Distributions (2)
Renewal process with interarrival times that follow an Erlang distribution:
proc[k_, λ_] := RenewalProcess[ErlangDistribution[k, λ]];DiscretePlot[#[proc[3, 4][12], x], {x, 0, 30}, PlotRange -> All, ExtentSize -> 1 / 2, PlotLabel -> #]& /@ {PDF, CDF, HazardFunction, SurvivalFunction}Compute the probability of an event:
Probability[x[27] < 5, xproc[2, 1 / 3]]ListPlot[Table[RandomFunction[proc[2, λ], {200}], {λ, #}], PlotLegends -> (StringJoin["λ = ", ToString[#]]& /@ #)]&@{.2, .5, .7}Renewal process with interarrival times that follow a gamma distribution:
proc[α_, β_] = RenewalProcess[GammaDistribution[α, β]];DiscretePlot[#[proc[.7, 2][12], x], {x, 0, 30}, PlotRange -> All, ExtentSize -> 1 / 2, PlotLabel -> #]& /@ {PDF, CDF, HazardFunction, SurvivalFunction}Compute the probability of an event:
Probability[12 < x[12] ^ 3 + 12x[12] < 120, xproc[.7, 2]]NProbability[12 < x[12] ^ 3 + 12x[12] < 120, xproc[.7, 2]]ListPlot[Table[RandomFunction[proc[.7, β], {200}], {β, #}], PlotLegends -> (StringJoin["β = ", ToString[#]]& /@ #)]&@{.5, 1, 2}Discrete Interarrival Distributions (2)
Renewal process with interarrival times that follow a Pascal distribution:
proc[p_] := RenewalProcess[PascalDistribution[3, p]]DiscretePlot[#[proc[.4][70], x], {x, 0, 15}, PlotRange -> All, ExtentSize -> 1 / 2, PlotLabel -> #]& /@ {PDF, CDF, HazardFunction, SurvivalFunction}Compute the probability of an event:
Probability[x[20] < 5, xproc[.7]]ListPlot[Table[RandomFunction[proc[p], {300}], {p, #}], PlotLegends -> (StringJoin["p = ", ToString[#]]& /@ #)]&@{.2, .5, .7}Renewal process with interarrival times that follow a Borel–Tanner distribution:
proc[p_] := RenewalProcess[BorelTannerDistribution[p, 2]]DiscretePlot[#[proc[.4][70], x], {x, 0, 30}, PlotRange -> All, ExtentSize -> 1 / 2, PlotLabel -> #]& /@ {PDF, CDF, HazardFunction, SurvivalFunction}Compute the probability of an event:
Probability[x[100] < 5, xproc[.7]]ListPlot[Table[RandomFunction[proc[p], {1, 300}], {p, #}], PlotLegends -> (StringJoin["p = ", ToString[#]]& /@ #)]&@{.2, .5, .7}Parameter Estimation (2)
sample = RandomFunction[RenewalProcess[BetaDistribution[1, 3]], {10 ^ 3}];Estimate the distribution parameters from sample data:
EstimatedProcess[sample, RenewalProcess[BetaDistribution[a, b]]]Approximate value of the mean:
ListLinePlot@Table[Mean[RenewalProcess[WeibullDistribution[1, .7]][t]], {t, 0.5, 8, .5}]ListLinePlot@Table[Mean[RenewalProcess[GammaDistribution[1, 4]][t]], {t, 0.5, 8, .5}]Applications (1)
Messages arrive at a communication line according to a two-phase hyperexponential distribution with phase probabilities 0.4 and 0.6. The average arrival times for the two phases are 4.8 milliseconds and 0.8 milliseconds, respectively. Simulate the process for 100 milliseconds:
communicationLine = RenewalProcess[HyperexponentialDistribution[{0.4, 0.6}, {1 / 4.8, 1 / 0.8}]];ListStepPlot[RandomFunction[communicationLine, {0, 100}], Filling -> Axis]Mean number of arrivals during the first 10 milliseconds:
Mean[communicationLine[10]]Compare with the value obtained from simulation of the time slice distribution:
Mean[RandomVariate[communicationLine[10], 10 ^ 4]]//NProperties & Relations (2)
RenewalProcess is a jump process:
path = RandomFunction[RenewalProcess[WeibullDistribution[2, .7]], {10}];
f = path["PathFunction"];
jumps = path["Times"];Plot[f[t], {t, 0, 10}, Exclusions -> jumps, ExclusionsStyle -> Red]RenewalProcess is not weakly stationary for any distribution:
WeakStationarity[RenewalProcess[GammaDistribution[a, b]]]WeakStationarity[RenewalProcess[PascalDistribution[n, p]]]Possible Issues (1)
Closed forms are not available for most properties:
𝒫 = RenewalProcess[BetaDistribution[1, 3]];Mean[𝒫[1]]Obtain approximate values using inexact input or simulation:
Mean[𝒫[1.]]Mean[RandomVariate[𝒫[1], 10 ^ 4]]//NRelated Guides
History
Text
Wolfram Research (2012), RenewalProcess, Wolfram Language function, https://reference.wolfram.com/language/ref/RenewalProcess.html.
CMS
Wolfram Language. 2012. "RenewalProcess." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/RenewalProcess.html.
APA
Wolfram Language. (2012). RenewalProcess. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/RenewalProcess.html
BibTeX
@misc{reference.wolfram_2026_renewalprocess, author="Wolfram Research", title="{RenewalProcess}", year="2012", howpublished="\url{https://reference.wolfram.com/language/ref/RenewalProcess.html}", note=[Accessed: 12-June-2026]}
BibLaTeX
@online{reference.wolfram_2026_renewalprocess, organization={Wolfram Research}, title={RenewalProcess}, year={2012}, url={https://reference.wolfram.com/language/ref/RenewalProcess.html}, note=[Accessed: 12-June-2026]}