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
5 views

Arimax Prediction : Using Forecast Package

The arimax function in the TSA package is to my knowledge the only R package that will fit a ...
1
vote
0answers
4 views

Zero Inflated Poisson Regression in R

I have a zero inflated Poisson regression model I want to run R, as my data are over-dispersed and I have many zeros, but I am not sure how to set up the model in R and interpret the output. First, is ...
0
votes
0answers
6 views

how one would implement a “modified maximum likelihood regression” method

I am looking to extend upon the ideas of the following study. In the paper, the author describes a Modified Maximum Likelihood Regression as proposed by Nakahira in 2000. So to follow the method from ...
0
votes
0answers
5 views

Output of fixed effects summary in lmerTest in R and post-hoc tests

I'm doing a two-factor ANOVA using the lmerTest package. Each factor has multiple levels. When one (or more) of the effects are significant, I would like to do a post-hoc test to determine which of ...
1
vote
0answers
10 views

Stationarity in Forecasting Weekly Retail Sales

I'm new to forecasting and trying to create a model to forecast one step ahead weekly sales from my company. The variables I've identified for the model are Lag 1 Markdown spend and lag 1 sales, and ...
-1
votes
0answers
9 views

Using Pelora function in R [on hold]

Although I examined the usage of pelora function which is under supclust library, I can't get why it uses. Is there anyone that explains how it works and what pelora function is ?
0
votes
0answers
5 views

clear function, useMethod from result [migrated]

how delete this function (x, ...) UseMethod("as.matrix") <bytecode: 0x000000000f13a888> <environment: namespace:base> from my result? my code: ...
0
votes
1answer
14 views

Create clusters subject to constraint

I have a set of nodes. Each node represents a building and it has some attributes. For example: x coordinate y coordinate population What I want is to create ...
1
vote
0answers
42 views

Square root link function

I'm running a glm that estimates gaussian variable of production in kilogrames using different independent variables. i found a problem of heteroskedasticity so i tried different transfromations of my ...
0
votes
0answers
12 views

R survreg returns reciprocal of incidence rate ratios?

I'm working through some poisson regression for survival data using R. As I understand it, the exponentiated output from survreg with dist = "exponential" should give the incident rate ratios. I ...
2
votes
0answers
13 views

How to do forecasting with detection of outliers in R? - Time series analysis procedure and Method

I have monthly time series data, and would like to do forecasting with detection of outliers. This is the sample of my data set: ...
0
votes
0answers
19 views

Implementation of ROC Curve in my data?

I have successfully created 2x2 contingency table.The contingency table is created through the following process: ...
0
votes
0answers
18 views

Evaluating logistic regression and interpretation of Hosmer-Lemeshow Goodness of Fit

As we all know, there are 2 methods to evaluate the logistic regression model and they are testing very different things Predictive power: Get a statistic that measures how well you can predict ...
0
votes
0answers
12 views

Improving my knn prediction

I am still somewhat new to R and machine learning. I am trying to predict the cases of diabetes in the PimaIndiansDiabetes data set. My results were not the best- and I am wondering what I could do to ...
0
votes
1answer
22 views

Numerically Solve for Parameters Characterizing Lognormal RV that's Truncated from Above

I am trying to numerically solve for parameters characterizing a lognormal distribution truncated from above with first moment = mean, second moment = ...
1
vote
2answers
29 views

Beta distribution GLM with categorical independents and proportional response

My data is percentage disease data of different varieties of plants that had been inoculated with disease from several different sources. having conducted two-way ANOVA in SPSS (using the log10+1 of ...
0
votes
0answers
9 views

“Market” testing instead of A/B testing?

I believe in the last couple of years some respected research organization wrote a blog post about rigorous statistical testing of variables across markets (e.g., countries) when you can't do per-user ...
0
votes
1answer
36 views

Random Forest Overfitting R

I used a two-step cforest in my model. the accuracy of the train set is 87%, yet the accuracy of the test set is 57%. This indicates the model is severely overfitting. How to solve this problem? ...
2
votes
1answer
77 views

Mathematically,how are the critical t values calculated?

I am willing to write a program for calculating the critical t values,preferably in R. I am aware of the functions used for the same. But, I am interested in getting aware of the mathematical ...
1
vote
1answer
36 views

Material for teaching R to medical students [on hold]

I am planing to start teaching R to undergraduate medical students. I am looking for any material that can help in doing so. These students did not have any experience using R. Another question: Is ...
0
votes
0answers
15 views

Finding most relevant variables for components in PCA / MCA [on hold]

I am conducting PCA and MCA analyses in datasets with 150+ variables using Factominer in R. My main goal is to create 2-dimensional views of these datasets. After doing so, I would like to explain the ...
0
votes
1answer
34 views

Greek letters in R legend [on hold]

How can I add greek letters like "gamma" in the legend? My program: ...
0
votes
1answer
39 views

R function survfit(): give newdata reference values for factors with >2 levels?

The goal is to plot separate survival curves for each level of a categorical variable in R, using a coxph model including the variable and a second categorical variable with >2 levels. If the second ...
0
votes
0answers
21 views

Interpreting circular-linear regression coefficient

I'm trying to use the circular package in R to perform regression of a circular response variable and linear predictor, and I do not understand the coefficient ...
0
votes
0answers
8 views

Prediction using r mlogit: How to reproduce “Frequencies of alternatives” from model summary

I am trying to reproduce results in Train & Croissant's vignette "Kenneth Train’s exercises using the mlogit package for R." Specifically, Train illustrates multinomial-logistic regression using a ...
1
vote
1answer
25 views

Interpret Wilcoxon-Mann-Whitney Result in R [duplicate]

I am doing a Wilcoxon rank sum test (aka Wilcoxon-Mann-Whitney) with R. There are two variables "ScorePerView", that is metric and "HasCodeElement", that is binary. So the syntax is like that ...
0
votes
0answers
17 views

t.test applied in pairs to more than two samples in R

I have 44 samples where I only have its sample size, average and 1 standard deviation. I asked about the possibility of doing a t.test and some of you returned the answer: ...
0
votes
0answers
7 views

Estimation of scale parameter and shape parameter of using 'maxLik' package [on hold]

I want to test the closeness of my sample data with the Generalized Exponential (GE) distribution. For that, I am using ks test in R. In the documentation of ...
0
votes
1answer
57 views

Why do sample ACF/PACF suggest different TS models after box-cox transformation?

I use auto.arima function in R to fit a TS model to a annual data composed of electricity demand. The series is twiced difference to eliminate the trend in the data. After that the data is transformed ...
0
votes
1answer
4 views

How to judge, which random effects to include in a mixed model?

I am currently writing my master thesis about the effect of an insecticide (clothianidin) on the microflora of bumblebees. I received the bumblebees from an experiment with a nested study design. 16 ...
1
vote
2answers
27 views

How to calculate the probability of success of a Logistic Regression model with a single continuous predictor?

I have a logistic regression model below, predicting a dichotomous variable type from a single continuous predictor fatigue. Using the coefficients below I can obtain the increase in the odds of a ...
0
votes
1answer
17 views

Flag Random or Nonsensical Entries in a data set using R

This is a variation on my question Machine learning to catch fraud I have a data set with about 5 million rows The data set contains names and addresses of companies The name of the company is free ...
0
votes
0answers
17 views

Descriptive statistics in my native language in R [migrated]

Regardless of R! variables: Sys.setlocale("LC_MESSAGES", 'pl_PL.UTF-8') Sys.setenv(LANG = "pl_PL.UTF-8") when I use function such as summary, output is given in ...
0
votes
0answers
17 views

R for lme complicated study design

I have a data set where siblings (brothers and sisters) had the quantity of a specific bacteria measured at specific intervals since their birth (2 weeks, 1 month, 2 months, 6 months, 1 year). Because ...
0
votes
0answers
10 views

R - cor.ci results interpretation

I am doing some survey analysis for my dissertation, and am interested in studying possible correlations between four variables, based on Technology Acceptance Model.The survey is built on a ...
0
votes
0answers
8 views

How to check stability condition of VEC estimates in R?

I am estimating a VEC model and need to check the stability of its parameters. The vars package has a function to do this on an object of class varest generated by ...
1
vote
1answer
19 views

Different output for cox regression in R vs SPSS

During my analysis that I perform with a collaborator we get consistently different ouputs in cox regression between R and SPSS. Betas and p-values differ minutely, usually somewhere in the third or ...
0
votes
0answers
12 views
0
votes
1answer
13 views

Confidence intervals for Tobit model in package AER in R

I use a Tobit model to predict censored data. I use the AER package in R. A toy example looks as follows: ...
0
votes
0answers
12 views

How do I check whether my data fits exponential or lognormal or weibull using chi square and ks tests in R?

I'm a newbie to R language and I need to quickly find whether a set of data I have in a CSV file fits the three distributions Generalized exponential, log normal and Weibull. I need to test it using ...
0
votes
0answers
21 views

How do I use MatLab or R data model in an machine learning application [on hold]

I want to use matlab or R to create a prediction model, once this is complete how would I use that model in a an application?
0
votes
0answers
12 views

How to extract values for only one particular fixed factor from lmer object in R? [on hold]

I have following output for lmer object constructed using coxme package: ...
0
votes
0answers
11 views

Correlation analysis and Clustering in R - Some suggestions [on hold]

I have a biological data, comprising of intensity values associated to 24 masses. Each mass file has 10,000 intensity values corresponding to 10,000 coordinates. So, the data matrix in R would look ...
-2
votes
0answers
14 views

Bayesian data analysis with R [on hold]

Lets say I have a big list of phone models I have bought through the years with several features listed as a "checkbox"-style: ...
0
votes
1answer
21 views

MLE of Gamma when only given observations [on hold]

i've been given 10 observations of X, where X is a random variable. the observations are 141 16 46 40 351 259 317 1511 107 567 and now given they are gamma distributed, how do you find the MLE using ...
1
vote
1answer
18 views

Bland & Altman plot for repeated measures using one measuring device

I am examining the reliability of a measuring device by measuring a force generated på test persons 7 days apart using R. In order to assess this properly I want to make a Bland & Altman plot ...
1
vote
0answers
51 views

How to detect different social groups in a small questionnaire dataset?

I work with a very small questionnaire dataset, in which perception of 10 respondents to the size of, say, an unseen population is quantified: The population is divided into eight sub-populations, and ...
0
votes
0answers
18 views

R ggplot2: overlaying multiple geom_ribbon objects in a single plot using a nested loop [migrated]

I have five line plots from which I'd like to output a shaded area that represents the region between their plotted upper and lower regions. I'm creating an R script (see below) as I have multiple ...
0
votes
0answers
9 views

How can I determine a transition point in a curve?

I have plotted a bar graph and knowing that the transition in the curve takes place around values 21 to 25 , I want to find such point in general. ...
0
votes
1answer
27 views

Comparing Height of Web Page against Scroll %

I have a collection of 108 data points in the following format: page height | % of users who scrolled upto 25% of page | 50% | 75% | 100% (full page) I'm trying ...