similar to: Equivalent AIC values

Hi all, I'm working to apply a function that will generate a matrix of results only when a specific criteria is met. I want my final results to be a matrix with both the values that meet the criteria (the results of the function), and those that to do in the same positions in the matrix (the original numbers). Here's a sample of what I would like to do: t.mean.1.c <- c(-15, -20,

Hello, I'm trying to develop a box plot of time series data to look at the range in the data values over the entire period of record. My data initially starts out as a list of hourly data, and then I've been using this code to make this data into the final ts array. # Read in the station list stn.list <- read.csv("/home/kbennett/fews/stnlist3", as.is=T, header=F) # Read in

Hi, I'm installing rgdal but I keep having failures because I have not been able to find a good source of information for the correct configuration settings when installing GDAL. My error from the R install.packages("rgdal") is below. Can someone point me to a good source to tell me how to set after ./configure when installing GDAL? I'd like to be able to work with raster

Dear friends, As we all know, the usual model selection criteria(e.g.deviance,AIC...) in GLMs isn't very good for selecting the best model when overdispersion exist, so we need to adjust the corresponding statistic,see(Fitzmaurice,G.M. (1997) Model selection with overdispersed

Hello, I am wondering if someone can help me. I have the following function that I derived using nls() SSlogis. I would like to find its derivative. I thought I had done this using deriv(), but for some reason this isn't working out for me. Here is the function: asym <- 84.951 xmid <- 66.90742 scal <- -6.3 x.seq <- seq(1, 153,, 153) nls.fn <- asym/((1+exp((xmid-x.seq)/scal)))

Hello I'm analyzing a dichotomous dependent variable (dv) with more than 100 measurements (within-subjects variable: hours24) per subject and more than 100 subjects. The high number of measurements allows me to model more complex temporal trends. I would like to compare different models using GLM, GLMM, GAM and GAMM, basically do demonstrate the added value of GAMs/GAMMs relative to

Could any of you provide me information on a good Polycom phone dealers to utilize. One who provides firmwares ..etc Thank you! Kenny ______________________________________________________ Click here to donate to the Hurricane Katrina relief effort. http://store.yahoo.com/redcross-donate3/

Hi R-help, I am working with US COOP network station data and the files are concatenated in single rows for all years, but I need to pull these apart into rows for each day. To do this, I need to extract part of each row such as station id, year, mo, and repeat this against other variables in the row (days). My problem is that there are repeated values for each day, and the files are fixed width

Well, I'm a recent convert from WBEL. My biggest concern with CentOS is that the community here seems to want to be more than a recompile of RHEL. But WBEL is floundering, what with Katrina and Rita, and there really being only 1 developer behind it, etc. I offer an automated shell script to switch from WBEL4 to CentOS4 (easy, it's hosted on my home DSL line!) It assumes that

Dear mailinglist members, I was wondering if there was a way to re-evaluate the instances of a tree (in the forest) again after I have manually changed a splitpoint (or split variable) of a decision node. Here's an illustration: library("randomForest") forest.rf <- randomForest(formula = Species ~ ., data = iris, do.trace = TRUE, ntree = 3, mtry = 2, norm.votes = FALSE) # I am

Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na + *nb - 1; for(i = 0; i < nab; i++) ab[i] = 0.0; for(i = 0; i < *na; i++) for(j = 0; j <

I'm very new to R and modeling but need some help with visualization of glms. I'd like to make a graph of my glms to visualize the different effects of different parameters. I've got a binary response variable (bird sightings) and use binomial glms. The 'main' response variable is a measure of distance to a track and the parameters I'm testing for are vegetation

Dear Catarina, I prefer to leave the AIC value as NA for the tweedie GLM family because it takes extra time to compute and is only occasionally wanted. It's easy to compute the AIC yourself using the dtweedie() function of the tweedie package. Best wishes Gordon At 03:05 AM 14/08/2007, Catarina Miranda wrote: >Dear Gordon; > >I have also sent this email to R help mailing list,

Hi, I installed Wine07252005, my printer is EPSON LX 300X, it's CUPS, the name is EPSON9Pin. I tested it with a simple Delphi program and have two problems: 1) If i print using a report builder like RBuilder, it prints quite well, but the font does look good in the print out. 2) If I use spooler to print, like following: if StartDocPrinter(PrnHandle, 1, @DocInfo) = 0 then begin

I apologize for the double post. I am curious as to what the usefullness is of the multiple line appearance feature on Polycom phones. I setup our phones to register one line per extension but I hear the IP501's can do three line appearances. Why and how could this feature be applied? Thanks again all. Kenny ______________________________________________________ Click here to

I have a lot (more than one hundred) of files with tables of the same kind (quadratic, same size) and I want to obtain some statistics for every position on them. Therefore, as a result I want another table. I import every table, and create an object “read.table” for it, then I have try to create a “list” or a “data frame“ and directly utilize some functions like sd( ) without success, because it

hi how does one control traffic of ftp using tc.I tried few things ,dont know how to do it.Could you please guide me in it. Thanks in advance Nampreet _________________________________________________________________ All that you wanted to know about Ms Beautiful Lips http://server1.msn.co.in/Profile/katrina.asp

Hi all I'm doing some things with a colleague comparing different sorts of models. My colleague has fitted a number of glms in Genstat (which I have never used), while the glm I have been using is only available for R. He has a spreadsheet of fitted means from each of his models obtained from using the Genstat "predict" function. For example, suppose we fit the model of the type

I'm very new to R and modeling but need some help with visualization of glms. I'd like to make a graph of my glms to visualize the different effects of different parameters. I've got a binary response variable (bird sightings) and use binomial glms. The 'main' response variable is a measure of distance to a track and the parameters I'm testing for are vegetation parameters

I have a zooreg object and I want to be able to generate a value for seasons and 11-day composites paste it onto my zoo data frame, along with year, month and days. Right now I have the following to work from: eg. dat.zoo.mdy <- with(month.day.year(time(dat.zoo)), cbind(dat.zoo, year, month, day, quarter = (month - 1) %/% 3 + 1, dow = as.numeric(format(time(dat.zoo), "%w")))) For