search for: iparraguirr

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 Tel: 028 9072 7350 Fax: 028 9031 9003 www.erini.ac.uk <http://www.erini.ac.uk/> [[alternative H...

...row in which all records are missing, as opposed to the first row in which no records are missing or the third one, in which only one record is missing. How can I get R identify the row in which all records are missing? Or, how can I get R delete/omit only this row? Thanks in advance, José José Iparraguirre Chief Economist Age UK T 020 303 31482 E Jose.Iparraguirre@ageuk.org.uk<mailto:Jose.Iparraguirre@ageuk.org.uk> Tavis House, 1- 6 Tavistock Square London, WC1H 9NB www.ageuk.org.uk<http://www.ageuk.org.uk> | ageukblog.org.uk<http://ageukblog.org.uk/> | @AgeUKPA<http://twitter...

...r=list(c(1,0),c(0,0),c(0,0),c(0,0),c(0,0))) d. So, I get estimates for the following coefficients: ma1; CO1-AR1; CO1-MA0; CO2-MA0; CO3-MA0; CO4-MA0; DUM-MA0 Can you spot any problems with the general modelling approach or in this specification? Many thanks in advance, José José Iparraguirre Chief Economist Age UK T 020 303 31482 E Jose.Iparraguirre@ageuk.org.uk<mailto:Jose.Iparraguirre@ageuk.org.uk> Twitter @jose.iparraguirre@ageuk Tavis House, 1- 6 Tavistock Square London, WC1H 9NB www.ageuk.org.uk<http://www.ageuk.org.uk> | ageukblog.org.uk<http://ageukblog.org.uk...

..., but I wouldn't get the corresponding values for z (ie 1 2 NA) Any suggestions about how to obtain the desired results efficiently (the actual dataset has millions of records and almost 50 columns, and I would apply the procedure on 12 of these columns)? Sincerely, Jose Luis Jose Luis Iparraguirre Senior Research Economist Economic Research Institute of Northern Ireland [[alternative HTML version deleted]]

Hi, i would like to plot a few hundred .csv files. Each file contains one curve with x,y values to plot. I have been searching for "gnu r read many files qplot" and similar words. I found for loops that use assign to generate one variable containing a dataframe. When i uesed the classic "plot' command i could add the curves with something like for... {

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

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 Ireland United Kingdom Tel: +44 (0)28 9072 7365 [[alternative HTML version deleted]]

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

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:

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

Hello, Is there a function in R by which one can run AR model with Bayesian Information Criteria (BIC)? To my knowledge, functions ar and ar.ols could select the order only by AIC. Thanks, Miao [[alternative HTML version deleted]]

Hi all, I want to test whether the MEAN of two different variables, (and different number of observations) are the same. I am trying to use the anova test but it doesn't seem to like that the number of observations are different: a=c(1:5) b=c(1:3) aov_test=aov(a~b) >>Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE) : variable lengths differ (found for

Dear colleagues. I wonder if anybody knows about a procedure in R to calculate the Berger-Parker index and the Brillouin Diversity index. I searched the net but did not find anything good about it. Thanks a lot for any help Mayo, Elvira. [[alternative HTML version deleted]]

Hello, Before I start to construct something inferior myself I would like to know if such a plot function is already out there. I have count data of a small survey with a few questions which use a scale from "Strongly disagree" to "Strongly agree" and similar things. What I would like to have for a slide is a plot which shows for each question a line over this scale and the

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

...the income figures for those individuals in poverty, I need to 'un-gross' the data to get back to n=2,263 -and table() simply does not do the trick, because of exactly the same situation in the example above. So, please, how can I retrieve z? Many thanks, Jose Mr José Luis Iparraguirre Senior Research Economist Economic Research Institute of Northern Ireland 2 -14 East Bridge Street Belfast BT1 3NQ Northern Ireland United Kingdom Tel: +44 (0)28 9072 7365 [[alternative HTML version deleted]]

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 =