similar to: vectorized approach to cumulative sampling

Dear All, I am hoping to implement a barchart using trellis graphics where all elements have a line width of 2. Using trellis.par.set(), I am able to make most elements lwd=2, but not all. In particular, the top of the box (above the upper most strip) and the left y-axis remain one point. Code with a barchart() example is below. Any advice would be greatly appreciated. Thanks for your

Dear R-Helpers, Can anyone please help me to interpret warning messages from zeroinfl (package pscl) while fitting a zero inflated negative binomial model? The console reports convergence and the parameters seam reasonable, but these <<Warning messages: 1: algorithm did not converge in: glm.fit(X, Y, family = poisson()) 2: fitted rates numerically 0 occurred in: glm.fit(X, Y, family =

# Hi All, # # I need to solve a somewhat complex equation at many parameter values for # a number of different parameters. # A simplified version of the equation is: 0= (d1/(h1^2))-(h2*(d2^2)) # I'd like to solve it across a parameter space of d1 and d2, holding # h1 and h2 constant. # It seems that uniroot() can do it, but I don't see how to vectorize it. #

Dear All, I need to produce a multi-panel plot where: 1) groups within panels receive distinct symbols; 2) a linear regression is fit to the entire panel, not the individual groups; 3) a common abline is plotted in each panel. Essentially I would like to merge the two plots below, but I can't seem to figure out how to call both lmline and abline, and also retain the grouping symbols.

Hi All, convhulln {geometry} computes the convex hull of a set of points in n- dimensions via .call, returning the hull itself, and also "unavoidably generates a diagnostic report on an Rterm console." See the example below. I need to access the results of the diagnostic report (specifically the computed hull volume) from within R, but I cannot see how. I've looked at

Dear All, I need to write an R function which computes values of Probabilities for the (standard) normal distribution, ¦µ(z) for z > 0 by summing this power series. (We should accumulate terms until the sum doesn't change). I also need to make the function work for vector of values z. The initial fomular is ¦µ(z) = ( 1/ sqrt(2¦Ð) )* ¡Ò e^(-t^2/2)dt (¡Ò is from -¡Þ, z)

Hi All, I am writing a simulation that examines the effects of species extinctions on ecological communties by sequentially removing individuals of a given species (sometimes using weighted probabilities) and replacing the lost individuals with species identities randomly sampled from the remaining individuals. Thus I use two dataframes. One contains all the individuals and their species

Dear R team I noticed a bug in the "medpolish" function in the "stats" package. When I run medpolish(x), with "x" the matrix in the attached .rda file, I get this rather cryptic error: "Error in if (converged) break : missing value where TRUE/FALSE needed In addition: Warning message: In sum(abs(z), na.rm = na.rm) : integer overflow - use

Hi All, I need to place lists or vectors within dataframes as single elements. However when I try this: df=data.frame(y=1, x=I(list(c("a","b"), c("f","c"), c("a")))) df df[1,'x']=I(c("a","d")) I get this error, even though I am using I(): Error in `[<-.data.frame`(`*tmp*`, 1, "x", value =

Dear R users, About a year ago Deepayan offered a suggestion to incorporate error bars into a dotplot using the singer data as an example <<http://finzi.psych.upenn.edu/R/Rhelp02a/archive/63875.html>>. When I try to utilize this code with a grouping variable, I get an error stating that the subscripts argument is missing. I have tried to insert them in various ways, but cannot

Full_Name: Frederic Schutz Version: 1.6.1 OS: Linux Submission from: (NULL) (128.250.252.193) The subject says everything... This has already been proposed in June 2001 (cf http://maths.newcastle.edu.au/~rking/R/help/01a/2333.html), but no reply followed. I attach the patch that I used on my system to add this feature (including an updated help page), if it can help. Sorry, the two

Hi, Infinite recursion in S3 methods seem to crash R on Windows 2000 (R terminating with the ("Rgui.exe has generated errors...") message, rather than throwing an error. This happens with both Rgui and Rterm. The following toy example triggers this: myf <- function(x, ...) UseMethod("myf") myf.default <- function(x, ...) myf(x) myf(1) ...R crashes... Which I

Obligatory: When will Centos 5.4 be ready? /me goes down into reinforced concrete bunker and locks all hatches and doors.

Hi, My company wants to put a SIP address on their website. The idea is that potential customers can call that address and will be forwarded to our main switchboard. It's fairly easy in theory because my asterisk server has a real IP address, so any calls to sip:<number>@asterisk-server.mycompany.com should connect just fine (except currently it will be blocked by the firewall). Our

I have a new device that takes measurements anywhere from every second, to every 15 minutes (depending on changes). The matrix has a date, time and Y column (Y is the measurement). For three days it is 25,000 rows. How do I average the measurements by every 30 minutes so my matrix is 48 rows per day? I have been working on this and cannot figure out a simple method. Any ideas? Thank you. ----- In

Dear list members, I am having some problems with a LDAP passdb authentication on Dovecot. Before I forget, the specs: it's a Ubuntu 7.10 server running Dovecot 1.0.5 connecting to 2 different machines running LDAP servers: gold with OpenLDAP 2.4.19 and extra with OpenLDAP 2.4.9 (extra is a replication slave of gold). The initial setup of dovecot was with a MySQL passdb, which is still the

I have a dataset of the form Year tosk.fai tosk.isd tosk.gr ....... tosk.total hysa.fai hysa.isd ... and so on. I want to sum all the columns using the first four letters in the columns label(e.g. 'tosk', 'hysa' etc.). How can you do that? Also, the sums should be without the '.total'column (e.g. 'tosk.total') as this serves as a check that

I'm trying to manipulate/change a formula prior to passing it to another function. A simplified example: User passes formula to my function: y~x My function does: lm(transform(y)~x) Here, transform() is added to the model's response. What is the best way to accomplish this?

Hi R! I am reading in a text file which has one column of alpha data and 5 columns of numeric data. There is a header row. I would like the alpha data column to just be character rather than factor. Is there a way to do this, please? I'm thinking that it might be I() but can't figure out exactly how. Thanks, Laura mailto: lauraholt_983 at hotmail.com R 2.0.1 Windows

hello, i wonder how is string represent in R. if i have a string s= "hello", how can i refer to first character in the string s? also if i have s1 = "hello.1", s2 = "ok.1", how can i separate the s1 into "hello" "1" and s2 into "ok" and "1"? I have tried to use the substring function, but i don't where i can get the index