similar to: Reading SPSS .por files

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

Dear all, Say I have the following dataset: > DF x y z [1] 1 1 1 [2] 2 2 2 [3] 3 3 NA [4] 4 NA 4 [5] NA 5 5 And I want to omit all the rows which have NA, but only in columns X and Y, so that I get: x y z 1 1 1 2 2 2 3 3 NA If I use na.omit(DF), I would delete the row for which z=NA, obtaining thus x y z 1 1 1 2 2 2

I am having difficulty installing the package "lattice". First I tried downloading it from the CRAN site (in the normal way) : the message said it worked but when I typed library I got an error message ("there is no package"). Second I tried installing it from the local zip. I have reproduced the result below. utils:::menuInstallLocal() package 'lattice'

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 --

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

Saben si el estimador VARHAC de Den Haan y Levin (1994) ha sido implementado en R? Solamente encontre rutinas para RATS y GAUSS en la pagina de Den Haan. Saludos, Jose Jose Luis Iparraguirre [[alternative HTML version deleted]]

Hola a todos. Estoy atascado en un sitio que no me he encontrado antes, debe ser fácil y rápido, a ver si alguien me puede ayudar. Puesto fácil, tengo: > aa <- matrix(11:19,3,3) > aa [,1] [,2] [,3] [1,] 11 14 17 [2,] 12 15 18 [3,] 13 16 19 Quiero la fila cuyo primer elemento es 12 y el segundo 15. > aa[aa[,1:2] == (c(12,15)),] Error: (subscript) logical

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

Hi, Imagine I have the following data frame: > a <- c(1,NA,3) > b <- c(2,NA,NA) > c <- data.frame(cbind(a,b)) > c a b 1 1 2 2 NA NA 3 3 NA I want to delete the second row. If I use na.omit, that would also affect the third row. I tried to use a loop and an ifelse clause with is.na to get R identify that row in which all records are missing, as opposed to the first

Hi all, After running kruskal.test I have got results (p<0,005) pointing to reject the hypothesis that the samples were draw from the same population. Howerver when I run the kruskalmc there are no significant differences in any of the multiple comparisons. Is that possible? Some clarification? Thanks, Humber <https://sites.google.com/site/humberandrade> [[alternative HTML version

Hello all, I tried to extract the p-value from adf.test in tseries; however, I got the error message such as > ht=adf.test(list.var$aa) > ht$p-value Error in ht$p - value : non-numeric argument to binary operator > ht Augmented Dickey-Fuller Test data: list.var$aa Dickey-Fuller = -2.3147, Lag order = 4, p-value = 0.4461 alternative hypothesis: stationary > ht$data [1]

Hi: I have tried to find the appropriate clust algorithm for mixed type of data. The suggested way I see is: 1. use daisy to get the dissimilarity matrix 2. use PAM/hclust by providing the dissimilarity matrix, to get the clusters but by following this, when the data set grows bigger say 10,000 rows of data, the dissimilarity matrix will be O(n^2), and out of memory will occur. I am

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

Hello everybody, I have got a data set with about 400 companies. Each company has a score for its enviroment comportment between 0 and 100. These companies belong to about 15 different countries. I have e.g. 70 companies from UK and 5 from Luxembourg,- so the data set is pretty unbalanced and I want to do an ANOVA. Somthing like aov(enviromentscore~country). But the aov function is just for

Hi! I've been looking for a way to calculate std errors of marginal effects when I use interaction terms, but with no success. I pretty much have two cases: continuous variable * continuous variable, and continuous variable * binary variable. In both cases, I know how to calculate the marginal effects, even with simulation. But I still can't figure out how to calculate the std errors of

Hi Just a quick problem that I hope is simple to resolve. I'm doing some work with zero inflated poisson models using the pscl package. I can build models using zeroinfl and get outputs fom them with no problem, but when I try to use the predict.zeroinfl function, I get Error: could not find function "predict.zeroinfl". I was using an older version of R, but still had the same

Hello, I wrote a class like so: > rcfdpsuperclass <- setRefClass( > Class="rcfdpsuperclass", > fields = list( > RcfpdVersion = "character"), > methods = list( > initialize = function(){ > 'Populates fields with defaults and lock as appropriate' > initFields( > RcfpdVersion =

Hello, R for Dummies. How can I exclude the first 1000 values of a vector (length 12000)? More generally all the values up to the ith? Thanks for your help, Dr Olivier Charansonney Cardiologue Centre Hospitalier Sud-Francilien, Corbeil-Essonnes, France [[alternative HTML version deleted]]

Dear all, I have just sent an enquiry but probably I hadn’t expressed myself properly. Could anyone help me with the following? When I run the code on my data I get a boxplot with outliers identified by numbers 200 & 201. However, what I would like is to label these outliers with their corresponding “DATA$num” values of the data frame. In this example, the outliers should be labelled as: