similar to: areaplot

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?

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

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)

can you recommend a good manual for R that starts with a data set and gives demonstrations on what can be done using R? I downloadedR Langauage definition and An introduction to R but haven't found them overly useful. I'd really like to be able to follow some tutorials using a dataset or many datasets. The datasets I have available on R are Data sets in package 'datasets':

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

I searched the forms (i.e., R Search)?and come up with the following suggested link: ?http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-create-rotated-axis-labels_003f I tried to implement what I believe was being implied by that URL and came up with the below: barplot(WorldPhones[1,], ??????? ylim=c(0, 50000), ??????? axes=FALSE, ann=FALSE, ???????

Hola: No consigo que la función sqldf () funcione dentro de una función. Alguien puede echarme una mano. En resumen, el problema es que cuando lo ejecuto fuera de una función no tengo ningún problema: ========================== > # install.packages("sqldf") > library(MASS) > library (sqldf) > data(Aids2, package="MASS") > options(digits=3) > table

Hola Carlos: Gracias por responder! Ayer caí en esta página al buscar el problema con google, pero no conseguir entender y implementar lo que sugiere. Sin entender nada, lo que hice fue: ======================== > library(MASS) > library (sqldf) # Per SQL > data(Aids2, package="MASS") > SQL_PROVA = function (XDADES, XWHE) + { + XDADES_SQL = sqldf ("select * from

For an application of ridge regression, I need to get the covariance matrices of the estimated regression coefficients in addition to the coefficients for all values of the ridge contstant, lambda. I've studied the code in MASS:::lm.ridge, but don't see how to do this because the code is vectorized using one svd calculation. The relevant lines from lm.ridge, using X, Y are:

Full_Name: Edward McNeil Version: 2.6.1 OS: Windows Submission from: (NULL) (203.170.234.5) I want to show students how the survival object looks like in R. Reproducible example: library(MASS) data(Aids2) attach(Aids2) status <- status=="D" stime <- death-diag surv <- Surv(stime, status) D <- data.frame(stime, status, surv) head(D,20) stime status x..i.. 1 176 TRUE

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, experts! I've tried to convert MASS.Rd, for instance: R CMD Rdconv -t txt MASS.Rd -o MASS.txt 2>MASS.log Comp was thinking for a long, about 30s on Linux 600MHz/128M, but after that R didn't produce the result MASS.txt file, but very large log. Here is its header: ******* Syntax error: mismatched or missing brackets in /----- \name{Aids2} \alias{Aids2} \title{

Thank you so much David, here is correction: d1=suppressWarnings(read.csv("/Users/anamaria/Downloads/B1.csv", stringsAsFactors=FALSE, header=TRUE)) d1$X <- NULL d2=as.matrix(sapply(d1, as.numeric)) pdf("~/graph.pdf") b<-barplot(d2, legend= c("SYCL", "CUDA"), beside= TRUE,las=2,cex.axis=0.7,cex.names=0.7,ylim=c(0,80), col=c("#9e9ac8",

Panasonic HC-X800 is expert HD camcorder with 3MOS system Pro, shooting 1080 60p/50p AVCHD footage at bitrates of up to 28Mbit/sec and exhibiting outstanding HD images with minimum noise. The HC-X800 records AVCHD MTS videos on P2 card. You may want to edit AVCHD MTS videos you shoot with iMovie. The problem is that iMovie often cannot import and recognize the 60p/50p footages, nor can it accept

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

MTS is a file extension for an AVCHD video clip format for high-definition video, s supporting 1080i (a high definition video format with 1080 horizontal scan lines, interlaced) and 720p (720 horizontal scan lines, progressive scan, rather than interlaced) in a relatively small file size. Recently more and more people have updated their Mac computer to Mac OS X 10.9 Mavericks, but video

Thanks to Andy and Thomas, Reading help(hclust) more carefully would have done it but sometimes you do not see the wood for the trees... So hc$merge does exactly what I want. I have never been aware of the command str to get the structure of an R-object. It seems pretty useful to me. Thanks, Arne > -----Original Message----- > From: Liaw, Andy [mailto:andy_liaw at merck.com] >

sorry for big posting direct me please on my errors samba 3.0.25b kerberos heimdal 0.8.1-p2 # kinit adminuser adminuser@USR.NW.MTS.RU's Password: # klist Credentials cache: FILE:/tmp/krb5cc_0 Principal: adminuser@USR.NW.MTS.RU Issued Expires Principal Jul 31 11:22:18 Jul 31 21:22:18 krbtgt/USR.NW.MTS.RU@USR.NW.MTS.RU # net -d 3 ads join ads -U adminuser

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 specification of the model. That's not to say that it *never* makes sense to allow singular fits (as in the situations you mentions). I'd favour setting

https://bugzilla.samba.org/show_bug.cgi?id=11609 Bug ID: 11609 Summary: Incorrect (or at least dangerous) behaviour of --append-verify Product: rsync Version: 3.1.1 Hardware: All OS: All Status: NEW Severity: normal Priority: P5 Component: core Assignee: