similar to: No subject

Dear Peter. I guess you know that Jim Lindseys code include nordr and ordglm in library gnlm - I attach the htmls which do various linear and nonlinear ordinal regressions - exemplified with just the data mentioned, McCullagh (1980) JRSS B42, 109-142. I had it work very fine. -----Oprindelig meddelelse----- Fra: Peter Malewski <p.malewski at tu-bs.de> Til: Troels Ring <tring at

min of an empty vector gives inf, should be NA / NULL ??? > min(c()) [1] 2147483647 Warning message: no finite arguments to min/max; returning extreme. in: min(..., na.rm = na.rm) > min(c(NA,NA),na.rm=T) [1] 2147483647 Warning message: no finite arguments to min/max; returning extreme. in: min(..., na.rm = na.rm) (the same for max) --please do not edit the information below--

Patrick Lindsey has made available a library devoted to ordinal models available at: http://www.luc.ac.be/~plindsey/publications.html Best wishes Troels Ring -----Oprindelig meddelelse----- Fra: Peter Malewski <p.malewski at tu-bs.de> Til: E. S. Venkatraman <venkat at biost.mskcc.org> Cc: r-help at stat.math.ethz.ch <r-help at stat.math.ethz.ch> Dato: 24. februar 2000 22:48

http://www.nabble.com/file/p24173983/groups.csv groups.csv http://www.nabble.com/file/p24173983/groups.r groups.r Hello, Points of an xyplot can be identified to a particular group according to 'sunflowers', a different one per group. However, this distinction does not appear in the legend generated by auto.key, all groups are expressed as a point or line, of different colours, but

Hallo, I've produced via R-0.62 several Life-tables 50*6 (mfrow=c(2,3)) as postscriptgraphs. The phd student, who asked for this, use word for windows--therefore i converted the whole stuff with pstoimg to gifs. Unlikly the result is not sufficient (I was told)... Is there a tool to convert ps to wmf/emf ? I looked for it but haven't found anything... Thanks for an answer Peter

Dear R-helpers, with R-0.99a, Linux, Red-Hat 6.1 derivat (Halloween Linux) I have the following Problem with this code: graphics.off() close.screen(all=T) temp <- matrix(c(0 , 1 , .4 , 1, 0 , 1 , 0 ,.4), ncol= 4, byrow=T ) split.screen(temp) plot(1:100)

In the following special case, seq fails to give the right answer (which is 0 ) > seq(0,0,1) Error in if (dd < sqrt(.Machine$double.eps)) return(from) : missing value where logical needed For any other equal from and to , it works: > seq(1,1,1) [1] 1 The error occurs in the statement dd <- abs(del)/max(abs(to), abs(from)) of seq.default -- for obvious reasons.

>>>>> "WSt" == Werner Stahel <stahel@stat.math.ethz.ch> writes: [to me privately] WSt> I do not want to call this a bug, but it bugged me. Please try the WSt> following: WSt> tform <- sqrt(RADAI) ~ sqrt(RADAI.e) + TAGE.ej + SPITAL + ARZT + DAS28 + WSt> SJC + TJC + DGA + HAQ + PGA + PAIN + YEAR.SYM + YEAR.DIA + WSt> ALTER +

symbols(iris[,1],iris[,2],rectangles=as.matrix(iris[,3:4]),inches=0.3, ylim=c(1,5)) Error in plot.default(NA, NA, type = "n", ylim = ylim, xlim = xlim, xlab = xlab, : formal argument "ylim" matched by multiple actual arguments Since plot is called inside symbols, with specified xlim and ylim, any specification through ,... hurts this call. I am pretty sure that you

I have just fitted a lm model containing a bs smooth term. I asked Martin how I could see the estimated function. He said that such a possibility -- which is available in S-Plus' plot.gam -- was on your list of desirable functions. This is a vote to give it high priority. (I appreciate all your efforts!) Werner Stahel ------------------ This message was sent by

Dear R-list, Is there any way I can extract part or whole column name from two variables and combine them? Not values, but column names. For example, I have data set with 70 variables. The column names are "Alfalfa", "Corn","Soybean",..."Sunflower". After I combine the two variable (say Alfalfa and Corn), I want to rename it as "AlfCor" or

Sehr geehrte Damen und Herren, An der Fachhochschule Winterthur wurde f??r mich im Rahmen einer Diplomarbeit (Prof. Ruckstuhl) ein R-Programm zur Auswertung von IR-Spektren (line-shape-Analyse zur quantitativen Analyse von Gasmischungen) entwickelt. Die Daten werden einer Excel-Tabelle entnommen. Bisher werden die spektroskopischen Daten ??ber den Inhalt von 64 Kan??len ??ber die

Dear fellow R users, I would like to draw a "sunflowerplot" because I have data (decade by month) that plots multiple times on the same x-y co-ordinates. Further I would like to colour each of the points/sunflower leaves on the plot according to the group they belong to (i.e. which type of event each represents within that decade and month). I thought that this would be relatively

I'm having a bit of trouble finding and understanding the correct function to make numeric vectors to feed the sunflowerplot function. I have 33k points to show and I want to do better than the standard scatter plot. I gather that I need two vectors (x and y) of the same count containing the "center" value of each bin. FYI - I have two pieces of data: 1 - x axis - time in days of

Hello all, I'm trying to graph a scatterplot of a large (5,000 x,y coordinates) of data with the caveat that many of the data points overlap with each other (share the same x AND y coordinates). In using the usual "plot" command, > plot(education, xlab="etc", ylab="etc") it seems that the overlap of points is not shown in the graph. Namely, there are

Hi all -- So I'm trying to use lattice graphics, but I want to use a sunflower plot, which doesn't seem to be part of lattice. No problem, I put together the following code, which mostly works -- *except* for the first graph it generates. If it opens the graphic device, then it draws the xygrid, clears the device, then draws the sunflowerplot. All subsequent output operations work

Hello, I wrote a wrapper for symbols() that produces a bivariate bubble plot, for use when plot(x,y) hides multiple occurrences of the same x,y combination (e.g. if x,y are integers). Circle area ~ counts per bin, and circle size is controlled by 'scale'. Question: how can I automatically make the smallest circle the same size as a standard plot character, rather than having to

We're using gluster-2.0.1 and a nufa volume comprised of thirteen subvolumes across thirteen hosts. We've found today that there are some files in the local filesystem associated with the subvolume from one of the hosts that are not being seen in the nufa volume on any gluster client. I don't know how or when this happened, but now we have to do some work to get this gluster volume

I have to admit that I'm a bit ignorant about regular expressions...I have a problem with "gsub": under unix I can : > cat tmp.txt F.123 F.123 F123 F123 sed 's/F\./d/g' tmp.txt d123 d123 F123 F123 that is: replacing "F.", however under R: gsub("F\.","d",c("F123","F.123"),extended=F) [1] "d23"

Does anybody has a function or idea how to calculate qt with a noncentrality parameter? I'm porting a SAS-macro for equivalence (Wellek, 1994 TT2ST) to R and it seems that I need this to calculate tinv (SAS-command, I hope that this is R's qt,...I get the SAS-macro manual next week). thanks Peter -- P.Malewski, Limmerstr.47, 30451 Hannover, 0511-2135008 At work: