similar to: (no subject)

R-helpers I have 120 small Excel sheets to read and I am using library(xlsReadWrite): one example below. I had hoped to read sheets by looping over a list of numbers in their name (eg Book1.xls, Book2.xls, etc). I thought I had seen examples which used eval-parse-paste in this way. However, I have not been able to get it to work.. 1. is this a feasible approach? 2. if not

Does anyone know how to read SPSS .por (ie Portable) files? The foreign package only deals with SPSS .sav files and not with those with a "por" extension. Thanks, José Mr José Luis Iparraguirre D'Elia Senior Research Economist Economic Research Institute of Northern Ireland Floral Buildings 2 - 14 East Bridge Street Belfast BT1 3NQ Northern Ireland United Kingdom

I have a data table with 712 cases (rows) describing young people’s activities for 72 months each case has been classified into one of 5 clusters. The first 72 columns are monthly activities coded 1 to 6 (e.g. school =1) and the 73rd column is the cluster number of the case. I wish to summarise the distribution of monthly activities by cluster e.g for cluster 1: 6 months school; 24 further

I have output from a program which produces a distance matrix I want to read into a clustering program in R. The output is a .txt file and is 'almost' lower triangular in the sense that it is just the triangle below the diagonal. So for example a 4-by-4 distance matrix appears as, 1 2 3 4 5 6 i.e. it looks like a lower triangular of a 3-by3. I thought I might be able

Dear all, Does anyone know if the PLM package (to run Panel Data Analysis) accepts quarterly data? The package vignette and documentation only use annual data -and the only time index available seems to work for years. José Mr José Luis Iparraguirre Senior Research Economist Economic Research Institute of Northern Ireland 2 -14 East Bridge Street Belfast BT1 3NQ Northern

Hi R-users,   I have this code here: library(numDeriv)   fprime <- function(z) { alp  <- 2.0165;   rho  <- 0.868;   # simplified expressions   a      <- alp-0.5   c1     <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)   c2     <- sqrt(rho)/(1-rho)   t1     <- exp(-z/(1-rho))   t2     <- (z/(2*c2))^a   bes1   <- besselI(z*c2,a)   t1bes1 <- t1*bes1   c1*t1bes1*t2 }   ## Newton

..::: Drunk Young :::... http://210.101.95.254/drunk-young ^^^^ ^^^^ Click here to Join Drunk Young will take you to the new level of excitement. It is a guarantee of Your pleasure, because we choose models very seriously & carefully, depending on you offers and suggestions. 'Drunk Young' is regularly updates to bring you more and more sweet nymphets 100%

I'm attempting to produce something like a violin plot to display how y changes with x for members of different groups (My specific case is how floral area changes over time for several species of plants). I've looked at panel.violin (in lattice), which makes nice violin plots, but is really set up to work on a single variable - the area trace represents the frequency of each value of x

Does anyone know how to retrieve the residuals after running an ERS unit root test on a time series (urca package)? There are no 'residuals' or 'res' (or similar) slots in the ur.ers-class, but a plot instruction - ie plot(ur.ers(....)) renders amongst other things a residual plot. If an augmented Dickey-Fuller test is run, an object of ur.df-class is generated which contains a

I am using the function prop.test (base package), which returns a list with class "htest". All I want to do is to assign one of its values to a variable, but I do not want R to print the results and added warning message whenever I invoke the function. How can I prevent R from printing on using prop.test? Regards, Jos? Luis Mr Jos? Luis Iparraguirre Senior Research Economist Economic

I could not find a package to run panel data unit root tests in R (there's a STATA module, though - PANELUNIT- and routines for Sarno and Taylor's MADF test, Levin-Lin-Chu's test, Im-Pesaran-Shin's Test, etc, to do this). Hence, I am toying with the idea of having a go at writing up one for R, but would like to know whether anyone is already working on this. José Mr

Hello R users, I would like to reduce the number of for loops in my code. I build these matrices (50000 times). The main diagonal are 1s and the two sides along the main diagonal mirror each other as follows: i<-5 fi<-matrix(0,nrow=i,ncol=i)#floral inhibition matrix for(r in 1:i){ for(c in 1:i){ if(r==c){ fi[r,c]<-1 }else if(r<c){ fi[r,c]<-1-runif(1)^.5 }else{ fi[r,c]<-fi[c,r]

Dear All, It's been several times I observed a very strange behaviour with R 1.2.0 under Windows NT 4.0 (with sp5): accuracies are lost in many cases. Things are back to normal when I exit R and restart it. Here is an example of my last encounter with this: > bt1 <- c(7.3, 8.1, 8.4, 7, 6, 7.2, 6.7, 8.3, 9, 7.4, 8.6, 6.2, 6.9, 7.6, 5.3, 5.6) > bt1 [1] 7 8 8 7 6 7 6 8 9 7 8 6 6 7 5 5

Hi, I am new to this list as a poster, but a reader for some time. I've using R for several weeks now, and I have a lot of questions about certain procedures. Here I go: I want to test if there are differences in the time spent by pollinators visiting flowers of a given plant species, according to a number of experimental manipulations made on those flowers. All experimental

Our ferret 0.10.13 index has been slowly growing on our debian server and has just got up over 14,000 records. Yesterday I randomly noticed that one search I did was suddenly giving whack, unexpected results. I have spent much time trying to track the problem. Tried ferret 0.10.9 - no change. Tried on a windows machine - where it works fine, and doesn''t give weird results (which just

Hi r-users,   I hope somebody can help me to understand the error message.  Here is my code; ## Newton iteration newton_gam <- function(z) { n   <- length(z)   r   <- runif(n)   tol <- 1E-6   cdf <- vector(length=n, mode="numeric")   fprime <- vector(length=n, mode="numeric")   f   <- vector(length=n, mode="numeric")     for (i in 1:1000)   {

Dear R users, Consider the first two columns of a data frame like this: > z[,1:2] x y 1 1 1 2 2 2 3 3 3 4 1 4 Imagine that y represents the times that the value x happens in a population. But z is not exactly a frequency table, because in z we have x=1 twice. So, the x=1 in the first line and the x=1 in the fourth are not the same, differing according to a third variable in

Hi, I have a Data Set x and I want to check with a Kolmogorow-Smirnow-Test, if x comes from a Subbotin Distribution, whose density function is: function(x,location,scale,tail) # Exponential power (=Subbotin) { const<- 2*scale*tail^(1/tail) *gamma(1+1/tail) z<- (x-location)/scale exp(-1/tail*abs(z)^tail)/const } How can i do this? Thank you --

Hallo. Here is a simplification of the automake stuff for nut 2.2.2. The idea is simple: libcommon.a contains stuff for the binaries libparseconf.la contains stuff for binaries and shared library Both are NOINST, automake understands perfectly our intentions. The rest is only a removal of the implicitly defined stuff. On top of the parallel build patch I got build success on all platforms

Hello, all. I've got a problem where we set up call pickup for a customer. If the Bob's extension rings and Bob is in Jim's office, Bob can press the button on his Snom 320 that says "Bob" and pick up his line. It works great for calls coming in from the outside but does not work for internal calls. Internal calls generate a app_directed_pickup.c:204 pickup_exec: No