similar to: Calculating confidence limits for the difference betweenmeans

Can anyone please tell whether there is any R function to act as "VEC" and "VECH" operator on Matrix? Yes of course, I can write a user-defined-function for that or else, I can put dim(mat) <- NULL. However I am looking for some R function. Your help will be highly appreciated. Regards, -- View this message in context:

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

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

Hi, I would like to use a weighted lm model to reduce heteroscendasticity. I am wondering if the only way to generate the weights in R is through the laborious process of trial and error by hand. Does anyone know if R has a function that would automatically generate the weights need for lm? Thanks, -- Tom [[alternative HTML version deleted]]

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 :

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

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

Hi, I need to correct for ar(1) behavior of my residuals of my model. I noticed that there are multiple gls models in R. I am wondering if anyone has experience in choosing between gls models. For example, how should one decide whether to use lm.gls in MASS, or gls in nlme for correcting ar(1)? Does anyone have a preference? Any advice is appreciated! Thanks, -- Tom [[alternative HTML

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

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

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

resending to provide a more informative subject line.... On Fri, Jan 2, 2009 at 3:21 PM, Kingsford Jones <kingsfordjones at gmail.com> wrote: > Hi David, > > A general answer to your question is: yes, R would be useful for such > analyses - particularly when interfaced with a GIS. For an > introduction to the types of spatial tools available in R see the Task > View for

R 2.8.1 Windows XP I am trying to analyze repeated measures data (the data are listed at the end of this Email message) and I need help to make sure that I have properly specified my model, and would like to know why lmer does not return a p value for Group, my fixed effect. My subjects are divided into two groups (variable GROUP), individual subjects are indicated by the variable SS, Value is

I hope that Marc doesn't mind, but I felt that part of his recent post was important enough to deserve it's own subject line rather then being lost in a 60-msg-long thread... On Sun, Jan 11, 2009 at 10:08 AM, Marc Schwartz <marc_schwartz at comcast.net> wrote: ... I strongly believe that the comments regarding R and the FDA are overly negative and pessimistic. The hurdles to

Colleagues, I am working on a problem at the edge of my knowledge of statistics. Basically, I have values for two groups - I am trying to calculate the 90% confidence limits for the difference between means. I can implement this without difficulty based on standard equations that are available in stat textbooks: (difference between means) / (pooled standard error) My problem arises

I just had a manuscript returned with the biggest problem being the analysis. Instead of using principal components in a regression I've been asked to analyze a few variables separately. So that's what I'm doing. I pulled a feather from young birds and we quantified certain aspects of the color of those feathers. Since I often have more than one sample from a nest, I thought I

Hello, I was wondering if anybody could help me with this? I have plotted an acf function for a time series and am very happy with it. Now I am interested in calculating for myself the two values for the confidence intervals that are plotted on the graph of the acf. The confidence intervals do not appear to be returned from the acf function (is this true?). So far I haven't managed to

Dear all, Myself Vijaykumar Muley working as senior research fellow. By training I am a computational biologist with not a strong knowledge of statistics. I have done some analysis which is explained as follows, I have 10340 (X) profiles of binary vectors with same length(N=845), I will call then "gene profiles" for example... v1 v2 v3 v4.....vN a 1 0 1 0 1 b 0

On 10/19/2016 01:54 PM, m.roth at 5-cent.us wrote: > Alice Wonder wrote: >> On 10/19/2016 11:34 AM, Leonard den Ottolander wrote: >>> Hello Gordon, >>> >> *snip* >>> >>> Personally I would be more concerned whether or not to enable ECDSA >>> algorithms (https://blog.cr.yp.to/20140323-ecdsa.html). >>> >> For web server ECDSA