search for: sd3

Dear all: Assume I have 3 distributions, x1, x2, and x3. x1 ~ normal(mu1, sd1) x2 ~ normal(mu2, sd2) x3 ~ normal(mu3, sd3) y1 = x1 + x2 y2 = x1 + x3 Now that the data I can observed is only y1 and y2. It is easy to estimate (mu1+m2), (mu1+mu3), (sd1^2+sd2^2) and (sd1^2+sd3^2) by EM algorithm since y1 ~ normal(mu1+mu2, sqrt(sd1^2+sd2^2)) and y2 ~ normal(mu1+mu3, sqrt(sd1^2+sd3^2)) However, I want to estimate mu1, mu2...

...       DEF = c(22,28,20,20,28,26,29,18,24,21),                                           GHI = c(32,27,32,36,37,37,34,23,25,32),                                           JKL = c(47,60,60,43,62,38,44,53,61,41)) sd_prices <- c(3.3483,3.9497,4.9721,9.3387)   # standard deviations say(sd1, sd2, sd3, sd4)                         weights <- c(0.10, 0.25, 0.20, 0.45)      # say (w1, w2, w3, w4) I need to calculate the standard deviation of the portfolio. The formula is stdev_portfolio = sqrt((w1*sd1)^2+(w2*sd2)^2+(w3*sd3)^2+(w4*sd4)^2 +                         2*w1*w2*sd1*sd2*correlation(A...

Hi, I''m running Opensolaris 2009.06, and I''m facing a serious performance loss with ZFS ! It''s a raidz1 pool, made of 4 x 1TB SATA disks : zfs_raid ONLINE 0 0 0 raidz1 ONLINE 0 0 0 c7t2d0 ONLINE 0 0 0 c7t3d0 ONLINE 0 0 0 c7t4d0 ONLINE 0 0

...get the power of the test for modified bartlett's test..but I'm not really sure that my coding is right.. #normal distribution unequal variance asim<-5000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) n1<-20 n2<-20 n3<-20 mu<-0 sd1<-sqrt(25) sd2<-sqrt(50) sd3<-sqrt(100) g1<-rnorm(n1,mu,sd1) g2<-rnorm(n2,mu,sd2) g3<-rnorm(n3,mu,sd3) x=c(g1,g2,g3) group=c(rep(1,n1),rep(2,n2),rep(3,n3)) N=60 k=3 v1=var(g1) v2=var(g2) v3=var(g3) #pooled variance A=((n1-1)*v1+(n2-1)*v2+(n3-1)*v3)/(N-k) #calculate B B=((N-k)*(log(A)))-((n1-1)*log(v1)+(n2-1)*log(v2...

...ior to be the new prior and now I cant figure out how to do the code below in a loop. The data below is dummy data, I used a starting mu of 1, and starting precision of 0. bayes.analysis.treat<-function(mu0,p0){ n1 = 5 n2 = 10 n3 = 15 ybar1 = 12 ybar2 = 13 ybar3 = 14 sd1 = 2 sd2 = 3 sd3 = 4 #posterior 1 var1 = sd1^2 #sample variance p1 = n1/var1 #sample precision p1n = p0+p1 mu1 = ((p0)/(p1n)*mu0)+((p1)/(p1n))*ybar1 sigma1 = 1/sqrt(p1n) #posterior 2 var2 = sd2^2 #sample variance p2 = n2/var2 #sample precision p2n = p1n+p2 mu2 = ((p1n)/(p2n)*mu1)+((p2)/(p2n))*ybar2 sigm...

Hi all, I have time series data in a matrix format 20 rows x 205 columns and have been trying to replace outliers with NA. My first column contains the outliers threshold (3 Standard deviations) for each row - here's a bit of the first row sd3 V1 V2 V3 V4 V5 V6 V7 V8 V9 1 13.03267 1797157 75 84 58 -1.958649 0.048775 2.056198 8.063622 3.071045 What I want is a statment that says if the absolute value of [,5:205] <= column 1 keep the value, otherwise replace with NA I've tried this...

hello, I would like to create the following barplot: I have 4 different data sets (same length + stddev for each data point) data1 sd1 data2 sd2 data3 sd3 data4 sd4 now, I'd like to plot in the following way: data1[1],data2[1],data3[1],data4[1] with it's sd-values side-by-side at one x-axis label (named "position 1") and each bar in different colors. data1[2],data2[2],data3[2],data4[2] at the next x-axis label (named "posi...

...) Residuals: Min 1Q Median 3Q Max -0.0049787 -0.0012948 0.0000703 0.0009653 0.0063192 Coefficients: Estimate Std. Error t value Pr(>|t|) sd1 -0.0002012 0.0007653 -0.263 0.793820 sd2 0.0013339 0.0007616 1.752 0.086378 . sd3 0.0007372 0.0007947 0.928 0.358348 Y.dl1 -0.2215246 0.1600532 -1.384 0.172875 C.dl1 0.9220846 0.1646314 5.601 1.08e-06 *** Y.dl2 -0.1600219 0.1273838 -1.256 0.215245 C.dl2 0.4712112 0.2124175 2.218 0.031401 * Y.l3 -0.3227708 0.0881079 -3.663 0.00...

...device r/s w/s kr/s kw/s wait actv svc_t %w %b fd0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd2 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd5 0.0 118.0 0.0 15099.9 0.0 35.0 296.7 0 100 sd6 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd7 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 0 sd...

Hi All: This was discussed in the R-developers list (see the thread starting at http://tolstoy.newcastle.edu.au/R/e3/devel/ 07/12/0560.html). It has to do with the behavior of sd() when the entire vector is NA. The behavior has changed between 2.6 and 2.7.1 as follows: Run in Version 2.7.1 > tt<-rep(NA, 10) > mean(tt, na.rm=T) [1] NaN >sd(tt, na.rm=T) Error in var(x,

I am needing to create a new identifier variable for a data set which has no ID variable in the original file. I am basically wanting to take the count of each row and add "sd" in front of it - so it would look like so sd1 sd2 sd3 sd4 sd5 etc....... I have no idea how to do this. I am a SAS user trying to learn R. This question may have been answered previously, but I could be searching for it incorrectly. Thank you for any assistance you can provide. HyunJee -- View this message in context: http://r.789695.n4.nabble.com...

...454.7 0.0 47168.0 0.0 0.0 5.7 12.6 0 74 sd4 440.7 0.0 45825.9 0.0 0.0 5.5 12.4 0 78 sd6 445.7 0.0 46239.2 0.0 0.0 6.6 14.7 0 79 sd7 452.7 0.0 46850.7 0.0 0.0 6.0 13.3 0 79 sd8 460.7 0.0 46947.7 0.0 0.0 5.5 11.8 0 73 sd3 426.7 0.0 43726.4 0.0 5.6 0.8 14.9 73 79 sd5 424.7 0.0 44456.4 0.0 6.6 0.9 17.7 83 90 sd9 430.7 0.0 44266.5 0.0 5.8 0.8 15.5 78 84 sd10 421.7 0.0 44451.4 0.0 6.3 0.9 17.1 80 87 sd11 421.7 0.0 44196.1 0.0 5.8 0.8 15.8 75 8...

...1Q Median 3Q Max -0.027598 -0.012836 -0.003395 0.015523 0.056034 Coefficients: Estimate Std. Error t value Pr(>|t|) ect1 -0.212955 0.064354 -3.309 0.00185 ** sd1 -0.057653 0.010269 -5.614 1.16e-06 *** sd2 -0.016305 0.009177 -1.777 0.08238 . sd3 -0.040859 0.008767 -4.660 2.82e-05 *** LRM.dl1 0.049816 0.191992 0.259 0.79646 LRY.dl1 0.075717 0.157902 0.480 0.63389 IBO.dl1 -1.148954 0.372745 -3.082 0.00350 ** IDE.dl1 0.227094 0.546271 0.416 0.67959 > sjd.vecm.rls$beta ect1 LRM.l2 1.000000 LRY....

...g*(exp((h*b^2)/2)) y<-(exp(g*c)-1)/g*(exp((h*c^2)/2)) z<-(exp(g*d)-1)/g*(exp((h*d^2)/2)) ######################## FIXED SYMMETRIC TRIMMED MEAN ############################### asim<-5000 pv<-rep(NA, asim) for(j in 1:asim) { print(j) set.seed(j) n1=15 n2=15 n3=15 n4=15 miu=0 sd1=1 sd2=1 sd3=1 sd4=1 a=rnorm(n1,miu,sd1) b=rnorm(n2,miu,sd2) c=rnorm(n3,miu,sd3) d=rnorm(n4,miu,sd4) ## data transformation g=0 h=0 w<-a*exp(h*a^2/2) x<-b*exp(h*b^2/2) y<-c*exp(h*c^2/2) z<-d*exp(h*d^2/2) mat1<-sort(w) mat2<-sort(x) mat3<-sort(y) mat4<-sort(z) alpha=0.15...

...sometimes it also has load/store for iteration variables in basic blocks containing "for.body". Therefore, this approach can not solve my problem. Any suggestion for my problem? Thanks, Henry ------------------------------------ [1]. Minjang Kim, Hyesoon Kim, and Chi-Keung Luk. 2010. SD3: A Scalable Approach to Dynamic Data-Dependence Profiling, MICRO2010 -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.llvm.org/pipermail/llvm-dev/attachments/20141211/0fa09980/attachment.html>

...SMART self-test Apr 27 09:30:03 weston sata: [ID 674233 kern.notice] SATA1 compatible Apr 27 09:30:03 weston sata: [ID 349649 kern.notice] queue depth 32 Apr 27 09:30:03 weston sata: [ID 349649 kern.notice] capacity = 145226112 sec tors Apr 27 09:30:03 weston scsi: [ID 193665 kern.notice] sd3 at ahci0: target 2 lun 0 Apr 27 09:30:03 weston genunix: [ID 936769 kern.notice] sd3 is /pci at 0,0/pci15d9, 8180 at 1f,2/disk at 2,0 Apr 27 09:30:04 weston genunix: [ID 408114 kern.notice] /pci at 0,0/pci15d9,8180 at 1f ,2/disk at 2,0 (sd3) online ........ Apr 27 09:30:04 weston genunix: [ID 408...

...> asim<-5000 > pv<-rep(NA, asim) > for(j in 1:asim) > { > print(j) > set.seed(j) > n1=15 > n2=15 > n3=15 > n4=15 > miu=0 > sd1=1 > sd2=1 > sd3=1 > sd4=1 > > a=rnorm(n1,miu,sd1) > b=rnorm(n2,miu,sd2) > c=rnorm(n3,miu,sd3) > d=rnorm(n4,miu,sd4) > > ## data transformation > g=0 > h=0 > > w<-a*exp(h*a^2/2) > x<-...

...uot;. Therefore, this > approach can not solve my problem. > > > > Any suggestion for my problem? > > > > Thanks, > > Henry > > ------------------------------------ > > > > > > > > [1]. Minjang Kim, Hyesoon Kim, and Chi-Keung Luk. 2010. SD3: A Scalable > Approach to Dynamic Data-Dependence Profiling, MICRO2010 > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.llvm.org/pipermail/llvm-dev/attachments/20141211/a15961c0/attachment.html>

Just look at this. I thought all the restarting resilver bugs were fixed, but it looks like something odd is still happening at the start: Status immediately after starting resilver: # zpool status pool: rc-pool state: DEGRADED status: One or more devices has experienced an unrecoverable error. An attempt was made to correct the error. Applications are unaffected. action: Determine

Hi all, just after sending a message to sunmanagers I realized that my question should rather have gone here. So sunmanagers please excus ethe double post: I have inherited a X4140 (8 SAS slots) and have just setup the system with Solaris 10 09. I first setup the system on a mirrored pool over the first two disks pool: rpool state: ONLINE scrub: none requested config: NAME