similar to: optim

Listmates, I have had a few boxes where the 'matchname' message is filling my syslog with one new message ever ~40 seconds or so. No performance problems, just the log annoyance. An example of the messages are: Jan 11 13:43:03 nemesis smbd[2246]: [2010/01/11 13:43:03, 0] lib/util_sock.c:1685(get_peer_name) Jan 11 13:43:03 nemesis smbd[2246]: Matchname failed on Adrienne-PC.rbpllc.com

To all, I'm working with code below to produce a map with station data plotted in points, but right now I'm having trouble with the mapping portion of this code states <- data.frame(map("state", plot=FALSE,xlim= c(-85,-75),ylim=c(33,37))[c("x","y")]) usamap<- ggplot(states)+geom_path(aes(x,y)) usamap When I plot this the problem is that the bounds of

I am making a barplot using barplot2 from gplots where each bar represents a specific tree species. I have formatted the species names on the x-axis so that the genus name is above the species name and have then rotated the labels 45 degrees to save room. This is my code: >columncolor<-c("grey20", "grey20", "grey20", "grey70",

I am trying to tweak how my categorical x-axis labels are formatted in my bar graph. Specifically, I would like to a) decrease the spacing between lines (e.g. spacing between Dialium and guianensis) b) right justify the text and c) have each species name align with the center of the corresponding bar graph. A simplified version of my code thus far is as follows. I am using RStudio v.

Greetings! I'm having an issue with the windrose produced by the windrose function from the circular package. For our weather stations in North Carolina I'm helping with a script which takes hourly wind speed and direction data to create windroses for our end users. One of the stations in the mountains frequently reaches wind speed of 40 to 60 mph and in storms can reach wind speed over

I am making a barplot using barplot2 from gplots and would like to format the names of my categorical variables (tree species) on the x- axis so that the genus name is above the species name (to save room). My code so far is: speciesnames<-c("Dialium guianensis", "Inga alba", "Tachigali versicolor", "Brosimum utile", "Caryocar

I am having problem loading packages in the newest version of R (2.13.0) on my Mac. I have tried to install various packages (e.g. lawstat, Rcmdr, car) and load them using the Package Installer and Package Manager menu options but I get the follow error: > library(lawstat) Loading required package: mvtnorm Error: package 'mvtnorm' could not be loaded In addition: Warning

Is there code in R to automatically assign letters to different groups that are found to be statistically significant using a post-hoc test following an ANOVA? For example, let's say I found that relative rates of tree growth were statistically significant between tree species using an ANOVA and then I ran a TukeyHSD post-hoc test to determine which species were statistically

Hi, I have a data frame with columns as follows: tree species (independent variable) and several dependent variables (e.g. carbon, nitrogen, phosphorus). Each row represents one tree sample. Some tree samples are unique species in the data frame while other species were replicated (i.e. rows 1,2,3 may be identical for the "tree species" column but have different values for the

Hi, There's probably some simple way of doing this, but I'm just not seeing it - How do I get the background to be white instead of grey when I have a cloud plot (using the lattices package)? par(bg="white") isn't working. I'm assuming par commands won't work on lattice plots. What should I use instead? Thanks, Adrienne [[alternative HTML version deleted]]

R-listers, I have noticed several posts on issues with difftime producing NA's but they have been for older versions of R. Here's the issue associated with difftime that I am dealing with in R 2.12.2. > preciptime = strptime("01/10/2007 14:00",format="%m/%d/%Y %H:%M") > class(preciptime) [1] "POSIXlt" "POSIXt" > # Now using difftime, this

Hi, If I have two sets of directional data (in radians) and want to compare them with a multifactorial anova. Is it even legitimate to compare circular data with an anova? The books I've picked up from the library don't really say, but it looks unlikely. If it is allowable, is my having stored the data as circular (X = as.circular(A)) something the aov() function will take into account,

Greetings all, Ran into a strange problem with the krige function from geoR. The problem that I am having is that while the krige function seems to work well, the resulting predicted values are all NAs. Given the size of the datasets I am working with can't attach it, but I can provide snippets of the datasets. > casedata station year month day obs mpe bias type

All, I'm working with ggplot2 to produce some graphics that will be part of my master's thesis. Here's what I would like to do here, (attached is an example of how far I've gotten). I'm working on mapping the values of several different parameters for each station in my study. So far I've been working with scale_color_gradientn to handle the coloring. What I would like

R-listers, I'm working on a project where I need to get the lowercase of the county variable in a dataset alike to the example below (only difference is that the full dataset has all the states and counties in the southeast United States). I keep getting this strange error with the lowerize function (which didn't occur the first few times I use the code below), and oddly enough the error

Hello everybody, I've got a bit of a problem with parLapply that's left me scratching my head today. I've tried this in R 2.11 and the 23 bit Revolution R Enterprise and gotten the same result, OS in question is Windows XP, the package involved is the snow package. I've got a list of 20 rain/no rain (1/0) situations for these two stations i and j, all the items in this list look

Dear Everyone, I am using R 2.6.2 on my Redhat AS4. I installed the RPM offered on the cran website. I would like to use R from my bash console with the following command: /path/to/R --vanilla --slave --args < /data/myscript.R this script invokes a range of functions and tools, among others to process a bunch of raster images in a loop using GDAL. The thing is that releasing this command in

> Sergio Gelato wrote; >> I see that the rest of that function has an "if (problem) goto out;" >> after >> every krb5 library call. Doesn't that also introduce measurable time >> differences? Interesting. > I wrote a test case with expect to measure the time difference for valid and invalid user with the same workaround as said before. It seems to

Dear R users, Since a long time now, I have the following error when I want to load the tcltk library in R 64 bit. Loading Tcl/Tk interface ...Error : .onLoad failed in loadNamespace() for 'tcltk', details: call: inDL(x, as.logical(local), as.logical(now), ...) error: unable to load shared object 'C:/Program Files/R/R-2.13.0/library/tcltk/libs/x64/tcltk.dll': LoadLibrary

Dear Rhelpers I have the problem with initial values, could you please tell me how to solve it? Thank you June > p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2))) Error in maxRoutine(fn = logLik, grad = grad, hess = hess, start = start, : NA in the initial gradient > p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2),method="BFGS")) Error in optim(start, func, gr =