similar to: bootstrapping lambda values from projections matrices

Dear list! I have encountered an odd problem with R-WinEdt. The program is mostly running fine and I had no problems with the installation. But several times a day I'm typing in R-WinEdt there appears a menu with Chinese characters where the curser is and I'm not able to type normally. I don?t know any other way to stop it except restarting R-WinEdt. After that it works normally for a

Dear R-users, I am trying to make a plot in R where x and y are plotted in a regular way, but the x axis corresponds to another set of values. For example I have x,y and T (all 29 values) x<- c( -1.31846232, -1.04744756, -0.87034853, -0.72883370, -0.60618971, -0.49501845, -0.39128988, -0.29250120, -0.19694055, -0.10334039, -0.01069355, 0.08185470, 0.17507665, 0.26971270, 0.36651292,

Dear R-users, I am trying to make a plot with ggplot-geom_tile(), but cannot remove some unwanted (white) lines through my plot. Below a reproducible example: ##### library(ggplot2) tot=as.data.frame(rep(seq(-50,50,5),each=21)) names(tot)="precip" temp=rep(seq(-5,5,0.5),21) tot$temp=temp disc=array(dim=c(21,21)) for(i in 1:21){ for(y in 1:21){ temp<-

Dear all, I have a dataframe which looks like this (dummy): date<-c("jan", "feb", "mar", "apr", "may", "june", "july", "aug","sep","oct","nov","dec") col1<-c(8.2,5.4,4.3,4.1,3.1,2.5,1.1,4.5,3.2,1.9,7.8,6.5) col2<-c(3.1,2.3,4.7,6.9,7.5,1.1,3.6,8.5,7.5,2.5,4.1,2.3)

Hi, Recently I started working with netcdf IPCC climate model data and I am using R to analyze these data. Some problems occur while doing this, as I am using the output of different climate models which are different in for example time calendar and grid coordinates. To analyze monthly means I have been using the function /chron/ for 365 day calendar time series. Originally the time is in

Hi all, I have used this library to create a (360 day)calendar for my rainfall data (which is divided over 9 gridcells): ## CODE## library(udunits) utInit() calendar <- att.get.nc(nc,'time','calendar') T <- var.get.nc(nc,"time") times.list <- utCalendar(T,"days since 1961-01-01",style='array',calendar=calendar) ## END CODE## **** To separate

Hi all, I have made a plot with ssplot, using a SpatialPointsDataFrame. The content is quite simple, as I have 9 grid points with lon/lat coordinates and 9 values attached to these coordinates. They are in a square area of 3 by 3 gridboxes. I would like to lay a map from maps() over these values, but when I try this, the grids of the maps (smaller) do not overlap with the grids I have

does anyone know how to extract Shrinkage intensity lambda and lambda.var values after run cov.shrink(x)? thanks, KZ [[alternative HTML version deleted]]

Hi, Erin: A cleaner way is to pass "n2" to "outer" as a "..." argument, as in the following modification of your code: boot1 <- function(y,method="f",p=1) { n1 <- length(y) n2 <- n1*p n3 <- n2 - 1 a <- 0.5*(outer(1:n3,1:n3,function(x,y, n2.){n2. - pmax(x,y)}, n2.=n2)) return(a) } y1 <- c( 9, 8, 7, 3, 6) boot1(y=y1,p=4)

Dear R and S People: I have run across something very strange. Here is a function that I wrote for R: boot1 <- function(y,method="f",p=1) { n1 <- length(y) n2 <- n1*p n3 <- n2 - 1 a <- 0.5*(outer(1:n3,1:n3,function(x,y){n2 - pmax(x,y)})) return(a) } and here is the R output: > y1 [1] 9 8 7 3 6 > source("boot1.R") > boot1(y=y1,p=4) [,1] [,2]

Hi all, I am using glampath package for L1 regularized logistic regression. I have read the article " L1 regularization path algorithm for GLM" by park and Hastie (2006). One thing I can't understand that how to find best lambda for my prediction. I want to use that lambda for the prediction not the entire set. thanks. -- View this message in context:

Usually when using raw glmnet I let the implementation choose the lambdas. However when training via caret::train the lambda values are predetermined. Is there any way to have caret defer the lambda choices to caret::train and thus choose the optimal lambda dynamically? -- Yang Zhang http://yz.mit.edu/

Hello, using the lda-method from MASS-package I was wondering whether there is a built-in method for figuring out Wilks' Lambda? Searching the Web I found in the r-help archive a thread form june 2002, but it didn't help me. I understand I can use manova and its summary-method to get Wilks' Lambda on the screen, but I don't see the connection to lda from MASS. Or does

Hello! As I remember there were questions concerning wilks' lambda lately. But I lost those e-mails. So I dare ask again. Is there a way of calculating wilks' lambda for lda? Christoph. -- Christoph Lange clange at epost.de Verhaltensbiologie, FU Berlin 838-55068 Haderslebener Str. 9, 12163 Berlin

Hi, I'm a beginner in the field, I have to perform the ridge regression with lm.ridge for many datasets, and I wanted to do it in an automatic way. In which way I can automatically choose lambda ? As said, right now I'm using lm.ridge MASS function, which I found quite simple and fast, and I've seen that among the returned values there are HKB estimate of the ridge constant and L-W

Dear useRs, I'm pleased to announce that version 1.1.0 of lambda.r is now available on CRAN (http://cran.r-project.org/web/packages/lambda.r/). This package provides a complete functional programming environment within R (and is backwards compatible with S3). Lambda.r introduces many concepts including: . Multipart function definitions . Guard statements to control execution of functions .

Dear useRs, I'm pleased to announce that version 1.1.0 of lambda.r is now available on CRAN (http://cran.r-project.org/web/packages/lambda.r/). This package provides a complete functional programming environment within R (and is backwards compatible with S3). Lambda.r introduces many concepts including: . Multipart function definitions . Guard statements to control execution of functions .

Hi all, Does anyone know of a function that return the best lambda in box-cox transformation? many thanks for any help. -- Cristiane S. Rocha Laboratorio Genoma Funcional - Bioinform?tica Centro de Biologia Molecular e Engenharia Genetica Universidade Estadual de Campinas Campinas - SP - Brasil Tel:(19)3788-1119

Greetings comRades: Has anyone worked out a procedure to calculate the estimated covariance between the beta and lambda estimates in the likelihood-based geostatistical model implemented in package geoR? If so, please help me out by letting me know. Do I have reasons to expect that this covariance term would be big or perhaps more happily, that it could be discarded without much effect?

Hi R People: Sorry about the off topic question. Does anyone know the reference for "-2 Log Lambda is approx dist. Chi square", please? It may be Bartlett, but I'm not sure.... thanks in advance! Sincerely, Laura Holt mailto: holtlaura at gmail.com