similar to: Influence measure + lme ?

Hi all, Suppose I have a 4 x 2 matrix A and I want to select the values in second column such that the value in first column equals to k. I gave the colnames as alpha beta, so I was trying to access the info using A$beta[A[,1]==k], however, I was told it's not a data frame, I can get the object by using dollar sign. I tried data.frame(A), but it didn't work. Any input

Hi all, I tried to do a variance components using nlme, but I got the following warning mesage ##################################################### not meaningful for factors in: Ops.factor(y[revOrder], Fitted) ###################################################### Can someone point out what is the meaning of this warning message? I tried to look at Ops.factor, but I don't

Hi all, If I want to write a for loop to loop through a set of factors, which are not coded in 1,2,3, for e.g in character, possible to write the for loop without changing it to categorical variables? I saw the manual mentions for loop will take a list, but I'm not sure how to create a list here. Any input will be appreciated. Thanks. Yen Lin [[alternative HTML version deleted]]

Hi all, I wonder how could I convert a matrix A to a dataframe such that whenever I'm running a linear model such lme, I can use A$x1? I tried data.frame(A), it didn't work. Should I initialize A not as a matrix? Thanks. Yen Lin [[alternative HTML version deleted]]

Hi all, I just wonder what singular convergence means. Thanks. Yen Lin Error in lme.formula(Data ~ 1, random = ~1 | Wafer/fie/loc, subset = Wafer == : singular convergence (7) [[alternative HTML version deleted]]

Hi all, I have read a data matrix of 304 x 404 using read.table. When I am trying to name the colnames, with first try colnames(L5)<-list(2:305), I keep getting error message such as Error in "colnames<-"(`*tmp*`, value = list(c(2, 3, 4, 5, 6, 7, 8, 9, : length of 'dimnames' [2] not equal to array extent and I don't know why. But, if I look at

influence.measures gives several measures of influence for each observation (Cook's Distance, etc) and actually flags observations that it determines are influential by any of the measures. Looks good! But how does it discriminate between the influential and non- influential observations by each of the measures? Like does it do a Bonferroni-corrected t on the residuals identified by

Hello: I am trying to understand the method 'hatvalues(...)', which returns something similar to the diagonals of the plain vanilla hat matrix [X(X'X)^(-1)X'], but not quite.  A Fortran programmer I am not, but tracing through the code it looks like perhaps some sort of correction based on the notion of 'leave-one-out' variance is being applied. Whatever the

Hello everyone, It seems that I was not successfully attached the code. Here is the code. I appreciate any help! Best, Yen ?? b88207001 at ntu.edu.tw: > Hello everyone, > > I have been working on a code which simply repeatedly appends a > number into a vector and write a file. However, it could not be > properly implemented when I use it. It works when I run it line by

All support on this list is voluntary, and support for old versions of R is not even necessarily on-topic here which is why you keep getting nudged to upgrade. Your "need" for support for an old version is definitely not "our" problem, so I suggest you start looking for a consultant if this issue is that important to you. Such is the nature of volunteer-developed open source

Thanks for the help. I have a reason to continue with R-3.0.3. I used maxLik to estimate econometric models and some of them are better handled with R-3.0.3 (but not later)----a sad reality I do not like. Here is what I did. I downloaded https://cran-archive.r-project.org/bin/windows/contrib/3.0/aod_1.3.zip and installed the zip file, which worked in both RStudio and R (without RStudio). In

Don't choose a mirror. That will override the repos choice. Do update R to a current version if you aren't able to debug this yourself. Duncan Murdoch On 08/10/2020 12:38 p.m., Steven Yen wrote: > Sorry Gentlemen and all. Now this is becoming a joke (to me). I repeated > what I did earlier, with and without the option to set repos suggested > by Duncan. Now it does not work.

Hello everyone, I have been working on a code which simply repeatedly appends a number into a vector and write a file. However, it could not be properly implemented when I use it. It works when I run it line by line. I wonder what is the problem and I appreciate anyone who is willing to help. The function and the example code is attached. Any advice is appreciated! Best, Yen

Barrett & Ling, JASA, 1992, v.87(417), pp184-191 define general classes of influence measures for multivariate regression models, including analogs of Cook's D, Andrews & Pregibon COVRATIO, etc. As in univariate response models, these are based on leverage and residuals based on omitting one (or more) observations at a time and refitting, although, in the univariate case, the

Sorry Gentlemen and all. Now this is becoming a joke (to me). I repeated what I did earlier, with and without the option to set repos suggested by Duncan. Now it does not work. I wonder whether it is dependent on the mirror I chose, but I do not remember the one I chose earlier when it work. I need your help, gentlemen, as I need to use R-3.0.3 for my task. >

Okay, so it's not an RStudio issue. However, I'd guess setting options(repos = "https://cran-archive.r-project.org") at the start of your session could make everything work. (I'm guessing you currently have it set to "http://cran.rstudio.com", which is the source of the last warning below, probably due to an R bug. But since you're using an obsolete

Hello, everyone, There's a problem about zero in R and I really need your help. I have a vector shown as x=c(0.1819711,0.4811463,0.1935151,0.1433675), The sum of this vector is shown as 1 in R, but when I type 1-sum(x), the value is not zero, but -2.220446e-16. I can accept that this value is quite small and could be seen as zero, but there would be a problem when it's not really

Hi, Just a quick question. I wonder if the package of "car" has been removed from the most current version of R. Or what can I do to renew it since I have to do the analysis with that? Thanks, yen

Hello everyone, I know where is wrong. I forget to specify the parameters in my function. Thank you for anyone who was trying to help me! Best, Yen ?? b88207001 at ntu.edu.tw: > Hello everyone, > > It seems that I was not successfully attached the code. Here is the > code. I appreciate any help! > > Best, > Yen > > ?? b88207001 at ntu.edu.tw: > >> Hello

Thanks. You gentlemen please tell me what this means. In R (outside of RStudio) I ran: install.packages("aod") Received a warning (and installation did not seem to go through). Then I tried install.packages("aod",repos='https://cran-archive.r-project.org') Received a warning but it went on to try