similar to: Q: extracting data from lm

Hi all, I wondered whether anyone has some advice on a stats-related 'sanity check', as I ran a nonparametric multivariate test (mulrank function as decribed by R. Wilcox, 2005) on both systems, but got different results (please see below for the system-specific outputs)! The functions I used are attached as well. Any advice would be much appreciated! Thanks in advance for getting back to

I'm trying to code more R-like and avoid loops. Here is an example I'm having a hard time getting away from loops with, and the real matrices are rather large, and the computation is too time-consuming. ### Dimensions N <- 2 M <- 3 P <- 4 ### Array and Matrices nu <- array(NA,dim=c(N,M,P)) Lambda <- matrix(1:12,P,M) F <- matrix(1:8,N,P) ### Loop to avoid for (i in 1:N)

I am seeing some curious behaviour with a 1.2.32 system, which I do not understand and so can't work out how to fix it. I have a PRI routed to context default. Here is the complete default context: [default] exten => _9X.,1,Dial(IAX2/m1peer/${EXTEN:1}) exten => _20XX,1,Dial(IAX2/sipeer/${EXTEN}) exten => _X.,1,Dial(IAX2/m1peer/${EXTEN}) exten =>

Dear R-help members, I have a question regarding how to use varComb function to specify a variance function for the "weights" in the gls. I need to fit a linear model with heteroscedasticity. The variance function is exp(c0+nu0*W +nu1*W^2) where W is a covariate. Initially I want to use varFunc to define my own variance function following the instruction in the Pinheiro and Bates

Dear R-help members, I have a question regarding how to use varComb function to specify a variance function for the "weights" in the gls. I need to fit a linear model with heteroscedasticity. The variance function is exp(c0+nu0*W +nu1*W^2) where W is a covariate. Initially I want to use varFunc to define my own variance function following the instruction in the Pinheiro and

Dear R people I want to contribute a new package (various noncentral hypergeometric distribution functions) using C++, but there are many errors when I try to build even a simple package. Here are the most important errors I have encountered so far: 1. Trying to build the DLL, it can't find some of the include files. I had to make the following change in Makefile and copy some .h files:

All, I have some data of animal movements that I'm plotting using xyplot() from lattice. I want to have the date (class POSIXct object) on the Y-axis and the animals longitude on X-axis. Eg. xyplot(date ~ longitude, groups = animal, data = my.data) with data like: animal ptt year month day lon date 125 03P0014 13273 2003 7 10 150.38 2003-07-10 14:03:48 126

On kde, gtk-window-decorator segfault everytime a minimised "Konqueror download window" close it self... Running with --sync fix the problem so no way to have a backtrace :( Cedric

Background: I administer a cluster of RedHat EWS 3U4 Linux workstations at a university. I built R 2.1.0 from source: ./configure \ --prefix=/sscc/opt/R-2.1.0 \ --with-blas=no \ 2>&1 \ | tee NUInstall.configure R is now configured for i686-pc-linux-gnu Source directory: . Installation directory: /sscc/opt/R-2.1.0 C compiler:

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I have a linux system with 4 ethernet interfaces, eth0 goes to the internet, eth1, eth2, and eth3 are NAT''d LANs. I want to use an ingress filter to prioritize bandwidth (downstream from internet) to various IPs. I want to sett it up something like this.... eth0 <--[NAT]--> dummy0 <---> dummy1 <---> eth1,eth2,eth3 dummy1

Dear R users The following function is R code for the main compuations in the article: M. Eliasziw, S Lorraine Young, M Gail Woodbury and Karen Fryday-Field (1994): Statistical Methodology for the Concurrent Assessment of Intrarater and Intrarater Reliability: Using Goniometric Measurements as an Example. Physical Therapy 74 (8); 777-788 The function gives the estimated inter- and intrarater

Dear Rusers, I have met a difficult problem on explaining the differences of principal component analysis(PCA) between R,S-PLUS and SAS/STATA/SPSS, which wasn't met before. Althought they have got the same eigenvalues, their coeffiecients were different. First, I list my results from R,S-PLUS and SAS/STATA/SPSS, and then show the original dataset, hoping sb. to try and explain it.

The behavior of the 'add user script' smb.conf option is rather weird: It is documented as an option to the login parts of the protocol, and used to add users dynamically during the logon process, if they don't exist locally. However, it is also used in the SAMR code when an admin explicitly creates a user. This is actually the more natural use for the parameter, but it is

Dear List, I am trying to estimate a 3 dimensional density through the DPpackage. For example # model sigma <- matrix(c(0.1,0.05,0.05,0.05,0.1,0.05,0.05,0.05,0.1), ncol=3) rnormm<- rmvnorm(n=100, mean=c(5,100,150), sigma=sigma) sigma2 <- matrix(c(10,0.05,0.05,0.05,10,0.05,0.05,0.05,10), ncol=3) rnormm2<- rmvnorm(n=100, mean=c(20,1,110), sigma=sigma) rnormm<-rbind(rnormm,rnormm2)

Hi, we are trying to install Junghann's quadBRI into Dell PowerEdge 2800 system without success. I don't know if the issue can be that Junghann's card fits 32-bit slot and Dell PE 2800 has only 3 PCI-X 64-bit slots. Can this be an issue? We get "CRC errors for HDLC frame" when the card is initialized. Any idea what can be wrong? 1/ We use latest bristuff packages. 2/

When you say that distribution shouldn't be used, do you mean within instcombine rather than some other pass? Or not all as an IR optimization? A dedicated optimization pass that looks for and makes factoring/distribution folds to eliminate instructions seems like it would solve the problems that I'm seeing. Ie, I'm leaning towards the proposal here: https://reviews.llvm.org/D41574

Dear R-users, I have to deal with (txt) files of the form 10092007 24.62 24.31 24.90 11092007 19.20 23.17 22.10 13092007 24.71 27.33 23.10 14092007 27.33 27.90 24.10 15092007 28.22 28.55 24.30 16092007 28.53 29.24 27.40 17092007 24.19 30.64 26.80 18092007 22.60 20.62 28.40 19092007 18.89 21.70 14.70 20092007 21.27 22.92 17.30 21092007 22.38 24.72 18.80 22092007 23.94 24.73 20.40 23092007 25.33

Full_Name: bill langdon Version: 2.4.1 OS: ubuntu Submission from: (NULL) (155.245.58.159) #WBL 22 Feb 2007 ubuntu R.version #platform i486-pc-linux-gnu #arch i486 #os linux-gnu #system i486, linux-gnu #status #major 2 #minor 4.1 #year 2006 #month 12 #day 18 #svn rev 40228 #language R

Dear R users, I have a txt file of the form 10092007 24.62 24.31 24.90 11092007 19.20 23.17 22.10 13092007 24.71 27.33 23.10 14092007 27.33 27.90 24.10 15092007 28.22 28.55 24.30 16092007 28.53 29.24 27.40 17092007 24.19 30.64 26.80 18092007 22.60 20.62 28.40 19092007 8.89 1.70 14.70 20092007 21.27 2.92 17.30 21092007 22.38 24.72 8.80 22092007 23.94 24.73 20.40 23092007 25.33 25.12 22.60

Dear List, Apologies for this off-topic post but it is R-related in the sense that I am trying to understand what R is telling me with the data to hand. ROC curves have recently been used to determine a dissimilarity threshold for identifying whether two samples are from the same "type" or not. Given the bashing that ROC curves get whenever anyone asks about them on this list (and