Running examples in 'statgenHTP-Ex.R' failed
The error most likely occurred in:

> ### Name: estimateSplineParameters
> ### Title: Extract estimates from fitted splines.
> ### Aliases: estimateSplineParameters
> 
> ### ** Examples
> 
> ### Estimate parameters for fitted P-splines.
> 
> ## Run the function to fit P-splines on a subset of genotypes.
> subGeno <- c("G160", "G151")
> fit.spline <- fitSpline(inDat = spatCorrectedVator,
+                         trait = "EffpsII_corr",
+                         genotypes = subGeno,
+                         knots = 50)
> 
> ## Estimate the maximum value of the predictions at the beginning of the time course.
> paramVator <- estimateSplineParameters(x = fit.spline,
+                                        estimate = "predictions",
+                                        what = "max",
+                                        timeMin = 1527784620,
+                                        timeMax = 1528500000,
+                                        genotypes = subGeno)
> head(paramVator)
  genotype plotId max_predictions max_timeNumber       max_timePoint
1     G151  c13r1       0.6916861         636000 2018-06-08 01:17:00
2     G151 c15r48       0.7030388         651200 2018-06-08 05:30:20
3     G151  c1r16       0.6953883         651200 2018-06-08 05:30:20
4     G151 c21r25       0.7324238         166400 2018-06-02 14:50:20
5     G151   c2r7       0.6942117         651200 2018-06-08 05:30:20
6     G151  c6r46       0.6897347         655200 2018-06-08 06:37:00
> 
> ## Create a boxplot of the estimates.
> plot(paramVator, plotType = "box")
> 
> ## Estimate the minimum and maximum value of the predictions.
> paramVator2 <- estimateSplineParameters(x = fit.spline,
+                                         estimate = "predictions",
+                                         what = c("min", "max"),
+                                         genotypes = subGeno)
> head(paramVator2)
  genotype plotId min_predictions min_timeNumber       min_timePoint
1     G151  c13r1       0.4256623         795200 2018-06-09 21:30:20
2     G151 c15r48       0.3575513         775200 2018-06-09 15:57:00
3     G151  c1r16       0.3769714         775200 2018-06-09 15:57:00
4     G151 c21r25       0.3013801         787200 2018-06-09 19:17:00
5     G151   c2r7       0.3558154         776800 2018-06-09 16:23:40
6     G151  c6r46       0.3613267         776800 2018-06-09 16:23:40
  max_predictions max_timeNumber       max_timePoint
1       0.7400419        1212800 2018-06-14 17:30:20
2       0.7367294        1207200 2018-06-14 15:57:00
3       0.7312961        1214400 2018-06-14 17:57:00
4       0.7337103        1214400 2018-06-14 17:57:00
5       0.7353195        1204000 2018-06-14 15:03:40
6       0.7334287        1214400 2018-06-14 17:57:00
> 
> 
> ### Estimate parameters for fitted HDM-splines.
> 
> ## The data from the Phenovator platform have been corrected for spatial
> ## trends and outliers for single observations have been removed.
> 
> ## We need to specify the genotype-by-treatment interaction.
> ## Treatment: water regime (WW, WD).
> spatCorrectedArch[["treat"]] <- substr(spatCorrectedArch[["geno.decomp"]],
+                                       start = 1, stop = 2)
> spatCorrectedArch[["genoTreat"]] <-
+   interaction(spatCorrectedArch[["genotype"]],
+              spatCorrectedArch[["treat"]], sep = "_")
> 
> ## Fit P-Splines Hierarchical Curve Data Model for selection of genotypes.
> fit.psHDM  <- fitSplineHDM(inDat = spatCorrectedArch,
+                           trait = "LeafArea_corr",
+                           genotypes = c("GenoA14_WD", "GenoA51_WD",
+                                        "GenoB11_WW", "GenoB02_WD",
+                                        "GenoB02_WW"),
+                           time = "timeNumber",
+                           pop = "geno.decomp",
+                           genotype = "genoTreat",
+                           plotId = "plotId",
+                           difVar = list(geno = FALSE, plot = FALSE),
+                           smoothPop = list(nseg = 4, bdeg = 3, pord = 2),
+                           smoothGeno = list(nseg = 4, bdeg = 3, pord = 2),
+                           smoothPlot = list(nseg = 4, bdeg = 3, pord = 2),
+                           weights = "wt",
+                           trace = FALSE)
Error: 'crossprod' is not an exported object from 'namespace:Matrix'
Execution halted