search for: deviation

I have a graph of residuals and I am attempting to get a list of the indexes of each time the residual is greater than 2 standard deviations or less than -2 standard deviations, but only the first point of the section. And then I'd also need the first point where the point returns to the range between +/- 2 standard deviations. So basically if my standard deviation=1 and my residuals=c(1, 2.1, 3, 4, 3, 1, 0, -4, -1) I want it to gr...

Hi, I have problem getting the standard deviation from the manova output. I have used the manova function: myfit <- manova(cbind(y1, y2) ~ x1 + x2 + x3, data=mydata) . I tried to get the predicted values and their standard deviation by using: predict(myfit, type="response", se.fit=TRUE) But the problem is that I do...

Hi everyone, I have just started using R, and I have a simple question. How can I get the Standard deviation for rows. basically I am looking for something like "rowMeans()" but for Standard deviation (I tried "rowSds()" didn't exist) Thanks, -- View this message in context: http://www.nabble.com/Standard-deviation-for-rows-tp19998106p19998106.html Sent from the R help mailing...

I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow. > integrate(dnorm, -1.96, 1.96) 0.9500042 with absolute error < 1e-11 > mean = 2.0 > sd = 1.0 > integrate(dnorm, -1.96, 1.96) 0.9500042 with absolute error < 1e-11 >...

Dear R users, This is more a statistical question rather than an R question. I'd appreciate it if you can give me some suggestions. I have a sample of a time series (sample size 500, fat tail in density). I am trying to calculate the Standard error of standard deviation of a sub-block-sample (sample size 250). I take 100 this kind of sub-block-sample, randomly. For these 100 subsamples, I use the following 3 methods to calculate the standard error of standard deviation: 1. From book "Handbook of applicable mathematics", Walter Ledermann (chief editor) V...

...tafor package. However I have some difficulties setting up my raw data to enter it to the meta-analysis models. I have a group of subjects that have been measured in two continuous variables (A & B). I have the information about the mean of the two variables, the group size (n) and the standard deviations of the two variables. Now I would like to average both variables (A & B) and get the mean and standard deviation of the merged variable (C). As for the mean this would be quiet easy: I would just take the mean of mean A and mean B to get the mean of C. However for the standard deviation this s...

...tAsymptote(sXY) NLSstClosestX(sXY, (a+d)/2) ------------------------ I think the problem arises when the target y value is exactly equal to one of the y values in sXY and can be fixed by trapping that situation thus: -------------------- NLSstClosestX.sortedXyData <- function (xy, yval) { deviations <- xy$y - yval if (any(deviations==0)) xy$x[match(0, deviations)] else { # new line inserted if (any(deviations <= 0)) { dev1 <- max(deviations[deviations <= 0]) lim1 <- xy$x[match(dev1, deviations)] if (all(deviations <= 0)) {...

I need some help for a curious question of a friend of mine. She usually does some experiments (3-5 repeats for each exp) and then she calculates mean and standard deviation. In microsoft excel she writes something like the following sample mean sd a 1.25 0.35 b 2.65 0.65 c 3.45 0.50 She can do a vertical barplot graph just giving mean value and specifying the standard deviation value for each bar in the graph. My friend wants to know if...

I'd like to use vectorization to take a 4 point moving window on standard deviation on the close column and create another variable (st.dev) in the dataframe. Here's the dataframe head(xyz) Date Close 1 2011-01-28 56.42 2 2011-01-27 57.37 3 2011-01-26 56.48 4 2011-01-25 56.39 5 2011-01-24 55.74 6 2011-01-21 55.46 So the first 3 elements to the new st.dev column woul...

In the manual page for prcomp(), it says that sdev is "the standard deviations of the principal components (i.e., the square roots of the eigenvalues of the covariance/correlation matrix, though the calculation is actually done with the singular values of the data matrix)." ?However, this is not what I'm finding. ?The values appear to be the standard deviations of a...

< Does anyone know why they're called random deviates, as opposed to random numbers?> Others will probably give you some technical reason about random numbers can be considered as random deviates from a mean (I think at least the 1875 Galton paper at http://www.mugu.com/galton/ uses similar terminology (I'm not claiming this is the earliest use - just the easiest to access at the

Hi, Could anyone give suggestions how to plot a scatter plot with 1 standard deviation for each point. To make it clearer, here is a simple example: the scatterplot is plot(X, Y), but I want to add 1 standard deviation according to the value of Z for each Y. X Y Z 1 3.5 1.1 . . . . . . . . . Thanks a lot in advance. FD [[alte...

Hi, I have two (different types of) measurements, say X and Y, resulting from the same set of experiments. So X and Y are paired: (x_1, y_1), (x_2, y_2), ... I am trying to calculate the standard deviation of Y with respect to X. In other words, in terms of the scatter plot of X and Y, I would like to divide it into bins along the X-axis and for each bin calculate the standard deviation along the Y results in that bin. (Though I am not totally sure, this seems to remind me of the conditional expectat...

Hi, I have the following data about courses (504) in a university, two attributes about the proportion of resources used (#resources_used / #resources_available), namely the average and the standard deviation. Thus I have: [1] n=504 rows [2] 1 id column and 2 attributes Here's a sample of the data: courseid,average,std 12741,1,0 17161,1,0 12514,1,0 12316,0.8666666692648178,0.26090261464799325 2467,0.8623188442510107,0.24920700355307424 3047,0.85,0.2314550249431379 1747,0.8481481481481481,0.2307844...

I don?t know if my question is answerable, but it is worth a try. I have a data set that I am trying to analyze in R for a course and the instructions were to get a standard deviation which I already computed in R and use that number and change it to a biased standard deviation?.(I have the two equations and I understand the difference between the two and how the unbiased has the degree of freedom?..I just do not know how use my standard deviation and transform it in R to a bias...

Hi R-helpers, I have been struggling with calculating row and column statistics, e.g. standard deviation. I know that > datac$Mean<-rowMeans(datac,na.rm=TRUE) will give me row means. I have tried to replicate those row means with the apply function: > datac$Mean2<-apply(datac,2,mean) so that I can replace the function argument with "sd" (instead of mean) to get standard deviat...

Is there a simple function in R to get a moving standard deviation (i.e. for the last x samples)? My goal is to plot bollinger bands around a moving average for price data. I use kernel smoothing for the moving average. cheers and thanks! over and out -- doktora

Hi, R seems to use the "1/n-1"-factor calculating the standard-deviation sd(). If i wat to get the "non-sample" standard-deviation i use sqrt(sd(x)^2*((n-1)/n)) Is there a parameter to get the sd()-function using the "1/n" factor directly? Or is there any other function to do so? Thank you in advance :-) Best, Eli

Full_Name: Maria Rizzo Version: 2.10.1 (2009-12-14) OS: Windows XP SP3 Submission from: (NULL) (72.241.75.222) platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor

The model matrix for the code at the end the email is shown below. Since the model matrix doesn't have -1, I think that it is made of dummy regressors rather than deviation regressors. I'm wondering how to make a model matrix using deviation regressors. Could somebody let me know? > model.matrix(aaov) (Intercept) A2 B2 B3 A2:B2 A2:B3 1 1 0 0 0 0 0 2 1 0 0 0 0 0 3 1 0 0 0 0 0 4 1 0...