similar to: Models for Discrete Choice in R

Hi, I'm trying to fit a model using variance components in R, but if very new on it, so I'm asking for your help. I have imported the SPSS database onto R, but I don't know how to convert the commands... the SPSS commands I'm trying to convert are: VARCOMP RATING BY CHAIN SECTOR RESP ASPECT ITEM /RANDOM = CHAIN SECTOR RESP ASPECT ITEM /METHOD = MINQUE (1) /DESIGN

I would like to buy a basic statistics book (experimental design, sampling, ANOVA, regression, etc.) with examples in R. Or download it in PDF or html format. I went to the CRAN contributed documentation, but there were only R textbooks, that is, textbooks where R is the focus, not the statistics. And I would like to find the opposite. Other text I am trying to find is multivariate data

Hi, I would like to have an index for a column in a matrix encoded in a cell of the same matrix. For example: x = matrix(c(11,12,13,1, 21,22,23,3, 31,32,33,2),byrow=T,ncol=4) In this case, column 4 is the index. I then access the column specified in the index by: > for (i in 1:3) print(x[i,x[i,4]]) [1] 11 [1] 23 [1] 32 > > for (i in 1:3) {x[i,x[i,4]] <- x[i,x[i,4]] + 5} > x

Hi, Why sum() on a 10-item vector produces a different value than its counterpart on a 2-item vector? I understand the problems related to the arithmetic precision in storing decimal numbers in binary format, but shouldn't the errors be equal regardless of the method used? See my example: > options(digits=22) > x=rep(.1,10) > x [1] 0.10000000000000001 0.10000000000000001

Hi, I am trying to learn time series, and I am attending a colleague's course on Econometrics. However, he uses e-views, and I use R. I am trying to reproduce his examples in R, but I am having problems specifying a AR(1) model. Would anyone help me with my code? Thanks in advance! Reproducible code follows: download.file("https://sites.google.com/a/proxima.adm.br/main/ex_32.csv

Hi, I am not sure if R-help is the right forum for my question. If not, please let me know. I have to do some discriminant validity tests with some constructs. I am using the method of doing a CFA constraining the correlation of a pair of the constructs to 1 and comparing the chi-square of this constrained model to the unconstrained model. If the chi-square difference is not significant, then I

Hi, I want to run t.test() for several variables among two groups, and I would like to skip the tedious process of collecting information to assemble a table, but I am not sure if the function I want already exists. Any suggestion would be appreciated. I have a working example, as required by the posting guide: my_swiss = swiss[-1,] my_swiss$facto = rep(1:2,nrow(my_swiss)/2)

Hi, I don't know much about the RGL package, and I have read the documentation and tried some parameters, with no luck... I would like to generate a movie from a 3D object (code below), where the vortex A is closer to the observer, and then the object rotates and the B vortex gets closer. I would like to capture this movie to a file. By the way, I am not being able to insert unicode text

Hi, I am trying to compare two 3D plots. For that, I am trying to use the "pch" parameter in the "points3d" function, but it is not working. Is it implemented? Any suggestion? Here goes a reproducible code. I wanted the second plot having different symbols for the points. x <- rbind(matrix(rnorm(100, sd = 0.3), ncol = 2), matrix(rnorm(100, mean = 1, sd = 0.3),

Sorry for the cross-posting, but I would like to know if anyone is aware of a package in R for this. ---------- Forwarded message ---------- From: John Antonakis Sent: Sunday, April 10, 2011 3:26 PM To: RMNET Subject: Meta-analysis of a correlation matrix (correct thread title) Hi: Does anyone know of good program that can do a meta-analytic multiple regression (with multiple correlated

Hi, I'm doing a k-means cluster with 6 clusters and 15 variables. Any suggestions on how to plot the results? I've tried the standard xy plot, but couldn't get much of it. Thansk in advance, Iuri.

Dear netters, Sorry for cross-posting this question. I am sure R-Help is not a research methods discussion list, but we have many statisticians in the list and I would like to hear from them. Any function/package in R would be able to deal with the problem from this researcher? ---------- Forwarded message ---------- From: Heidi Bertels Date: Tue, Jun 5, 2012 at 4:31 PM Subject: How to best

Hello, Thanks for the reply Rui and for pointing out that I forgot to attach my code. Please find attached in this email my code and data. Thanks in advance. Best regards, Pedro Gerhardt Gavronski. On Fri, Feb 23, 2024 at 5:50?AM Rui Barradas <ruipbarradas at sapo.pt> wrote: > > ?s 16:34 de 22/02/2024, Pedro Gavronski. escreveu: > > Hello, > > > > My name is Pedro

The data came through fine, the program was a miss. Can you paste the program into a ".txt" document like a notepad file and send that? You could also paste it into your email IF your email is configured to send text and NOT html. TIm -----Original Message----- From: R-help <r-help-bounces at r-project.org> On Behalf Of Pedro Gavronski. Sent: Friday, February 23, 2024 5:00 AM To:

Hello, I am working in a longitudinal study, with a categorical variable as a dependent variable (alcohol consumption: no use, use, abuse and dependence) with repeated measures (baseline, 1 year, 2 years). Besides I have another variable with two groups: control and experimental. I would like to analyze the evolution in each group, I have though in lme, but I am not sure if I can do an lme with

Three related questions on LMEs and GLMMs in R: (1) Is there a way to fix the dispersion parameter (at 1) in either glmmPQL (MASS) or GLMM (lme4)? Note: lme does not let you fix any variances in advance (presumably because it wants to "profile out" an overall sigma^2 parameter) and glmmPQL repeatedly calls lme, so I couldn't see how glmmPQL would be able to fix the dispersion

?s 16:34 de 22/02/2024, Pedro Gavronski. escreveu: > Hello, > > My name is Pedro and it is nice to meet you all. I am having trouble > understanding a message that I receive when use function ar.ols from > package stats, it says that "Warning message: > In ar.ols(x = dtb[2:6966, ], demean = FALSE, intercept = TRUE, > prewhite = TRUE) : > model order: 2

For the following model, library(vcd) arth.polr <- polr(Improved ~ Sex + Treatment + Age, data=Arthritis) summary(arth.polr) where Improved is an ordered, 3-level response I'm looking for a *simple* way to test the validity of the proportional odds assumption, typically done via a score test for equal slopes/effects over the predictors. I do find a po.test= option in the repolr package

Hello Prof. Harrell and dear R-help mailing list, I wish to perform a non-parametric repeated measures anova. If what I read online is true, this could be achieved using a mixed Ordinal Regression model (a.k.a: Proportional Odds Model). I found two packages that seems relevant, but couldn't find any vignette on the subject: http://cran.r-project.org/web/packages/repolr/

Hello! I need some help with the linear discriminant analysis in R. I have some plant samples (divided into several groups) on which I measured a few quantitative characteristics. Now, I need to infer some classification rules usable for identifying new samples. I have used the function lda from the MASS library in a usual fashion: lda.1 <- lda(groups~char1+char2+char3, data=xxx) I'd