Well, I would point out to you that you only have 9 values of h which certainly is a problem in calculating 10 colmans(?) -- more generally, I'd suggest you learn more about the difference between vectorized operations and operations on vectors: things like `%*%` and outer() will help you greatly. If you are perhaps coming from a Matlab background, keep in mind that R's default behavior is to include the "." before all operations. E.g.,
R> x = 1:3 R> x^2 [1] 1 4 9 R> x * x [1] 1 4 9 R> x %*% x [,1] [1,] 14 R> outer(x,x) [,1] [,2] [,3] [1,] 1 2 3 [2,] 2 4 6 [3,] 3 6 9 Michael On Fri, Sep 9, 2011 at 2:32 PM, Halabi, Anan <anan.hal...@hp.com> wrote: > What I mean by the metric not created,**** > > Desired metric is of 1100 values by 110 rows (t), and 10 Colman (h)**** > > Thank, **** > > ** ** > > *From:* R. Michael Weylandt [mailto:michael.weyla...@gmail.com] > *Sent:* Friday, September 09, 2011 10:25 PM > *To:* Halabi, Anan > *Cc:* r-help@R-project.org > *Subject:* Re: [R] Reliability metric**** > > ** ** > > Your code works fine -- the warning message is just a heads up that the > object z/eta (which has length 111) does not quite match up with 1-h which > has length 10: usually this is a sign that the user defined something > improperly so R gives you a nice warning, but it doesn't halt the > calculation. > > As to what you mean about your metric not being created -- if it's "R" (odd > choice for a variable name given your platform) it sure seems to > exist....especially since you just gave us 111 values for it. > > Though I think you don't know what you mean by "[the] function R". Perhaps > you should read about how user defined functions are set up in R (the > fantastically useful insightful warning message giving language): > http://www.statmethods.net/management/userfunctions.html > > Michael Weylandt**** > > On Fri, Sep 9, 2011 at 2:14 PM, Halabi, Anan <anan.hal...@hp.com> wrote:** > ** > > > 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.**** > > ** ** > [[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.