similar to: VEC Operator in R

Sender: r-help-bounces at r-project.org On-Behalf-Of: kingsfordjones at gmail.com Subject: Re: [R] Calculating confidence limits for the difference between means Message-Id: <2ad0cc110810231137n66760e43i61bdd0c89622d040 at mail.gmail.com> Recipient: ngottlieb at marinercapital.com -------------------------------------------------------- This information is being sent at the

I know, this is a forum about R. But I am so desperate of this problem (BTW, anyone knows any good Statistics/Math forum to post question like this?): A and B are both n x n positive definite matrix. Denote A > B, if A - B is positive definite. I know this is true: if A > B, then A^{-1} < B^{-1}. But how to prove this? I tried to diagonalize A and B, but since they can have different

Hello, I want to add text annotation about correlation on "pairs" plots. I found that I could pass a function to the "panel" argument of pairs : panel.annot <- function(x, y, ...) { points(x, y, ...) c <- cor.test(x, y) legend("topleft", legend=substitute(rho == r, list(r=sprintf("%.2f", c$estimate))), bty="n") } And then :

This might not be the correct forum for this question for there might be some flaws in my logic so the R function I'm looking for might not be the correct, but I know there?s a lot of smart people in this forum so please correct me if I'm wrong. I have been googling and searching in this forum for something useful but so far I'm out of luck. This is the background to my problem. I

Good Morning, I am using xyplot to show two variables for a large number of subjects. Currently I am using xyplot(y~x|as.factor(ID), layout=c(7,9)) to give me a little plot of the x and y factors for each ID. The purpose of the plot is to let the PI quickly look through the data for irregularities. The good news is that they like the layout of the xyplot, but they don't want to have to

odfWeave is throwing up a cryptic error, after successfully processing a file with a large number of figures (~30) and many \Sexpr{} calls. The error is (at least to me) cryptic. I am not sure where to look to correct this problem. Here's the error: Post-processing the contents Error in .Call("RS_XML_Parse", file, handlers, endElementHandlers, as.logical(addContext), :

hello R users, I didn't find a solution for a special problem. I have two dataframes. dataframe1: X value row col ID 1 8.973498062 5512625 3460000 1 2 11.656658570 5501625 3464000 2 3 11.121777570 5495625 3473000 3 4 9.310465964 5508625 3477000 4 5 8.883483845 5515625 3496000 5 dataframe2: X value

just to add to ted's explanation in case it helps to fix the email server problem: below are the two Recipients and one of the two is always contained in the repeated emails that I'm receiving. Also, of course i don't mean to claim that the names contained in below are doing anything to cause the problem . VEC Operator is just an example. There are 3 or 4 others that repeat

Hi folks, Both plotrix and prettyR have been updated. plotrix now has the "Matlab" style dotplot that was requested not long ago on the list, provided by that intrepid duo, Barry Rowlingson and Rolf Turner. Mike Cheetham has contributed an implementation of the Piper diagram, also in response to a request. prettyR has finally learned how to display SPSS value labels in the xtab and

Sender: r-help-bounces at r-project.org On-Behalf-Of: comtech.usa at gmail.com Subject: Re: [R] [R-SIG-Finance] forecasting earnings, sales and gross margin of a company... Message-Id: <b1f16d9d0810231239k506d582i7ecb908b84bc1642 at mail.gmail.com> Recipient: ngottlieb at marinercapital.com -------------------------------------------------------- This information is being sent at the

Using R 2.6.2, say I have the following list of lists, "comb": data1 <- list(a = 1, b = 2, c = 3) data2 <- list(a = 4, b = 5, c = 6) data3 <- list(a = 3, b = 6, c = 9) comb <- list(data1 = data1, data2 = data2, data3 = data3) So that all names for the lowest level list are common. How can I most efficiently access all of the sublist items "a" indexed by the outer

Hello everyone, Considering the following code sample : ---- indexes <- function(vec) { vec <- which(vec==TRUE) return(vec) } mat <- matrix(FALSE, nrow=10, ncol=10) mat[1,3] <- mat[3,1] <- TRUE ---- Issuing apply(mat, 1, indexes) returns a 10-cell list, as expected. Now if I do: ---- mat[1,3] <- mat[3,1] <- FALSE apply(mat, 1, indexes) ---- I would expect a

I need to process some datasets where the configuration information was stored in .INI-like files, i.e., text files with sections like this: [Section1] var1=value1 var2=value2 [Section2] A=value3 B=value4 ... >From Google and other searches I haven't found any package, or function within a package, that reads .INI files into an R list, or other data structure. Any suggestions, or

What should findInterval(x,vec,all.inside=TRUE) return when length(vec)<=1, so there are no inside intervals? R-3.3.0 gives a decreasing map of x->output when length(vec)==1 and -1's when length(vec)==0. Would '0' in all those cases be better? > findInterval(x=c(10, 11, 12), vec=11, all.inside=TRUE, rightmost.closed=FALSE, left.open=FALSE) [1] 1 0 0 >

Let say i have a square matrix and applied the 'vech' operator to stack the lower triangular elements into a vector: > Mat <- matrix(1:25, 5) > Mat ? ? ?[,1] [,2] [,3] [,4] [,5] [1,] ? ?1 ? ?6 ? 11 ? 16 ? 21 [2,] ? ?2 ? ?7 ? 12 ? 17 ? 22 [3,] ? ?3 ? ?8 ? 13 ? 18 ? 23 [4,] ? ?4 ? ?9 ? 14 ? 19 ? 24 [5,] ? ?5 ? 10 ? 15 ? 20 ? 25 > Mat[lower.tri(Mat)] ?[1] ?2 ?3 ?4 ?5 ?8 ?9 10 14

Hi everyone, I'm using the Johansen framework to determine a VEC using package urca. I have estimated the corresponding VEC using likelihood ratio test for restrictions on alpha, beta or both and I have generated objects of the class cajo.test. Now I want to diagnostic tests in the model, like heteroskedasticity test, residuals normality and serial autocorrelation, but I cannot find the way

Dear R-list members, Just wondering if there is any way to compute the duplication matrix in R. I tried to search for it but only found functions "xpnd" and "vech". Basically for a symmetric n by n matrix A, the duplication matrix D_n is a matrix of dimension n^2 by n(n+1)/2 such that D_n vech(A)= c(A), where c(A) just vectorizes A. The duplication matrix is defined on page

summary: Specifically, how does one do stack/FIFO operations in R? Generally, how does one code functions with side effects in R? details: I have been a coder for years, mostly using C-like semantics (e.g., Java). I am now trying to become a scientist, and to use R, but I don't yet have the sense of "good R" and R idiom (i.e., expressions that are to R what (e.g.) the Schwartzian

On Mon, Nov 16, 2020 at 12:55 PM David Blaikie <dblaikie at gmail.com> wrote: > I will say I'm not a huge fan of adding even more names for things in > this fairly core/common use case (now we'll have even more vector > names to pick from) - can we use default template arguments so we can > write SmallVector<T> instead of SmallVector<T, N> and would that >

I believe eigen(), svd() and qr() can all do it. Andy > From: Jonathan Baron > > On 06/17/04 19:04, Fred wrote: > >Dear R-listers: > > > >I am trying to test an algorithm on a set of linearly > independent vectors > >{x1,x2,...,xn}. > > Well, here's an idea, for 10 vectors of length 10, > as columns of a matrix m1. The 11th seems to be needed.