similar to: Fwd: Documenting data sets with many variables

Hello, I am just studying the following example from vignette: strucchange-intro, contineousely ending up in an error. This is the given code: 1. library(strucchange) 2. data(USIncExp) 3. if (!"package:stats" %in% search()) library(ts) 4. USIncExp2 <- window(USIncExp, start = c(1985, 12)) A.Modelling: coint.res <- residuals(lm(expenditure ~ income, data = USIncExp2))

Hi All, I have some data where I am doing fairly simple calculations, nothing more than adding, subtracting, multiplying and dividing. I’m running into a problem when I divide one variable by another and when they’re both 0 I get NaN. I realize that if you divide a non-zero by 0 then you get Inf, which is, of course, correct. But in my case I never get Inf, just NaN because of the structure

I am interested in utilizing this so called "double hurdle" model in my study. We can write the model in the following way: if (z'a + u > 0 & x'b + e > 0) y = x'b + e, else y = 0 In the model, consumption y is the (left-) censored dependent variable. e and u are the normally distributed error terms. z'a is the participation equation and x'b is the

Dear All, I am economist and working on poverty / income inequality. I need descriptive statitics like the ratio of education expentitures between different income quintiles where each household has a different weight. After a bit of google search I found 'Hmisc' and 'quantreg' libraries for weighted quantiles. The problem is that these packages give me only weighted quintiles;

Hi, has anybody an example of an Engel curve analysis in R or does there exist a package to estimate and plot Engel curves from expenditure / income data in R? Thanks a million for your hints, Werner

Hi, i try to plot my graph into different device using x11(), but when i do this comes up: > x11(print(plot(A5e$ECAB,A5e$EXPEND,type='p',main='Per capita expenditure > against economic + ability index without outliners',xlab='ECAB',ylab='EXPEND',xlim=c(0,150),ylim=c(150,400), + col='red',col.axis='blue'))) NULL Warning message:

Hi, Problem: I have the ususal linear regression model y=Xb but I cant observe y (expenditure for a particular good) for some y <c where c is fixed for all observations. But I can observe the corresponding covariates x (incom). I used the survivla5 library (survreg(Surv(y,y<c,type="left")~x) for a censored regression model (thanks again for this tip). Is there a way to estimate

On Thu, 24 Oct 2024, Chris Rapier wrote: > Have people given thought to the private key encryption methods in light of > potential quantum attacks? While the recent paper about breaking 50bit RSA > doesn't pose a threat I've been thinking about future harvest now, decrypt > later attacks against CC20 and AES. Are there post quantum ciphers that can > effectively replace

Hi, I have a data.frame with 100 variables and I have assigned a "label", "units" and "category" attribute to each variable. I would like to reorder the variables in the data.frame by the "category" attributes but can't find a way. For example, the first variable is: > attributes(hh$aez) $levels [1] "coastal" "forest"

Dear all, I have uploaded a new package called micEcon (version 0.1-3) to CRAN (an early version of this package has been already presented at useR! 2004). It contains tools for microeconomic analysis and microeconomic modeling. These are for instance: - tools for demand analysis with the 'Almost Ideal Demand System' (AIDS): e.g. econometric estimation, calculation of price and

Dear all, I have uploaded a new package called micEcon (version 0.1-3) to CRAN (an early version of this package has been already presented at useR! 2004). It contains tools for microeconomic analysis and microeconomic modeling. These are for instance: - tools for demand analysis with the 'Almost Ideal Demand System' (AIDS): e.g. econometric estimation, calculation of price and

I need some guidance from someone who is familiar/has some experience with the survey package. The data that I am using is from the Medical Expenditure Panel Survey (www.meps.ahrq.gov <http://www.meps.ahrq.gov/> ). The STRATA and PSU variables are varstr01 and varpsu01 respectively. When I try to specify them with the svydesign function I get an error message. An excerpt of my session

Hi, I would like to estimate a linear expendire with Systemfit package. (method: "SUR") As someone could show me how to define the equations? Thanks. -- Think before you print ! ********************************************************************** Disclaimer: This e-mail may contain confidential informa...{{dropped:9}}

Dear all, Using the data set and code below, I am interested in modelling how egg temperature (egg.temp) is related to energy expenditure (kjday) and clutch size (treat) in incubating birds using the lme-function. I wish to generate the F-distribution for my model, and have tried to do so using the anova()-function. However, in the resulting anova-table, the parameter kjday has gone from being

How do I go about generating a WEIGHTED mean (and standard error) of a variable (e.g., expenditures) for each level of a categorical variable (e.g., geographic region)? I'm looking for something comparable to PROC MEANS in SAS with both a class and weight statement. Thanks. Marc [[alternative HTML version deleted]]

Hello all, I'm an R novice, recently trying to implement R in my research. Using the data frame below, I want to construct a repeated measures model, with energy expenditure (kjday) as dependent of treatment (code) using mass as a covariate. ind mass kjday code 79 15.8 45.216 3 42 16.5 44.64 3 10 14.85 45.216 3 206 15.75 45.216 3 23 12.15 42.336 3 5 14.6 51.264 3 .... 79 16.9

I apoligize in advance for this newbie question on what I perceive as a mostly advanced level list... I did some searching, but would like some of your expert opinions. I'm building an asterisk server to be used in the home, both to learn, and as proof of concept of applying this solution in a home. To keep costs down, I'm considering one X100P ($25 ebay clone) card to connect the

I'm setting up a small office PBX on asterisk, and I've got to the point where I have to decide on what to do about FAX. The current situation is: PRI (E1) connected to an E100 and sip hard phones on their own network. I'd like to add some very limited fax capability to the system. Basically we need to send a couple of faxes a month, and receive about the same number. My options

http://www.nabble.com/file/p20332968/spe.dat spe.dat Hi, i have a problem in sorting out some command in R and i am really hoping some expert can help me out please! i have the spe.dat file which i upload here, and when u read into R u got something like this: A5 <- read.table('spe.dat' ,header=TRUE) A5 EXPEND ECAB MET GROW YOUNG OLD WEST STATE 1 256 85.5 19.7 6.9 29.6

I've been translating an economic model from Python into R, and I thought members of the list would like to see a presentation I've written about it. I've blogged this at http://www.j-paine.org/blog/2018/03/r-taxben-a-microsimulation-economic-model-in-r.html , and the presentation itself is a slideshow at http://www.j-paine.org/rtaxben/R/reveal/rtaxben.html . The slideshow is written