similar to: column permutation of sparse matrix

I came across a problem when one of our clients was not able to access any of the servers on our network. This person has never connected to us before and now for this first time was trying to do it from his home is Houston, TX using earthlink cable service provided by Time Warner. All this information, I think, is important because when I started examining my shorewall logs I found out

Dear all, I often import data sets from other programs that contain user defined missing values. As an example consider a variable MXPLOCO with a user defined missing value of as -9. > summary(MXPLOCO) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's -9.000 4.750 5.000 5.349 6.000 8.500 2.000 The following command successfully converts the value -9 to NA. However,

Dear useRs, I am trying to model a time series with a transfer function. I think it can be put into the ARMA framework, and estimated with the 'arima' function (and others have made similar comments on this list). I have tried to do that, but the results have so far been disappointing. Maybe I am trying to make 'arima' do something it can't... The data are time series of

With rw2010dev I get a strange protect(): protection stack overflow error with a small data frame which otherwise is usable: If anybody wants to have a look I can provide an RData file with the problematic data frame. Doesn't seem to be necessary, the following simulated example generates the error: > testmat <- matrix(1:80, 20,4) > dim(testmat) [1] 20 4 > str(testmat) int

hallo, i have a problem reading in the above mentioned kind of a file. does anybody know an easy way how to read it in? i can save it as a text file that looks like: Datum des Fixings Restlaufzeit 1 Jahr 2 Jahre 3 Jahre 4 Jahre 5 Jahre 6 Jahre 7 Jahre 8 Jahre 9 Jahre 10 Jahre 12 Jahre 15 Jahre 13.01.06 2.819 2.983 3.073 3.137 3.194 3.247 3.302 3.355 3.409 3.459 3.548 3.649 12.01.06 2.847

I am having trouble using the "resize" function (in the package EBImage) with matrices containing double values. Let's say I have a matrix 'b' with these values:   > a = c(1.25,1.555,1.245,1.265,1.656,1.235,1.848,1.959,1.548,1.245454,1.989,1.24) > b = matrix(a,4,3) > b       [,1]  [,2]     [,3] [1,] 1.250 1.656 1.548000 [2,] 1.555 1.235 1.245454 [3,] 1.245 1.848

Hi everyone, I am fitting a simple linear regression model in R. My command is j=lm( Y ~ Sex + begsal + time + int) Call: lm(formula = Y ~ Sex + begsal + time + int) Coefficients: (Intercept) Sex begsal time int 191.916 -241.805 3.969 5.003 3.040 Now I wish to access the values of these coefficients for other purposes

Hi Everyone, I just don't know how to extract the information I want from the summary of a linear regression model fitting. For example, I fit the following simple linear regression model: results = lm(y_var ~ x_var) summary(results) gives me: Call: lm(formula = y_var ~ x_var) Residuals: Min 1Q Median 3Q Max -5.9859 -1.5849 0.4574 2.0163 4.6015 Coefficients:

I have just compiled and installed R-patched on a Sparc 20 running Solaris 2.5.1. The compilation and installation went fine except for the message: Bare word found where operator expected at /usr/src/R/R-patched/share/perl/R/Rdconv.pm line 1491, near "$txtout txt_header" (Missing operator before txt_header?) I get similar messages for lines 1390, 1405, 1420, and 1448 in Rdconv.pm.

Hello, my test.R file contains two huge arrays (>3000 entries), from which R needs to calculate the Pearson Correlation, if I look at the file the numbers look correct. if I run R R < test.R --no-save I see things like this: 0.723, 0.838, 1.002, 0.364, 0.357, 0.227, 0.982+ , 0.963, 0.535, 1.214, 1.270, 0.832, 1.033, 0.632, 2.482, 1.239, 0.743, 1.077, 0.962, 1.052, 1.075, 1.427, 1.395,

Hi all, I'm working with a bunch of large graphs, and stumbled across something useful. Probably many of you know this, but I didn't and so others might benefit. Using pch="." speeds up plotting considerably over using symbols. > x <- runif(1000000) > y <- runif(1000000) > system.time(plot(x, y, pch=".")) user system elapsed 1.042 0.030 1.077

I really would like to be able to read about this in a document but I cannot find my way around the documentation properly Given the code... M <- matrix(runif(5*20), nrow=20) colnames(M) <- c('a', 'b', 'c', 'd', 'e') ind <- c(1,2,3,4) dep <- 5 I can then do... l2 <- lm(M[,dep]~M[,ind]) ## Clearly not useful! summary(l2) I am not sure what

Hi, I am analyzing a data set with greater than 1000 independent cases (collected in an unrestricted manner), where each case has 3 variables associated with it: one, a factor variable with 0/1 levels (called XX), another factor variable with 8 levels (X) and a third response variable with two levels (Y: 0/1). I am trying to see if X1 has an effect on the relationship between X2 and the

Full_Name: Henrik Bengtsson Version: 1.3.1 OS: Sun Solaris 8 Submission from: (NULL) (130.235.3.49) I have been trying to install [R] v1.3.1 locally on a Sun Solaris 8 machine and I ran into the (already known) perl problems that express theselves as: Bare word found where operator expected at /usr/matcent/hb/R/R-1.3.1/share/perl/R/Rdconv.pm line 1390, near "$txtout mywrap"

Dear Rusers, I am now using R and SAS to fit the piecewise linear functions, and what surprised me is that they have a great differrent result. See below. #R code--Knots for distance are 16.13 and 24, respectively, and Knots for y are -0.4357 and -0.3202 m.glm<-glm(mark~x+poly(elevation,2)+bs(distance,degree=1,knots=c(16.13,24)) +bs(y,degree=1,knots=c(-0.4357,-0.3202

Hi there - I am using R1.5.1 on WinNT and the latest MASS (Venables and Ripley) library. Running the following code: >minimod<-glm(miniSF~gtbt*f.batch+log(mxjd),data=gtbt,family="binomial") >summary(minimod,cor=F) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.91561 0.32655 2.804 0.005049 ** gtbtgt 0.47171

Dear ReXperts, I have the below text file output. I need to extract the T, QC, QO, QO-QC and WT columns for the data between T = 10 and T=150. Any ideas? Thanks in advance. ======================================================================================== 1 D C ---CAT-- T THETA QC QO QO-QC QC/QO WT FSD 8 1 0 1.0000E+01

#I am trying to understand how R fits models for contrasts in a #simple one-way anova. This is an example, I am not stupid enough to want #to simultaneously apply all of these contrasts to real data. With a few #exceptions, the tests that I would compute by hand (or by other software) #will give the same t or F statistics. It is the contrast estimates that R produces #that I can't seem to

Here is the result of running 20 threads on 20 physical cores (40 virtual cores). 19002.081139 task-clock (msec) # 2.147 CPUs utilized ( +- 2.88% ) 23,006 context-switches # 0.001 M/sec ( +- 2.24% ) 1,491 cpu-migrations # 0.078 K/sec ( +- 22.50% ) 2,607,076 page-faults # 0.137 M/sec

Summary: Any user can gain root privileges on a Intel Linux system with suidperl 5.003 (having the suid bit, of course) even if "SUIDBUF" and "two suidperl security patches" have been applied. Non-Intel / non-Linux platforms may be affected as well. Quick fix: chmod u-s /usr/bin/sperl5.003 (what else?) Details: There is a nasty bug in mess() (util.c): it is possible to