Use this tag for any *on-topic* question that (a) involves `R` either as a critical part of the question or expected answer, & (b) is not *just* about how to use `R`.

learn more… | top users | synonyms

1
vote
0answers
9 views

How to predict & plot quasipoisson GLM in R

I have a set of complex survey data with sampling weights. I am using the svyglm() function from the survey package in R to ...
0
votes
1answer
11 views

How to calculate distance between 2 lists in r with only general attribute information?

I need to calculate distance between two simple lists in R and would like a little logic help please. The issue is that I do not have a data set, rather just a range for each attribute. Here's what I ...
1
vote
1answer
12 views

Incorrect class displayed when using 'apply' function in R

I see a problem with the apply function in R. 1. This does not show the right class. ...
1
vote
1answer
12 views

What does the acf of the log of data tell us?

I'm working with time series data. And there is a use of acf ( auto correlation function ) being applied to the log of a ...
0
votes
0answers
3 views

How to plot identical slopes for an ANCOVA (1 categorical variable(2 levels) with 2 covariates) using ggplot2

I am trying to compare and plot regression lines of an ANCOVA when their slopes are the same and no interaction effect There are two sets of ANCOVAs plots that I am trying to create: a) 1 ...
0
votes
0answers
9 views

R and Eviews differences in ADF test

So I've implemented the ADF test in a large number of variables using eviews and R (urca package, ur.df). Althought the criteria used to select lags and the deterministic variables included were the ...
1
vote
1answer
28 views

What is the most appropriate way to transform proportions when they are an independent variable?

I thought I understood this issue, but now I'm not as sure and I'd like to check with others before I proceed. I have two variables, X and ...
0
votes
0answers
13 views

Metafor package: rma - How to get agreement with metan in stata when the confidence intervals are not symetric about the effect estimate

I would be extremely grateful for help with the following. I perform meta-analyses using effect estimates and confidence intervals. However, I can only get agreement between rma in R and metan in ...
1
vote
0answers
8 views

R Caret - Repeated CrossValidation, finalModel and ROC curves

I got a problem understanding the meaning of the finalModel when using a repeated CV. ...
0
votes
1answer
11 views

Test same coefficient from two different models in R

I am working with two ivreg models. The only difference between them is the choice of instruments on the other side of the line. Example: ...
0
votes
0answers
21 views

question about GBM [on hold]

...
1
vote
0answers
4 views

How to implement a SVAR with sign restrictions

I am trying to estimate a bi-variate sign-restricted SVAR with daily oil and stock prices and two shocks (demand and supply). The ultimate goal is to explain how much of the recent fall in oil ...
0
votes
0answers
13 views

If statement doesnt work [on hold]

When I run my code, then nothing happens. My data isnt modified as it should. My code: ...
0
votes
1answer
33 views

Multiple Correspondence / Correlation Analysis

I'm a little bit new to statistics, so I'm not sure what I really need. I have a table like the following with some categorical/nominal values (like Gender and Age Group) and a ratio scaled value ...
-1
votes
0answers
9 views

How to create a data subset by using a dynamic column name / variable? [on hold]

Why are dynamic column names not working in brackets? ...
0
votes
0answers
14 views

Statistical analysis for interaction of genes in multisamples

Say I have 120 human individuals. I have done gene-gene interaction in those individuals and found the frequency of all the interacting genes as shown in the data below. I want to do some analysis to ...
0
votes
0answers
15 views

Random effects in GAM [on hold]

I am now trying to use random effects in GAMs developed by Professor Simon Wood. Prof Wood uses s(...,bs="re") to account for the random effects. Random intercepts models or random slopes models are ...
1
vote
1answer
24 views

How to calculate confidence interval using parametric bootstrap

I was able to calculate the age-specific hazard ratio (HR) for breast cancer at age 40 and 60 years old. The HR was a continuous, piecewise linear function of age which was constant before 40 years, ...
0
votes
0answers
16 views

Clustering Matrix with binary values

I have a large matrix. It consists of about 10.000 rows (each row one document) and 10.000 columns (each column one word). The binary value indicates if a word exists (1) in the particular document or ...
0
votes
0answers
18 views

Face detection in R [on hold]

I'm trying to get started with this area (face detection) and I see a lack of information with R language. Are there any example or algorithm for face detection in ...
0
votes
0answers
19 views

beta regression creating a wild residual vs. fitted plot- whats going on?

I recently ran a beta regression model in R using the betareg package. I am modeling a continuous dependent variable (a fraction out of 1) that is bound between 0 ...
0
votes
1answer
22 views

R & propensity matching [on hold]

I am new to R & propensity matching. I doing a quasi-experimental study in which a group of 320 students received treatment and a group of 280 students did not receive treatment. I would like ...
0
votes
0answers
15 views

add legend to ggplot geom points [on hold]

I've created a ggplot indicating 3 different activities for male bats in three different colours: ...
0
votes
1answer
18 views

Negative BIC in k-means

Probably a simple question but I'm trying to interpret BIC for k-means. I have some k-means clustering and calculating BIC gives me a negative value, with a plot something like this: ...
0
votes
0answers
16 views

How to export my residual data? [on hold]

I am having a lot of trouble simply trying to export my residual data. I'm new to R so not sure what I'm doing wrong. I've tried doing library(xlsx) write.xlsx(Residuals, "c:/ResidualsData.xlsx") ...
0
votes
0answers
10 views

Is it possible to use color.scale to create color graded polygons in a polar.plot with plotrix? [on hold]

I have been trying to map sound levels in multiple directions from a single sound source. I have average dB readings from 45 degree intervals around the source. I have plotted these using the ...
0
votes
0answers
10 views

What is the best way to write an R data frame to a SQL Server database [migrated]

I have a data frame that I am trying to write to a table on SQL Server. The data frame contains only 8 of about 12 columns contained in the table and the columns in the data frame are not in the same ...
0
votes
0answers
13 views

TBATs - setting seasonal periods

I am looking to do some forecasting with multiple seasons and after a lot of research it looks like I'm going to test 4 different approaches. The first two will be Arima. One with xreg = dummy where ...
0
votes
0answers
3 views

Extracting significance score from AER coeftest results in R [migrated]

Calling coeftest(m) returns z test of coefficients: ...
-2
votes
0answers
18 views

How to apply a running function in R [on hold]

I need to develop a function in R which can be used to compute time varying variable P which is defines as: p(t)=k(t)*p(t-1)+j(t). ...
0
votes
0answers
23 views

How to update R version using code instructions? [on hold]

I want to write a program which automatically updates R to the latest version. Is there such instruction I can use for this task?
-1
votes
0answers
29 views

Wishart PDF for 1000 x 8 [on hold]

I have a dat aset which is composite of 1000 row and 8 column , how can i calculate the wishart pdf for covariance matrix ? what are the argument for the Wishart pdf in this case ?? ...
0
votes
0answers
13 views

Cross-validation results differ from package

I am trying to test out a hand-rolled cross-validation procedure for ridge regression. I've run the glmnet package which gives me an MSE bottoming out around ...
0
votes
0answers
11 views

What is the difference between the “+” operator in ggplot2 and the “%>%” operator in magrittr? [migrated]

What is the difference between the "+" operator in ggplot2 and the "%>%" operator in magrittr? I was told that they are the same, however if we consider the following script. ...
1
vote
0answers
23 views

How to estimate ratio of variances?

Suppose we have two time series generated by \begin{align} x_t &= \sigma_t\,e_t \\ y_t &= c\,\sigma_t\, f_t \end{align} where $e_t$ and $f_t$ are IIDs with mean zero and unit variance. We ...
1
vote
0answers
33 views

Testing 3-way interaction of categorical variables in a meta analysis in R

I'm currently working on a meta-analysis using the metafor package for R. I have several categorical moderators, a few of which are multi-category, that I'm including in the analyses. When using one ...
2
votes
0answers
30 views

Fit a scatter plot spline function

I am trying to fit a cubic spline through some data, while trying to reproduce results found in Generalized Additive Models: an introduction with R by S.Woods p125 (pg 138 of pdf), but I am not sure ...
0
votes
0answers
26 views

Dickey-Fuller test: conflicting results between SAS and R

I'm doing a Dickey-Fuller test on Airplane data in R. I'm using adf.test in "tseries" package and also doing this test in SAS. But SAS calculates p-value 0.2 and R ...
0
votes
0answers
6 views

McFadden pseudo R squared on GLM.NB

I have used the McFadden pseudo R squared on a negative binomial glm (family=poisson). I was wondering if it's appropriate. I have found informations suggesting McFadden is useful for poisson glm but ...
0
votes
0answers
12 views

Bonferroni correction after hypergeometric test

I have a small doubt regarding a Bonferroni after running an hypergeometric test on genomic data. Let say I have a population of 379876 output data(integer numebers), and this population is clustered ...
0
votes
0answers
4 views

How to compare two incidence ratios statistically- the rate parameter test

I want to compare two standardized incidence ratios (SIR) obtained in stratified analyses. I read in Mocci et al (2013) that they used the rate parameter test and cited the book Statistical models in ...
1
vote
0answers
8 views

What is the appropriate statistical test to perform a difference of means within a difference of means?

Let's say I am interested in measuring the effect of three different types of interventions on test scores. Let's say the intervention is book quality, of which there are three kinds of book quality ...
1
vote
2answers
33 views

auto.arima not giving the best model according to information criteria

I'm using auto.arima to get the best model for the MASS dataset deaths. However, ...
0
votes
0answers
14 views

How do I compute the posterior distribution of the success probability with uniform priors?

This is a question from a homework assignment A clinical trial is conducted to compare the effectiveness of three drugs. 100 patients are randomly assigned to each drug (300 total patients), and Y_1 ...
0
votes
0answers
34 views

Statistics help required [on hold]

I am doing a project in A/B testing.I want statistical approach to evaluate my work. A/B testing is comparing two versions of a web page/ taglines /landing pages to see which one performs better. You ...
0
votes
0answers
26 views

How to fit a curve for a histogram [on hold]

I am trying to fit a curve (gaussian) for my distribution represented here with an histogram, like the red curve represented in the attached shema. X = rchisq(1000, df = 5) hist(X) Thank you
0
votes
0answers
11 views

Planned contrasts in ezANOVA

I have a data set that is 3 x 5 and I can easily run a repeated measures ANOVA using ezANOVA in R. I now wish to have planned contrasts going into the ANOVA. The data consists of one control ...
1
vote
1answer
60 views

Statistical significance in time series (R)

Stats newb here, I have to determine if two time series are really different instead of being part of the same population with noise in the samples. The data is a comparison between two algorithms ...
0
votes
1answer
16 views

SOM (Kohonen) using the term document matrix [on hold]

Language: R Package: kohonen Function: som I have a term document matrix (tdm) with 64 terms (row) and 1017 documents (columns). I want to use the self-organized-map to cluster the terms on 25 ...