similar to: BradleyTerry "subscript out of bounds"

Dear all, I have carried out a pairwise comparison study that I want to analyze using the BradleyTerry package to establish a rank order of my stimuli. However, BT does not handle ties between stimuli, so I need to find those in my data before I can use that model. The code below goes from the format of my result file(s) to a data frame suitable for BT, but as you can see, there are some ties. I

Hello, I am trying to make a data frame from many elements after running a function which creates many elements, some of which may not end up being real elements due to errors or missing data. For example, I have the following three elements p1s, p2s, and p3s. p9s did not generate the same data as there was an error in the function for some reason. I currently have to delete p9s from the

How do I install and load the BradleyTerry add on package in R 2.5.1 in MSWindowsXP environment? Kalyan Roy Indian Market Research Bureau (IMRB) International New Delhi, India [[alternative HTML version deleted]]

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

I am writing a rails app with my colleague who wrote a lot of the tests. He had to take a leave of absence due to a death in his family, and I am needing helps with changing the model so that the tests will pass in our spec file. Here is the model: class Worker < ActiveRecord::Base attr_accessible :first_name, :last_name has_many :writings end Here is the spec:

I was uncertain if the following is too commercial for the FreeBSD-announce mailing list, but on the other hand, things that promote the use and value of FreeBSD seem like a good thing too. So if this is too commercial for the list you can let me know at brian@mediacity.com. ----------------------------------------------------------------------- Talisman Sequence P1S $299 PCI Video Capture

Hi, I wish to analyze with R the results of a perception experiment in which subjects had to recognize each stimulus among three choices (this was a forced-choice design). The experiment runs under two different conditions and the data is like the following: N1 : count of trials in condition 1 p11, p12, p13: proportions of choices 1, 2, and 3 in condition 1 N2 : count of trials in

Hi, I'm trying the brlr function in a penalized logistic regression function. However, I am not sure why I am encountering errors. I hope to seek your advice here. (output below) Thank you! Your help is truly appreciated. Min-Han #No error here, the glm seems to work fine >

Hello I do not think I have fully grasped how to use Amelia to deal with missing data. For instance, suppose I have a data.frame variable with 4 columns (year, mon, ssn, dev) = (year, month, measurements, standard deviation of the measurement). Of course, there are some random missing values on columns 3 and 4. The measurements are an almost periodic time-series contaminated by noise.

A user logged in with a captial "L" for his username: Jul 7 12:30:31 postamt dovecot: auth(default): cache(Loser,10.47.64.227): miss Jul 7 12:30:31 postamt dovecot: auth(default): shadow(Loser,10.47.64.227): lookup Jul 7 12:30:31 postamt dovecot: auth(default): shadow(Loser,10.47.64.227): unknown user Jul 7 12:30:31 postamt dovecot: auth(default): cache(Loser,10.47.64.227): miss Jul

I am looking for hints on how to estimate ratings for competitors in an ongoing pairwise competition using R... my particular area of interest being the game of Go, but the idea of identifying ratings (on a continuous scale) rather than relative rankings seems easily generalized to other competitions so I thought someone might be studying something related already. I presume the rating of a

I am fairly new at using R/programming in general so I apologize if I am leaving crucial parts of the puzzle out, but here goes. First and foremost this is the error I am receiving: Error in muPriors[priors[, 1:2]] <- priors[, 3] : NAs are not allowed in subscripted assignments This occurs only when I am using observational priors and some number of nominal variables, it does not

Hi, I am trying to maximize a utility function using optim. I have a problem with optim, since some of the values in the caw, mom, me and btm matrixes in the code bellow are missing. Is there a handy way just to skip the missing values in the loop? g <- 5 retp <- NULL object <- function (x) { b1 <- x[1] b2 <- x[2] b3 <- x[3] for(i in 1:nrow(ret)){ for(j

Hello Everyone, ? I need to perform a sample size calculation using a drop the loser/keep the winner design. I've been searching for examples of how to do this using R but haven't found much. ? The most promising possibility thus far is an R?function called DrpLsrNRst. This appears in the book "Adaptive Design Theory and Implementation Using SAS and R". I found the function

Dear R-Help List, I'm trying to implement Firth's (1993) bias correction for log-linear models. Firth (1993) states that such a correction can be implemented by supplementing the data with a function of h_i, the diagonals from the hat matrix, but doesn't provide further details. I can see that for a saturated log-linear model, h_i=1 for all i, hence one just adds 1/2 to each count,

Hi, we try to do a logistic regression with the function glm. But we notice that this function don't give the same results as the SAS proc catmod (differents estimate given). We try to change the contrast on R system with: > options(contrasts=c(unordered="contr.SAS",ordered="contr.poly")) We also try with brlr and logistf functions. Unfortunately, the estimate

Hello all, I would like to plan and analyse a study with "k" treatments (one of which is a "control"), with some binary outcome, in order to find the "best" treatment (e.g: the one with a high number of "successes"). If this was done with a fixed sample size, the analysis is well known. However, I would rather be able to "drop" treatment(s), if

Hello R-sig-debian and (hopefully) Dirk: On Debian wheezy, I have the R packaging that CRAN (you) provide. I run into a little trouble while trying to fiddle with alternative BLAS. I know you and I went around on this last year and I think perhaps I've found something wrong in the framework, or I've just done something wrong. I installed the packages openblas-base and openblas-dev, and

Please always reply-all so the mailing list can record the answer along with your question. I am cc'ing this time. Thanks for the unusual case. -- Sent from my phone. Please excuse my brevity. On March 1, 2018 12:46:53 AM PST, zn l <esetlzn at gmail.com> wrote: >It is my mistake. I find that I input # character in the Chinese. The >R.app won?t display any error messages. >

I don't think there"most". It is normal to have to call the "library" function for each package you need devices to during an R session. "chakri_amateur" <chakri2sai at yahoo.co.in> wrote: > >In my Windows XP, I just have to load package every time I open R, but not >Install !! > >There must be some problem during installation ! >