search for: 0.0211

I have two matrices (see example below) and I want the differences for the matching row numbers, but the row numbers are not identical in the two matrices. There are probably many ways to do this. Anyone know of any easy way to do this? I could loop over them, but you know what they say about for loops... Thanks, Roger > out.r[1:5, 1:3] 1 2 3 1100 -0.0992 -0.0802 -0.0653 1200 -0.1242

Hi I have a very big data.frame : str(fslu12) 'data.frame': 277200 obs. of 11 variables: $ V1 : num 304 304 304 304 304 ... $ V2 : num 300 300 300 300 300 300 300 300 300 300 ... $ V3 : num 10 10 10 10 10 10 10 10 10 10 ... $ V4 : num 79.1 52 30.6 37.3 63 79.1 52 30.6 37.3 63 ... $ V5 : num 0.005 0.00498 0.00496 0.00491 0.00484 ... $ V6 : num 0.01 0.00997 0.00988 0.00969

Hello All, Thank you for taking my question. I am looking for information on how R handles interaction terms in a multiple regression using the ?lm? command. I originally noticed something was unusual when my R output did not match the output from JMP for an identical test run previously. Both programs give identical results for the main test and if the models do not contain the interaction

Hi, I'm using metafor to perform a mixed-effects meta-analysis. I'd like to test whether the effect is different for animals and plants/whether "group" (animal/plant) influences the effect size, but am having trouble interpreting the results I get. I've read previous posts about QM in metafor, but I'm still a bit confused. I've dummy-coded the factors:

Hi, I'm sure this is simple, but I haven't been able to find this in TFM, say I have some data in R like this (pasted here: http://pastebin.com/raw.php?i=sjS9Zkup): > head(df) gender age smokes disease Y 1 female 65 ever control 0.18 2 female 77 never control 0.12 3 male 40 state1 0.11 4 female 67 ever control 0.20 5 male 63 ever state1 0.16

Hello! I'm puzzled by the following problem. It occurs while trying to predict responses in a test-dataset using a linear model fitted with regFit from the rMetrics "fRegression"-package. All goes well when I call "predict" using the training dataset. However, a call using the test-dataset retuns an error message - telling me that the latter dataset provides variables

Hello, I have used R in the past to conduct multiple comparisons on standard linear models, but am a bit confused as to how to go about doing it with a mixed effects model. I am conducting a bioindication study using carabid beetles in which I have four treatment types (forest harvest types with varying levels of canopy structure retention), and am using canopy closure percent as a covariate in

Hi all I have another GLIM question. I have been using R as well as Genstat (version 6) in order to fit GLIM models to the data (displayed below). The same models are fitted but the answers supplied by the two packages are not the same. Why? Can anyone help? A discription of the data and the type of model/s fitted can be found below. Regards Allan The

Dear all, I'm using R 2.8.1 under Windows Vista on a dual core 2,4 GhZ with 4 GB of RAM. I'm trying to reproduce a result out of "Analysis of Financial Time Series" by Ruey Tsay. In R I'm using the fGarch library. After fitting a ar(3)-garch(1,1)-model > model<-garchFit(~arma(3,0)+garch(1,1), analyse) I'm saving the results via > result<-model

Hello R-users, I have been trying to fit what I think is a simple mixed-effects model using lmer (from lme4), but I've run into some difficulty that I have not been able to resolve using the existing archives or Pinheiro and Bates (2000). I am measuring populations (of birds) which change with time at a number of different sites. These sites are grouped into regions. Sites are not measured

Dear all R-users, I wanted to fit a Garch(1,1) model to a dataset by: >garch1 = garch(na.omit(dat)) But I got a warning message while executing, which is: >Warning message: >NaNs produced in: sqrt(pred$e) The garch parameters that I got are: > garch1 Call: garch(x = na.omit(dat)) Coefficient(s): a0 a1 b1 1.212e-04 1.001e+00 1.111e-14 Can any one

Hi, I recently changed from uw imap to dovecot on the sound recommendation of a friend and have mostly succeeded in getting all of my clients up and running, but am really stuck with the iPhone which is failing to make connections. I run certificates on all of my clients and thunderbird happily connects both locally and remotely. I installed the certificate on the iPhone after great pain (pk12

I've searched the Web with Google and do not find what might cause this particular error from an invocation of cenboxplot: cenboxplot(cu.t$quant, cu.t$ceneq1, cu.t$era, range=1.5, main='Total Recoverable Copper', ylab='Concentration (mg/L)', xlab='Time Period') Error in if (n > 0) (1L:n - a)/(n + 1 - 2 * a) else numeric() : argument is of length zero I do

On Fri, Sep 14, 2007 at 11:42:18PM -0700, Tanya Lattner wrote: > The 2.1 pre-release (version 1) is available for testing: > http://llvm.org/prereleases/2.1/version1/ > > [...] > > 2) Download llvm-2.1, llvm-test-2.1, and the llvm-gcc4.0 source. > Compile everything. Run "make check" and the full llvm-test suite > (make TEST=nightly report). > > Send

Are these 225 compile time regressions real? It sure looks bad! Ciao, Duncan. On 01/12/11 09:39, llvm-testresults at cs.uiuc.edu wrote: > > bwilson__llvm-gcc_PROD__i386 nightly tester results > > URL http://llvm.org/perf/db_default/simple/nts/380/ > Nickname bwilson__llvm-gcc_PROD__i386:4 > Name curlew.apple.com > > Run ID Order Start Time End Time > Current 380

LLVMers, The 2.1 pre-release (version 1) is available for testing: http://llvm.org/prereleases/2.1/version1/ I'm looking for members of the LLVM community to test the 2.1 release. There are 2 ways you can help: 1) Download llvm-2.1, llvm-test-2.1, and the appropriate llvm-gcc4.0 binary. Run "make check" and the full llvm-test suite (make TEST=nightly report). 2) Download