similar to: Reading data from a dataframe

Thanks for the replies to my prior question. My problem is that R always says object not found when I enter a variable name into a command. I converted a Stata file into an Rdata file by first loading the foreign package by entering require(foreign) Then I asked R to read the Stata file by entering pol572a1<- read.dta("C:\\alex\\Graduate Coursework\\Pol 572\\pol572a1.dta") So

I have recently started learning R. I converted a Stata datafile into an R image file, but I am unable to do anything with the data. For example, I am unable to calculate the mean of any variable or create a scatterplot of some of the data. I always get a message saying, "Error in plot(rgnpc, incmean) : object "rgnpc" not found" when I attempt to make a simple plot of the

Dear r-devel list members, On a couple of occasions I've encountered the issue illustrated by the following examples: --------- snip ----------- > mod.1 <- lm(Employed ~ GNP.deflator + GNP + Unemployed + + Armed.Forces + Population + Year, data=longley) > mod.2 <- update(mod.1, . ~ . - Year + Year) > all.equal(mod.1, mod.2) [1] TRUE > > f <-

Dear Terry, Even the behaviour of lm() and glm() isn't entirely consistent. In both cases, singularity results in NA coefficients by default, and these are reported in the model summary and coefficient vector, but not in the coefficient covariance matrix: ---------------- > mod.lm <- lm(Employed ~ GNP + Population + I(GNP + Population), + data=longley) >

[This question is hopefully straight-forward, but difficult to provide reproducible code.] I'm doing a multivariate bootstrap, using boot::boot(), where the output of the basic computation is a k x p matrix of coefficients, representing a tuning constant x variable, as shown in the $t0 component from my run, giving a 3 x 6 matrix > lboot$t0 GNP Unemployed Armed.Forces

I want to apply Nyblom-Hansen test with the strucchange package, but I don't know how is the correct way and what is the difference between the following two approaches (leeding to different results): data("longley") # 1. Approach: sctest(Employed ~ Year + GNP.deflator + GNP + Armed.Forces, data = longley, type = "Nyblom-Hansen") #results in: # Score-based CUSUM

Hello dear R help group, I have the following setup to analyse: We have about 150 subjects, and for each subject we performed a pair of tests (under different conditions) 18 times. The 18 different conditions of the test are complementary, in such a way so that if we where to average over the tests (for each subject), we would get no correlation between the tests (between subjects). What we wish

>>>>> Fox, John <jfox at mcmaster.ca> >>>>> on Thu, 14 Sep 2017 13:46:44 +0000 writes: > Dear Martin, I made three points which likely got lost > because of the way I presented them: > (1) Singularity is an unusual situation and should be made > more prominent. It typically reflects a problem with the > data or the

>>>>> Martin Maechler <maechler at stat.math.ethz.ch> >>>>> on Thu, 14 Sep 2017 10:13:02 +0200 writes: >>>>> Fox, John <jfox at mcmaster.ca> >>>>> on Wed, 13 Sep 2017 22:45:07 +0000 writes: >> Dear Terry, >> Even the behaviour of lm() and glm() isn't entirely consistent. In both cases,

Dear all, I would like to suggest changing the pairs.formula command such that a command like pairs(GNP ~ . - Year - GNP.deflator, longley) would behave in a similar fashion as lm(GNP ~ . - Year - GNP.deflator, longley) i.e., make a pairwise scatterplot of GNP and all other variables in the (longley) dataframe except for Year and GNP.deflator. The above command, with the

Hello, I have posted this mail a few days ago but I did it wrong, I hope is right now: I have the following doubts related with lm.ridge, from MASS package. To show the problem using the Longley example, I have the following doubts: First: I think coefficients from lm(Employed~.,data=longley) should be equal coefficients from lm.ridge(Employed~.,data=longley, lambda=0) why it does not happen?

Dear all I'm familiarising myself with Ridge Regressions in R and the following is bugging me: How does one get p-values for the coefficients obtained from MASS::lm.ridge() output (for a given lambda)? Consider the example below (adapted from PRA [1]): > require(MASS) > data(longley) > gr <- lm.ridge(Employed ~ .,longley,lambda = seq(0,0.1,0.001)) > plot(gr) > select(gr)

I can plot a time series with: gnp <- ts(cumsum(1+round(rnorm(100), 2)), start=c(1954,7), frequency=12) plot(gnp) But I want to plot more time series in the same plot. How can I do it? =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Changing library(RODBC,boot) to library(RODBC) library(boot) seems to have solved the problem. _____ _____________________________________________ From: Gary Smith [mailto:gary.smith28 at comcast.net] Sent: Saturday, January 31, 2009 12:55 PM To: 'r-help at R-project.org' Subject: can't get package boot to load Hi, I am new to R and I'm totally confused by this problem.

Hello!!, for linear models fit I use Gretl, but now I'm starting to use R, I would like to know if is there some function to obtain a extended summary like in Gretl. I will write a example in Gretl Modelo 1: MCO, usando las observaciones 1968-1982 (T = 15) Variable dependiente: Invest Coeficient St error t-ratio p-value const 377,631 35,0955 10,7601 <0,00001 *** GNP

Hola!! A ver si alguien puede ayudarme!! Para ajuste de modelos lineales normalmente uso Gretl. Ahora estoy empezando a hacerlo en R. Me gustaría saber si existe alguna función que haga un summary extendido como el de Gretl. Os pongo un ejemplo del summary de Gretl. Modelo 1: MCO, usando las observaciones 1968-1982 (T = 15) Variable dependiente: Invest              Coeficiente   Desv. Típica

I'm trying to make a set of contour plots of bivariate kernel density estimates, showing three such plots overlaid, similar to this plot http://euclid.psych.yorku.ca/SCS/Private/Test/ridge-boot2.pdf except that I would like to have the contours *filled* (using transparent colors). To make this reproducible, I've saved the results of KernSmooth::bkde2D() in the following file:

We recently observed a problem with Solaris HVM domains. The bug was seen was seen with a higher number of VCPUs (3 or more), and always had the same pattern: some memory was allocated in the guest, but the first reference caused it to crash with a fatal pagefault. However, on inspection of the page tables, the guests'' view of the pagetables was consistent: the page was present.

Hi, I'm trying to upload two CSV files into R, and I'm having some trouble. I've used the instructions I've found on the Web, but they haven't helped so far. I'm worried this could be because I'm using a Mac. The two data sets are here: http://data.un.org/Data.aspx?d=WDI&f=Indicator_Code%3aNY.GNP.PCAP.PP.CD http://data.un.org/Data.aspx?d=UNODC&f=tableCode%3a1

>>>>> Fox, John <jfox at mcmaster.ca> >>>>> on Wed, 13 Sep 2017 22:45:07 +0000 writes: > Dear Terry, > Even the behaviour of lm() and glm() isn't entirely consistent. In both cases, singularity results in NA coefficients by default, and these are reported in the model summary and coefficient vector, but not in the coefficient covariance