search for: ith

...gt; 023 2 5 > 023 4 6 > 023 5 7 > 412 2 5 > 412 3 4 > 412 4 6 > 412 7 9 > 220 5 7 > 220 4 8 > 220 9 8 > ...... > Now I want to compute the following double summation: sum_{i=1}^k sum_{j=1}^{n_i}(x_{ij}-mean(x_i))*(y_{ij}-mean(y_i)) i is from 1 to k, indexing the ith subject id; and j is from 1 to n_i, indexing the jth observation for the ith subject. in the above expression, mean(x_i) is the mean of x values for the ith subject, mean(y_i) is the mean of y values for the ith subject. Is there a simple way to do this in R?

...ver time, using multiple sites: y[ij]=mu+b[j]+a[i]+w[j]*(beta +t[i])+c[ij] where: y[ij]= a response variable at site i and year j mu = fixed intercept Beta=fixed slope w[j]=constant representing the jth year (covariate) b[j]=random effect of jth year, iid N(0,sigma2[b]) a[i]=random effect of the ith site, iid N(0, sigma2[a]) t[i]=random effect of ith site, iid N(0, sigma2[t]) c[ij]=random error associated with ith site and jth year I would like to assume that an unstructured relationship applies to a[i] and t[i] (i.e., I would like to assume that the random effects a[i] and t[i] are drawn fro...

..."BM" [9] "BF" "IM" "IF" "AM" "AF" "Yr" how do i loop over a subset of variables in this frame to create a vector of length equal to the number of variables in the subset such that the vector's ith element is the result of an aggregate fncn applied to the ith variable? eg, sumvar<-c(sum(WM),...,sum(AF)) if i try for(i in 4:13){ sumVar<- sum(popA[,i]) sumVar2<-c(sumVar) } it returns > sumVar2 [1] 287567 > length(sumVar2) [1] 1 which is only the value at the last spot. i...

...adopt a value as a function of variables from the previous iteration. When I comment out the ?if else? statement, let W = Wo, and adjust the rest of the file so I am reading the first row of the variables specified instead of entire columns in my calculations, the calculations all work nicely without errors - I get the numbers I ought to in one row of ?bio?. When I activate the loop and try to run the file with ?the if .else?, it looks like it wants to work, but here?s where it stalls; Error in "[<-"(*tmp*, i, value = ((bio[i - 1, 1] * bio[i - 1, 9])/Ef)) : nothing to rep...

Dear all, I'm trying to get openmp to work with llvm, either with static compiler or with jit, however with no success. I've followed this: http://blog.llvm.org/2015/05/openmp-support_22.html [http://llvm.org/img/DragonSmall.png]<http://blog.llvm.org/2015/05/openmp-support_22.html> LLVM Project Blog: OpenMP Support OpenMP support i...

Hello everyone. I recently boght a Nokia BH900 headset and USB bluetooth dongle and I'd like to use them to make calls from a sofphone. I managed to this with boxe XTen-Lite and the Zoiper - but they both see the device as a simple sound card through the BlueSoleil drivers. While this is allmost usable, the headphone seems to be kept in "transmission" mode all the time and I get a constant "hissss" in the headset when I'm not actu...

Hi, I am facing a serious problem with file_column plugin. Though it is absolutely fine while uploading images, and the way it uploads the image on the file system and stores the file name in a database table is awsome. But while I try to get the image by passing its ID, surprisingly it returns the full absolute path rather than the ima...

Hi all, If I have a huge matrix/ dataframe and I want to create a new matrix/ dataframe with every second (or third, or fourth etc.) row of the original matrix, how can I do it? Any help much appreciated. Thanks, Umesh [[alternative HTML version deleted]]

hello everybody I use the following code for my programming & it runs with the error as specified below.Any help that would disolve the error will be highly appreciated. Thanks in advance my code looks like this #### R programme for simulating the power of the two sample t test vs various #### non-parametric alternatives sim.size <- 200 sample.size <- 10 set.seed(...

...ion about sapply and split combination. I have a big dataframe (40000 observations, 21 variables). First variable (factor) is "date" and it is in format "8.29.97", that is, I have monthly data. Second variable (also factor) has levels 1 to 6 (fractiles 1 to 5 and missing value with code 6). The other 19 variables are numeric. For each month I have several hunder observations of 19 numeric and 1 factor. I am normalizing the numeric variables by dividing val1 by val2, where: val1: (for each month, for each numeric variable) difference between mean of ith numeric variable in f...

Sorry for the repeat. The PRESS statistic is defined as sum(y-yhat(i))^2, where yhat(i) denotes the ith predicted value using all the data except the ith case (as used typically in linear models). Thanks again Jacob Jacob L van Wyk Department of Mathematics and Statistics Rand Afrikaans University P O Box 524 Auckland Park 2006 South Africa Tel: +27-11-489-3080 Fax: +27-11-489-2832 _____________...

Hello, fellow R enthusiasts. Ok, I've been racking my brain about this small issue, and between searching the help archives and reading through the plot-related documentation, I can't figure out how to achieve my desired endpoint without some ugly, brute force coding. What I would like to do is make a plot in which only a subset of my data are plotted, but in regular intervals, such as every 5th point along the sequence. Is anyone aware of a built-in function in plot or a related graphing family that can do this, or alternative...

Hi all! I have a list of matrices and I want to multiply the ith element of the list with the ith element of a another vector. That is, > LL <- list(A=diag(3),B=diag(3),C=diag(3)) > vec <- 1:3 > for(i in 1:3) + { + LL[[i]] <- LL[[i]]*vec[i] + } > LL $A [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 1 0 [3,] 0 0 1 $B...

...im. I figured one place to start would be some of the really cool applications that Asterisk has that -- generally, at least -- don't require licensing. Some of my favorites are follow-me, meetme, voicemail-to-e-mail and fax-to-e-mail. What are some of your favorite features/applications, be ith native or third-party? Thanks, -Ken -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

I have a df with a vector v. For each element of the vector, I want to know whether the i-2nd element is the same as the ith element. For example: given v=c(A,C,D,C) the result should be: FALSE,FALSE,FALSE,TRUE. I attempted something using indexing in a for loop such as (bad, incorrect example): for (i in v){ i...

Dear List I am having a bit of a problem getting a program to work. For each of i=1 to n persons I have a matrix (different for each person) with m rows. What I want to do, is create m new data sets such that the first is made up of the first row for each person from the original matrices, the second contains the second row for each person from the original matrices etc etc up to the mth new set containing the mth row from the original...

Hello, I have two vectors x1 and x2 both in increasing order. I want to select the x1[j]th entry which is the max min of the x2[i]th entry. I can do this using if and for statements but is there a quick way to do it without running a loop? Thank you in advance, Pantelis

Dear list; How can I speed up the run of following code (illustrative) #======================================================================== con<-vector("numeric") for (i in 1:limit) { if(matched data for the ith item found) { if(i==1) {con<-RowOfMatchedData } else {con<-rbind(con,matchedData)} } } #======================================================================== each RowOfMatchedData contains 105 variables, when "i&quot...

...o R code. I would like to use R to create the C[] matrix required for a car.proper model in WinBUGS, but I am having a difficult time negotiating the coding. The C matrix provides normalized weights for each pair of spatial areas. So the WinBUGS example is as follows: # of the weight matrix with elements Cij. The first J1 elements of the C[] vector contain the # weights for the J1 neighbours of area i=1; the (J1+1) to J2 elements of the C[] vector contain # the weights for the J2 neighbours of area i=2; etc. # To set up this vector, we need to define a variable cumsum, which gives the v...

...ike to generate 33 x 1000 newd values again using newa and newb values and the same c values, eg: Code #3 newa <-matrix(seq(16,24, by = sign(24-16)*0.25)) newb <-matrix(seq(24,36, by = sign(36-24)*0.25)) c <-runif(1000,50,60) newd <-exp(-c*newa)+exp(-c*newb)   Next I would like to sum d(ith) and newd(ith) according to their by sequence location in the matrices (or whatever type of objects we are working at this point) based on a and newa . As an example, Code #1 should generate a of 0,0.25,0.5, and so on with respective d for the a value 0, 0.25, 0.5, and so on..., while Code #3 shoul...