Hi, Below is the code I write, I am trying to create a metric of h and t while the values is out of function R. First I have message warning Second, the metric is not created
> h <- seq(0.1, 0.9, by=0.1) > t <- seq(0,11000, by=100) > z <- cbind(t) > eta=10000 > beta=2 > R <- array (1:1100, dim= c(110,10)) > R= exp(-(z/eta*(1-h))^(beta*(1-h))) Warning messages: 1: In z/eta * (1 - h) : longer object length is not a multiple of shorter object length 2: In (z/eta * (1 - h))^(beta * (1 - h)) : longer object length is not a multiple of shorter object length > R t [1,] 1.0000000 [2,] 0.9995586 [3,] 0.9974647 [4,] 0.9919725 [5,] 0.9801987 [6,] 0.9572081 [7,] 0.9141353 [8,] 0.8341665 [9,] 0.6833618 [10,] 0.9892125 [11,] 0.9825766 [12,] 0.9727664 [13,] 0.9583520 [14,] 0.9370675 [15,] 0.9051390 [16,] 0.8559267 [17,] 0.7769472 [18,] 0.6423128 [19,] 0.9629361 [20,] 0.9521001 [21,] 0.9382259 [22,] 0.9201104 [23,] 0.8958341 [24,] 0.8622051 [25,] 0.8136275 [26,] 0.7395534 [27,] 0.6175830 [28,] 0.9246320 [29,] 0.9127580 [30,] 0.8982795 [31,] 0.8800821 [32,] 0.8564152 [33,] 0.8244048 [34,] 0.7790512 [35,] 0.7109198 [36,] 0.5996181 [37,] 0.8767646 [38,] 0.8671159 [39,] 0.8550329 [40,] 0.8394499 [41,] 0.8187308 [42,] 0.7902253 [43,] 0.7493487 [44,] 0.6874269 [45,] 0.5854263 [46,] 0.8215803 [47,] 0.8170930 [48,] 0.8098560 [49,] 0.7988930 [50,] 0.7827045 [51,] 0.7588540 [52,] 0.7231046 [53,] 0.6673758 [54,] 0.5736597 [55,] 0.7611968 [56,] 0.7642511 [57,] 0.7637415 [58,] 0.7588504 [59,] 0.7482636 [60,] 0.7297788 [61,] 0.6994875 [62,] 0.6498159 [63,] 0.5635906 [64,] 0.6975913 [65,] 0.7098972 [66,] 0.7174446 [67,] 0.7196221 [68,] 0.7153381 [69,] 0.7026455 [70,] 0.6779572 [71,] 0.6341555 [72,] 0.5547793 [73,] 0.6325678 [74,] 0.6551296 [75,] 0.6715521 [76,] 0.6814181 [77,] 0.6838614 [78,] 0.6771931 [79,] 0.6581383 [80,] 0.6199984 [81,] 0.5469392 [82,] 0.5677232 [83,] 0.6008650 [84,] 0.6265229 [85,] 0.6443861 [86,] 0.6537698 [87,] 0.6532207 [88,] 0.6397566 [89,] 0.6070646 [90,] 0.5398727 [91,] 0.5044222 [92,] 0.5478565 [93,] 0.5827146 [94,] 0.6086282 [95,] 0.6250023 [96,] 0.6305684 [97,] 0.6226043 [98,] 0.5951486 [99,] 0.5334377 [100,] 0.4437834 [101,] 0.4967086 [102,] 0.5404019 [103,] 0.5742126 [104,] 0.5975006 [105,] 0.6091056 [106,] 0.6065193 [107,] 0.5840942 [108,] 0.5275280 [109,] 0.3866770 [110,] 0.4478912 [111,] 0.4997913 > length(R) [1] 111 Anan Halabi Reliability Eng, R&D HP Scitex Tel: 972-9-8924648 mobil: 972-52-6624231 [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.