similar to: Different results under Windows and Linux

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! 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

Hi all, is the Matrix package no longer available for download via install.packages? When I try to install it (from any mirror), I get the following error message: install.packages("Matrix",lib="/home/kn52/R-2.5.0/library") Warning message: package 'Matrix' is not available in install.packages("Matrix",lib="/home/kn52/R-2.5.0/library") Other

Hi, I wanted to reply messages (e.g. to the cross-build-env in centos-devel) but my replies get posted as a new thread. I'm using Emulution 2.8 and I think Evolution eats the in-reply-to information. Is there anyone who also uses Evolution (especially in conjunction with a Scalix imap server)? So please tell me the right settings to acting not so disturbingly. Thanks. Brian.

After retrieving a URL using an Ajax.Request... function onComplete (request, json) { // request.responseText holds the entire page Only the content contained in elementID = ''source'' is intended for display At this time, the following code works, but it disturbingly copies the entire page from the source URL into the target DIV $(''target'').innerHTML =

----- Quick summary: In the File:Open dialog, please change "S files (*.q)" to "S files (*.q, *.ssc, *.S)" and show the corresponding files (including .SSC and .s files). ----- Background This is motivated by the following query to R-help: >Date: Thu, 16 Mar 2006 22:44:11 -0600 >From: "xpRt.wannabe" <xprt.wannabe at gmail.com> >Subject: [R] Is

R-Users, I had a strange occurrence this AM I could use some help with. 'var' seems to have ceased to work for me when I used JGR 1.3 interface. So, I ran a quick test (following pg 57 of Dalgaard) in both the basic RGui and JGR interface (see below). I have re-installed R 2.2.1 and JGR 1.3 but had no luck so far figuring this out? Anyone have a suggestion? Or did I just make a

R-help, Sorry for posting the same question (almost) once again. I'm using a lm model for some data. The code is below: fitData <- glm(pm ~ age + mw , family = binomial(link = logit), maxit = 30, data = inputModel1) "age" is a factor and "mw" is a continous variable. The summary function does not provide the parameter estimates for the first level of the

Trying to install the Matrix package with install.packages fails for me on Linux, it is trying to fetch the wrong version 0.99-4, downloading the tgz which is 0.99-6 and using R CMD INSTALL works fine The output from the failed install was > install.packages(c("Matrix")) trying URL 'http://cran.uk.r-project.org/src/contrib/Matrix_0.99-4.tar.gz' Error in

Legend titles work in linear plots: curve(1/x, xlim = c(0, 1)) legend(x = 'topright', inset = 0.04, legend = '1/x', lty = 1, title = 'Legend Title') But when you change to a log plot on either dimension things get screwy: curve(1/x, xlim = c(0, 1), log = 'y') legend(x = 'topright', inset = 0.04, legend = '1/x', lty = 1,

About a week and half ago I upgraded from 4.0.12 to 4.9.6. Overall, things are functioning. However, I have come across several strange behaviors and wandered if anyone else has noticed similar behavior on 4.9.6 or has any suggestions of what might be occurring. As background information, I have 3 DCs (dc3, dc4 and dc5) -- all running the same version (4.9.6) and all have the same configuration;

Dear List, The following code produces a list, which is what I what: set.seed(123) tmpf <- function() { x <- rpois(rpois(1,4),2) } n <- 3 m <- replicate(n,tmpf()) m [[1]] [1] 3 2 4 [[2]] [1] 0 2 4 2 2 5 2 [[3]] [1] 2 0 4 1 0 Now I need something that would to extract iteratively (or as many times as the size of 'n') the values that are greater than 2 in each component of

> version _ platform x86_64-unknown-linux-gnu arch x86_64 os linux-gnu system x86_64, linux-gnu status major 2 minor 2.1 year 2005 month 12 day 20 svn rev

Dear R users, My situation: (1) I have limited workspace for my work harddisk (about 10 GiB). (2) I have a lot of data files in R workspace (*.RData) which most of them > 200 MiB. For some reason I zip some of them, for instance "filename.RData (250 MiB)" to "filename.zip (3MiB)". In this work I have a lot of more space of my harddisk. Normally, If I want to use

Dear R-devels, I'd like to create a plot method for a class of objects that passes the '...' argument to both plot() and legend(), e.g., x <- list(data = rnorm(1000)) class(x) <- "foo" plot.foo <- function(x, legend = FALSE, cx = "topright", cy = NULL, ...){ dx <- sort(x$data) plot(dx, dnorm(dx), type = "l", ...) if (legend)

Given the following, one of the things I am trying to see is what % of draws are below a certain number: lambda <- 3 rate <- 5 n <- 5 set.seed(123) v <- replicate(n, rexp(rpois(1,lambda), rate)) vv <- unlist(v) cat("% of draws below 0.1:", round(length(subset(vv, vv < 0.1))/length(vv)*100,0), "%\n") In actuality, my lambda, rate, and n are 26, 10, 1000000,

Any help on the following would be appreciated... > R.version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 2.1 year 2005 month 12 day 20 svn rev 36812 > language R > set.seed(2) > x=matrix(ncol=4, nrow=10) > for(i in 1:4) x[,i]=rnorm(10) > x=as.data.frame(x) > x$V1=c(rep(0,5), rep(1,5)) >