similar to: distance between two matrices

>From: Sundar Dorai-Raj <sundar.dorai-raj at PDF.COM> >Reply-To: sundar.dorai-raj at PDF.COM >To: Zhen Pang <nusbj at hotmail.com> >CC: r-help at stat.math.ethz.ch >Subject: Vectorising and loop (was Re: [R] optim "a log-likelihood >function") >Date: Wed, 29 Sep 2004 18:21:17 -0700 > > > >Zhen Pang wrote: > >> >>I also use

For the following code below, the x-axis ticks are 1,2,3,4,5,6,7 when I was expection them to be 1,2,8,9,10,11,12. Please help me figure out where is the mistake. library(lattice) testdata <- as.data.frame(t(structure(c( 1,2005,9.24,6.18,634, 2,2005,8.65,6.05,96, 8,2004,6.81,6.51,16, 9,2004,9.0,7.29,8, 10,2004,8.84,6.18,524, 11,2004,8.54,6.35,579, 12,2004,9.97,6.3,614, 12,2005,8.75,5.84,32,

I've just encountered a segfault when using DierckxSpline::percur function. Below is the minimal example which triggers the error: --- library(DierckxSpline) x <- 1:10 y <- rep(0, 10) pspline <- percur(x, y) --- *** caught segfault *** address (nil), cause 'memory not mapped' Traceback: 1: .Fortran("percur", iopt = as.integer(iopt), m = as.integer(m), x =

Hi When I use plot(hclust(dist..)...)...) etc to create a dendrogram of a hierarchial cluster analysis, I end up with a vertical tree. What do I need to do to get a horizontal tree? Also, my users are used to seeing trees who's leaves all "end" at the same place (eg. Like in minitab). Is this possible in R? Thanks Mick Michael Watson Head of Informatics Institute for Animal

Full_Name: Sundar Dorai-Raj Version: 2.8.0 OS: Windows Submission from: (NULL) (76.220.41.126) The following code crashes R: library(tcltk) tt <- tktoplevel() tc <- tkcanvas(tt, yscrollcommand = function(...) tkset(ts, ...)) > sessionInfo() R version 2.8.0 (2008-10-20) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United

Dear useRs, I'd like to produce some scatter plots where N units on the X axis are equal to N units on the Y axis (as measured with a ruler, on screen or paper). This approach x <- sample(10:200,40) ; y <- sample(20:100,40) windows(width=max(x),height=max(y)) plot(x,y) is better than plot(x,y) but doesn't solve the problem because of the other parameters (margins etc). Is

Hi all, Here's a curosity I ran into since upgrading to 1.6.2 (precompiled for Windows). When using the chm help I get the following warning. I saw a recent post regarding this as a new warning (``dyn.load warning message in R1.6.2 on Windows XP'' dated 1/28/03), but not in the context of the help system. The warning only appears once and does not prevent the chm file from

I contacted John Fox about this first, because parts of the file are attributed to him. He says that he didn't write rstandard.glm(), and suggests asking r-devel. As it stands, rstandard.glm() has summary(model)$dispersion outside the sqrt(), while in rstandard.lm(), the sd is already sqrt()ed. This seems to follow stdres() in VR/MASS/R/stdres.R. Of course for the c("poisson",

Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason Dr. Iasonas Lamprianou Department of Education The University of Manchester Oxford Road, Manchester M13 9PL, UK Tel. 0044 161 275 3485 iasonas.lamprianou@manchester.ac.uk ___________________________________________________________ now. [[alternative HTML

I just noticed the following in update.formula and I'm wondering if this behavior is the intention of the developers. Here's an example: update(`a: b` ~ x, ~ . + y) Note now that the response has no backticks and is interpreted as a:b (i.e. ":" is now an operator). This is because in update.formula the call to terms.formula uses simplify = TRUE. I'm working with data that

Hi, Anyone can tell me how to generate "A01", "A02", ..., "A99"? paste("A", 1:99, sep="") generates "A1", "A2",..., "A99". This is not what I want. Thanks for the help. -MY [[alternative HTML version deleted]]

Hi, All, I'm writing a wrapper for stop that produces a popup window using tcltk. Something like: error <- function(...) { msg <- paste(..., sep = "") if(!length(msg)) msg <- "" if(require(tcltk, quiet = TRUE)) { tt <- tktoplevel() tkwm.title(tt, "Error") tkmsg <- tktext(tt, bg = "white") tkinsert(tkmsg,

I'm new to lattice and can't figure out what the problem is with the following example: ######################### > library(lattice) Loading required package: grid Attaching package `lattice': The following object(s) are masked _by_ .GlobalEnv : xyplot The following object(s) are masked from package:base : levels > test.data <- data.frame(x=rnorm(100), +

Hello, i know that i have to use optim, but i'm confused how its possible maximize the sum over all l[i] and get the optimized max(LL), r and alpha? LL <- function(trans,time){ for(i in 1:length(trans){ l[i] <- log(lgamma(r+trans[i] - gamma(r+1)*(alpha/alpha+t[i]))**r)*(t[i]/alpha+t[i]))**trans[i] } return(sum(l)) } i'm confused how i have to set r and alpha and i

Hi group, I'd like to do spatial analysis of my data using the spdep package. It appears that a file of class nb is necessary, but I do not find what that should look like and if there is a function that creates such file for me. How can I create a nb-object of my data points? Thanks, René Eschen. [[alternative HTML version deleted]]

I have two data sets: File1.txt: Name id1 id2 id3 ... N1 0 1 0 ... N2 0 1 1 ... N3 1 1 -1 ... ... File2.txt: Group id1 id2 id3 ... G1 1.22 1.34 2.44 ... G2 2.33 2.56 2.56 ... G3 1.56 1.99 1.46 ... ... I like to do: x1<-c(0,1,0,...) y1<-c(1.22,1.34, 2.44, ...)

Hi all, I have a question that I don't have a good answer for (note the word "good"; I have an answer, but I consider it not "good"). Take the following data in a single tab-delimited text file: <text> A Labels Value SE 2.5% 97.5% R90 0.231787 1.148044 0.035074 1.531779 R0 0.500861 0.604406 0.185336 1.353552 B Labels Value SE 2.5% 97.5% (Intercept) 1.367514

Hi, all, I'm looking a utility for removing a directory from within R. Currently, I'm using: foo <- function(...) { mydir <- tempdir() dir.create(mydir, showWarnings = FALSE, recursive = TRUE) on.exit(system(sprintf("rm -rf %s", mydir))) ## do some stuff in "mydir" invisible() } However, this is assumes "rm" is available. I know of

Hello, Is there anyone who could give me an example of how to import an ascii grid (i.e. ArcGIS exported raster) into R. I want to use it with gstat but don't know the appropriate import routine. Thanks very much for your help. Regards, femke Femke Reitsma Graduate Student (ABD) Geography Department 2181 LeFrak Hall University of Maryland College Park, MD 20742 Phone: 301-405-4121

Hi all, I have a data.frame with the following colnames pattern: x1 y11 x2 y21 y22 y23 x3 y31 y32 ... I.e. I have an x followed by a few y's. What I would like to do is turn this wide format into a tall format with two columns: "x", "y". The structure is that xi needs to be associated with yij (e.g. x1 should next to y11 and y12, x2 should be next to y21, y22, and