similar to: Confidence Interval

Hi, so I don't speak computer and I have no idea what this code is telling the program to do, but I apparently need to be able to find and isolate influencial observations. Problem, I have no idea what the error means and where it may be from in the code. error I get is below the code { ## OLS results NameC<- lm(gpanew~female+female:lastinit+agenew+canadian+mom_ed+yearstudy) ## default:

Dear List-members, I found a strange behaviour in the lqs function. Suppose I have the following data: y <- c(7.6, 7.7, 4.3, 5.9, 5.0, 6.5, 8.3, 8.2, 13.2, 12.6, 10.4, 10.8, 13.1, 12.3, 10.4, 10.5, 7.7, 9.5, 12.0, 12.6, 13.6, 14.1, 13.5, 11.5, 12.0, 13.0, 14.1, 15.1) x1 <- c(8.2, 7.6,, 4.6, 4.3, 5.9, 5.0, 6.5, 8.3, 10.1, 13.2, 12.6, 10.4, 10.8, 13.1, 13.3, 10.4, 10.5, 7.7, 10.0, 12.0,

I'm trying to use least trimmed squares using ltsreg with nsamp="exact". When I use the following: rg <- ltsreg(x,y,nsamp="exact") I get: Error in lqs.default(x, y, nsamp = "exact", method = "lts") : NAs in foreign function call (arg 10) In addition: Warning message: NAs introduced by coercion Incidentally, there are no missings in x or y,

I have found that specifying different "sizes" in the sample command has a funny effect on the random sampling. The code below is a condensed version of a function I wrote to simulate a bootstrap method. For simplicity, I eliminated the internal bootstrap loop, but kept a statement to draw one bootstrap sample, because this is where the problem occurs. The output (mean(y)^2) should be

Hello, I've a problem with a self written routine taking a lot of memory (>1.2Gb). Maybe you can suggest some enhancements, I'm pretty sure that my implementation is not optimal ... I'm creating many linear models and store coefficients, anova p-values ... all I need in different lists which are then finally returned in a list (list of lists). The input is a matrix with 84 rows

Hi, I have a complex-valued vector X in C^n. Given another complex-valued vector Y in C^n, I want to find a permutation of Y, say, Y*, that minimizes ||X - Y*||, the distance between X and Y*. Note that this problem can be trivially solved for "Real" vectors, since real numbers possess the ordering property. Complex numbers, however, do not possess this property. Hence the

Hello all, I looking at package dse or vars or mAr I know how to simulate a VAR(p) process, my problem is that most of those processes are unstable (not weakly stationary). Do anybody know how to generate a random VAR (or VARMA even better) process that is weakly stationary? Thanks -- View this message in context: http://r.789695.n4.nabble.com/simulating-stable-VAR-process-tp4261177p4261177.html

Hi, I'm completely new to R, I am all at sea with the interface and the confusing help files, so would appreciate some help to do a simple task. Need to present the mean and variance of 100 different samples of poisson distributions (N=1000, with fixed lambda) in a file in two columnns, and then produce histograms. So far I have figured out: > N <- 1000 > x <- rpois(N, 3.1) ,

Hello all I have a feeling this is very simple......but I am not sure how to do it My boss has two variables, one is an average of 4 numbers, the other is an average of 3 of those numbers i.e var1 = (X1 + X2 + X3 + X4)/4 var2 = (X1 + X2 + X3)/3 all of the X variables are supposed to be measuring similar constructs not surprisingly, these are highly correlated (r = .98), the question is how

Hi, I' am doing a stats project using R to work out the size of a t-test and wilcoxon test depending on the distribution and sample size. I just can't get it to work - I want to put my results from the function size() into an array.At the moment I keep getting the error message:Error in res[distribution, test, samplesize] <- results : subscript out of boundsCan anyone tell me where

Hi, I have been using the lqs function with method='lms'. However the results I get are a little different from the results noted by Rousseeuw & Leroy (Robust Regression and Outlier Detection) and I was wondering how to use these results for outlier detection. I'm using the stackloss dataset, for which the original Rousseeuw et al. program points out that observations 1,2,3,4

hello all i am busy testing some code that will be used for parallel computations on a Moab cluster. I have written a C++ program and am linking it with R using the .C command - this works correctly. Extracts of the r code (as well as the *.sh file) is included below. I am able to run the code successfully if I remove all of the doMPI commands. i.e. outs=.C("prog1", all of the

Hi, I want to calculate the R-squared between two variables. Can you advice me how to identify and remove the outliers before performing R-squared calculation? Thanks, Kan --------------------------------- SBC Yahoo! DSL - Now only $29.95 per month! [[alternate HTML version deleted]]

Hi All I have been struggling with this model for some time now and I just can't get it to work correctly. The messages I get when running the code is: DLSODA- Warning..Internal T (=R1) and H (=R2) are such that in the machine, T + H = T on the next step (H = step size). Solver will continue anyway. In above message, R = [1] 0 0 DINTDY- T (=R1) illegal In above message, R = [1]

Hi there! I appear to be having problems with trusted domains using 3.0a21. Using Samba as a member server, as user "nirs" on domain "Exanet-il" I am unable to log on to a member server ("snow") who is a member of the "Exanet-qa" domain. The "exanet-qa" is supposed to trust "exanet-il", and it all worked with 2.2 . This is the relevant

On Fri, Dec 18, 2015 at 01:14:25AM +0900, HAT wrote: > Hi, > > >> not sure, but the colon in xattr names is probably not handled > >> correctly at this point. Can you please file a bugreport so we can > >> track this? > > > > Bug 11162 > > https://bugzilla.samba.org/show_bug.cgi?id=11162 > > I tested samba 4.3.1 on Fedora 23 and rawhide.

Hi all, I have a huge matrix of TRUE/FALSE table like following, and I want to count the number of TRUEs in each row. Instead of looping through each row and do length(Z[Z==TRUE]), I am wondering if there is an easier way of doing this. [,1] [,2] [,3] [1,]TRUE FALSE FALSE [2,]FALSE TRUE TRUE Thank you in advance. Wendy -- View this message in context:

> not sure, but the colon in xattr names is probably not handled > correctly at this point. Can you please file a bugreport so we can > track this? Bug 11162 https://bugzilla.samba.org/show_bug.cgi?id=11162 -- HAT

Hello! I have created a bplot-figure using this code: *file <- "2dcali_red.ttt" ux<-as.matrix(read.table(file, dec = ",")) mode(ux)<-'numeric' vel<-ux[,1] ang<-ux[,2] x<-ux[,3] y<-ux[,4] dat<- data.frame(ang=ang, x=x,y=y) require(rms) ddist2 <- datadist(dat) options(datadist="ddist2") fitn <- lrm(ang ~ rcs(x,4) +

Hello, I have read over all of the FAQ's I could find, read all of the documnetation that I could find and still I cannot find an answer to the following problem. When I use the smbclient to check to see if Samba is working correctly I get the following response. This is running on Solaris 7 on a Sun box. If you need any other info please let me know. Any help would be very much