simulate_etube: Simulate Random Elliptical Tubes (ETReps)

In ETRep: Analysis of Elliptical Tubes Under the Relative Curvature Condition

Simulate Random Elliptical Tubes (ETReps)

Description

Generates random samples of ETReps based on a reference tube with added variation.

Usage

simulate_etube(
  referenceTube,
  numberOfSimulation,
  sd_v = 10^-10,
  sd_psi = 10^-10,
  sd_x = 10^-10,
  sd_a = 10^-10,
  sd_b = 10^-10,
  rangeSdScale = c(1, 2),
  plotting = TRUE
)

Arguments

referenceTube	List containing ETRep information as the reference.
numberOfSimulation	Integer, number of random samples.
sd_v	Standard deviations for various parameters.
sd_psi	Standard deviations for various parameters.
sd_x	Standard deviations for various parameters.
sd_a	Standard deviations for various parameters.
sd_b	Standard deviations for various parameters.
rangeSdScale	Numeric range for random scaling.
plotting	Logical, enables visualization of samples (default is FALSE).

Value

List of random ETReps.

References

Taheri, M., Pizer, S. M., & Schulz, J. (2024). "The Mean Shape under the Relative Curvature Condition." arXiv. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.48550/arXiv.2404.01043")}

Taheri Shalmani, M. (2024). "Shape Statistics via Skeletal Structures." University of Stavanger. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.13140/RG.2.2.34500.23685")}

Examples

# Load tube
data("colon3D")
#Set Parameters
sd_v<-sd_psi<-1e-03
sd_x<-sd_a<-sd_b<-1e-04
numberOfSimulation<-4
random_Tubes<-
  simulate_etube(referenceTube = colon3D,
                 numberOfSimulation = numberOfSimulation,
                 sd_v = sd_v,
                 sd_psi = sd_psi,
                 sd_x = sd_x,
                 sd_a = sd_a,
                 sd_b = sd_b,
                 rangeSdScale = c(1, 2),
                 plotting = FALSE)
# Plotting
## Not run: 
plot_Elliptical_Tube(random_Tubes[[1]], add = FALSE)
plot_Elliptical_Tube(random_Tubes[[2]], add = TRUE)
plot_Elliptical_Tube(random_Tubes[[3]], add = TRUE)
plot_Elliptical_Tube(random_Tubes[[4]], add = TRUE)
 
## End(Not run)

ETRep documentation built on Nov. 5, 2025, 6:34 p.m.