similar to: Installing the Matrix package

Dear list, I have a series of data points which I want to approximate with exactly two linear functions. I would like to choose the intervals so that the total deviation from my fitted lines is minimal. How do I best do this? Thanks! Kamila The information transmitted in this electronic communication...{{dropped}}

Hi all! I'm afraid I programmed something totally non-sensical and inefficient, but I can't figure out how to do it better. I have a list of ~ 40 000 characters. I want to take each element at a time, map it to a large data frame with hit=which(data.frame$column==elementFromList), then compute some statistic on data.frame[hit,] and return a result that consists of either 1) a list

Hello, I have been following the thread dated Monday, October 9, 2006 when Kamila Naxerova asked a question about plotting elliptical shapes. Can you explain the equations for X and Y. I believe they used the parametric form of x and y (x=r cos(theta), y=r sin(theta). I don't know what r is here ? Can you explain 1)the origin of these equations and 2) what is r? Sincerely, Mary A. Marion

Is there a way to plot elliptical shapes? symbols() only provides circles... Thanks! Kamila

Hi everyone, disturbingly, a PCA yields different results when I run the analysis under Windows and Linux. I open the same workspace (which is on an independent partition) and perform the exact same function (by using the history file). Interestingly, PC1 is the same, but PC2 comes out totally different. Thanks! Kamila Windows: platform i386-pc-mingw32 arch i386 os

Hi! I don't understand this: layout(matrix(c(1:10),5,2),heights=c(1,rep(2,4))) plot(1,1) error in plot.new() : plot region too large Why?????? Thanks! Kamila

this is from a fresh svn checkout: [stvjc at stvjc1 R-devel]$ ls -tl Makefile -rw-rw-r-- 1 stvjc stvjc 12833 Jun 26 07:24 Makefile [stvjc at stvjc1 R-devel]$ svn up At revision 45988. [stvjc at stvjc1 R-devel]$ bin/R R version 2.8.0 Under development (unstable) (--) Copyright (C) The R Foundation for Statistical Computing ISBN 3-900051-07-0 .. > sessionInfo() R version 2.8.0 Under

Short form: How do I transform the output of grid.locator() (or grid.locator(unit='npc')) to the native (or npc) coordinates of a viewport other than the top-level viewport? Thanks in advance. -Ben Long form: I would like the user to be able to click anywhere on the R graphics window. I would then like to figure out which leaf viewport the user clicked in and compute the native

Hi, I have a seemingly common problem but I can't find a proper way to approach it. Let's say we have 5 samples (different size) of IC circuits coming from 5 production lines (A, B, C, D, E). We apply two different non-destructive QA procedures to each sample, producing to sets of binary outcomes (passed: no/yes). So, we have two groups of proportions: QA1 QA2 no/yes no/yes A

Hi, I am trying to make use of "approximate" option to obtain null distribution through Monte-Carlo resampling, as described in coin library documentation. Unfortunately, this does not play well with my data -- permutation process swallows astonishingly large amounts of RAM (4-5Gb) and runs far too long (30 min for B=10). Apparently, this is caused by the size of my dataset (see

Hello, I would like to display a heatmap with dendrograms and allow the user to click on either dendrogram to choose a sub-dendrogram. The function heatmap() uses layout() and image() and the plot method of dendrogram to plot the heatmap and dendrograms, but if locator() is then called and the user clicks, the coordinates returned from locator() are with respect to the last thing plotted (e.g.,

Hi all, I am trying to create simulated data for exploring reclassfication measures in a logistic setting with two continuous predictors and I would like to set the average population probability of outcome rather than the logistic beta0. Is there a way to find a beta0 that will generate the desired average population probability of outcome given x,y and their odds ratios? #Here is an

Hello, the man pages for openssh-2.3.0p1 look bad on our Tru64 Unix V4.0E systems. (Secure Shell) is a program for logging into a remote machine and for executing commands on a remote machine. It is intended to replace rlogin and rsh, and provide secure encrypted communica- tions between two untrusted hosts over an insecure network. X11 connections and arbitrary TCP/IP ports can also be

Hi Guys, Great example of how some of the OpenMoko guys are getting the word out there. If you don't know about OpenMoko you can review it here. http://deancollinsblog.blogspot.com/2006/11/fic-gta001.html http://deancollinsblog.blogspot.com/2006/11/open-phones.html http://deancollinsblog.blogspot.com/2007/06/home-brew-startrek-communica tor.html Otherwise check out some of the links

After installation of the latest version of R (2.6.0) I have a problem with loading the nlme library. When I use RGui under Windows XP, and try to install the nlme package, I get a warning that installation was unsuccessful due to a previously installed version of nlme. But the nlme library does not appear on the list of libraries that are available for loading under RGui. This happens after a

Hi, I am trying to replace (lattice) standard xyplot with xYplot variant from Hmisc package to be able to add error bars to my plots. However, this does not work, e.g: library(lattice) d <- data.frame( SKU=gl(3, 1, 21, labels=c("a", "b", "c")), Weekday=gl(7, 3, 21), QCRate=runif(21)) xyplot(QCRate ~ Weekday | SKU, data=d) (this plots nice 3 panels as per

I am receiving a list of name collisions when I launch R as seen below. I'm new to R and any suggestion or help with how I can go about getting rid of these collisions would be greatly appreciated. % R R version 2.5.1 (2007-06-27) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to

Dear all, it looks like fix() changes the class of mts objects, here is a reproducible example (tested both on WinXP and Linux): > x <- ts(cbind(1:5,1:5)) > x Time Series: Start = 1 End = 5 Frequency = 1 Series 1 Series 2 1 1 1 2 2 2 3 3 3 4 4 4 5 5 5 > class(x) [1] "mts" "ts" > edit(x)

Doing help("R_LIBS") brings up a help file (the same one as help(library)), but the help file doesn't mention R_LIBS. It does have a link to .libPaths, which does document R_LIBS. The quickest fix would be for help("R_LIBS") to bring up the .libPaths help file. --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = i386 os = mingw32

I was a bit puzzed by: > formatC(6.65,format="f",digits=1) [1] "6.6" So I experimented and found: > formatC(6.6500000000000001,format="f",digits=1) [1] "6.6" > formatC(6.650000000000001,format="f",digits=1) [1] "6.7" > round(6.6500000000000001,1) [1] 6.7 > round(6.650000000000001,1) [1] 6.7 > version