search for: 0.482

I have a data frame, lets call it dat, with 3 columns ( mc, yr, ret) which represent market cap, year, and return. mc is a factor, mc, and ret are real numbers. I want to add a column to the data calculated as follows. For each year, I want to split the data by mc decile, then calculate the mean ret within that mc decile, and finally subtract that year's decile mean from the raw return. Then

Hola a tod@s, Hace poco en la lista se discutieron algunas aproximaciones para determinar la distribución de probabilidad que potencialmente podrían haber generado los datos (ver [1]) y una de ellas fué el AIC. Haciendo uso del programa enviado por Pablo Verde (ver [1]), estos son los resultados para mis datos: weibull , AIC = 69839.44 exponential , AIC = 79488.77 gaussian , AIC = 69413.03

Dear R community, I'm currently using R to analyze functional Magnetic Resonance Imaging data. Each analysis involves running ~120,000 repeated-measures ANOVAs. I would like to know if there is any automatic way to access the F- and P-value data that are associated with each of these 120,000 ANOVAs. For example, if the summary output (for the 1st ANOVA of 120,000) shows the following value

Dear R community members, I think I am asking a very simple question, but I really looked up in the faqs and manuals and found nothing helpful. I am trying to plot a data frame with the following structure (this is just a small extract): glo conc odor line series X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 1 0 AIR LN1 UP -0.488

I don't understand which is the problem Could you help me please? Thanks in advance Marta > str(data_indiciN2) 'data.frame': 200 obs. of 36 variables: $ gruppo: Factor w/ 2 levels "0","1": 2 2 2 2 2 2 2 2 2 2 ... $ I001N2: num 19.32 8.22 28.35 7.24 14.7 ... $ I002N2: num 2.92 2.54 0.11 1.6 7.12 ... $ I003N2: num -22.362 -0.222 -19.291

Dear forum members, What is the formula to calculate denominator degrees of freedom (den df) for nonlinear mixed-effect models with covariates? My model is similar to a CO2 uptake example from Pinheiro and Bates (2000, page 376). In this CO2 dataset, there are two treatments and two types (84 observations in total), but den df for each parameter of the model is 64. Isn’t it too high? Your

On Wed, 13 Dec 2023 09:04:18 +0100 Hilmar Berger via R-devel <r-devel at r-project.org> wrote: > Still, I feel that default partial matching cripples the functionality > of data.frame for larger tables. Changing the default now would require a long deprecation cycle to give everyone who uses `[.data.frame` and relies on partial matching (whether they know it or not) enough time to

Hello, could anyone explain what am I doing wrong. When I use colldiag function from package perturb I get different Variance Decomposition Proportions matrix in R than in SAS, although the eigenvalues and indexes are the same. Thanks for your attention. Results: in R: eigen(cor(indep2)) $values [1] 4.197131e+00 6.674837e-01 9.462858e-02 4.070314e-02 5.323022e-05 colldiag(indep2,c=T)

Dear All, I am trying to calculate a 95% confidence interval for the difference in two c statistics (or equivalently D statistics). In Stata I gather that this can be done using the lincom command. Is there anything similar in R? As you can see below I have two datasets (that are actually two independent subsets of the same data) and the respective c statistics for the variables in both cases.

Dear All, Apologies if you have a seen a question like this from me before. I am hoping that if I re-word my question more carefully someone may be able to offer more help than the last time I asked something similar. I am using R 2.9.2 and Windows XP. I am trying to determine if there is a statistically significant difference between two c-statistics (or equivalently D statistics). In Stata

I have used the greedy.wilks to stepwise discriminant analysis, but it doesn't work with my dataset. I don't understand which is the problem Could you help me please? Thanks in advance Marta > str(data_indiciN2) 'data.frame': 200 obs. of 36 variables: $ gruppo: Factor w/ 2 levels "0","1": 2 2 2 2 2 2 2 2 2 2 ... $ I001N2: num 19.32 8.22 28.35 7.24

Hi Hilmar and Ivan, I have used your code examples to write a blog post about this topic, which has figures that show the asymptotic time complexity of the various approaches, https://tdhock.github.io/blog/2023/df-partial-match/ The asymptotic complexity of partial matching appears to be quadratic O(N^2) whereas the other approaches are asymptotically faster: linear O(N) or log-linear O(N log N).

In addition to the major points that others made: if you care about speed, don't use compression. With today's fast disks it's an order of magnitude slower to use compression: > d=lapply(1:10, function(x) as.integer(rnorm(1e7))) > system.time(saveRDS(d, file="test.rds.gz")) user system elapsed 17.210 0.148 17.397 > system.time(saveRDS(d,

Dear All, I'm not sure of the interpretation of interactions with contrasts. Can anyone help? I do an ANCOVA, dryweight is covariate, block and treatment are factors, c4 the response variable. model<-aov(log(c4+1)~dryweight+treatment+block+treatment:block) summary(model); Df Sum Sq Mean Sq F value Pr(>F) dryweight 1 3.947 3.947 6.6268 0.01076 *

Dear all, I'm working on some data from an experiment on the breeding behavior of birds. In short, I have been measuring how the time spent on performing a certain task (variable 'mean_on_active') differs over time (variable 'day', 2 levels) across three experimental categories (variable 'treat'; levels 'C', 'R', 'E'). The model shows a

Deepayan Sarkar and I have a source package of release candidate 1 of the 0.6 series of the lme4 package available at http://www.stat.wisc.edu/~bates/lme4_0.6-0-1.tar.gz This package requires Matrix_0.8-6 which has been uploaded to CRAN and should be available in a few days. A copy of the source package is available as http://www.stat.wisc.edu/~bates/Matrix_0.8-6.tar.gz

Deepayan Sarkar and I have a source package of release candidate 1 of the 0.6 series of the lme4 package available at http://www.stat.wisc.edu/~bates/lme4_0.6-0-1.tar.gz This package requires Matrix_0.8-6 which has been uploaded to CRAN and should be available in a few days. A copy of the source package is available as http://www.stat.wisc.edu/~bates/Matrix_0.8-6.tar.gz

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

This is a summary and extension of the thread "GLMM (lme4) vs. glmmPQL output" http://maths.newcastle.edu.au/~rking/R/help/04/01/0180.html In the new revision (#Version: 0.4-7) of lme4 the standard errors are close to those of the 4 other methods. Thanks to Douglas Bates, Saikat DebRoy for the revision, and to G?ran Brostr?m who run a simulation. In response to my first posting, Prof.

Hi, I am dealing with very large datasets and it takes a long time to save a workspace image. The options to save compressed data are: "gzip", "bzip2" or "xz", the default being gzip. I wonder if it's possible to include the pbzip2 (http://compression.ca/pbzip2/) algorithm as an option when saving. "PBZIP2 is a parallel implementation of the bzip2