similar to: creating a 2d plan graphic and projecting points in it

I created a 2 D plan: > plot(range(-2.5,0.95),range(0.00,1.00),type="n",axes=TRUE) I made a projection of points with their coordonates (X,Y) in that plan > fullpoints=read.csv2("fullpoints.csv",h=T) > plot(fullpoints) The points are listed in that .csv file it is organized this way: 0,48875 0,142857143 0,409 0,142857143 0,45611 0,25 0,49833 0,222222222 0,61158

I have a bivariate plot of axis2 against axis1 (data below). I would like to use different size, type and color for points in the plot for the point coming from different region. For some reasons, I cannot get it done. Below is my code. col <- rep(c("blue", "red", "darkgreen"), c(16, 16, 16)) ## Choose different size of points cex <- rep(c(1, 1.2, 1), c(16, 16,

Hi, I would like to plot a graph on the jpeg device for example and to write a table associated to this plot where I have the pixel coordonates of each plotted point so that I could include the jpeg image onto a html page and get all the information about each point when the mouse button is pressed. The indentify() can do it but on the window device... Thanks for any idea. Philippe --

Hi, I wonder if anybody could help me in converting this easy SAS program into R. (I'm still trying to do that!) PROC IMPORT OUT= WORK.CHLA_italian DATAFILE= "C:\Documents and Settings\carleal\My Documents\REBECCA\stat\sas\All&nutrients.xls" DBMS=EXCEL2000 REPLACE; GETNAMES=YES; RUN; data chla_italian; set chla_italian;

?for doesn't return anything help.search("for") doesn't return anything- Is the for loop so prevelant in computer programing that the documentation is implicit or is R paradigm to discourage the use of the for loop. I will post data probably tonight, but here is my problem. I have preformed an MDS on a set of data. I have the scores of the four axes that are the optimal

Hello Followed the doc's to setup a postscriptprinter but get the following message when starting wine: err:psdrv:PSDRV_FindPrinterInfo Error 2 getting PPD file name for printer 'Epson_360_dpi_' fixme:winspool:AddPrinterW DocumentProperties fails err:psdrv:PSDRV_FindPrinterInfo Error 2 getting PPD file name for printer 'Epson_360_dpi_' err:winspool:AddPrinterW

Hi everyone. There's been several threads on baseline hazard in Cox model but I think they were all on cumulative baseline hazard, for instance http://tolstoy.newcastle.edu.au/R/help/01a/0464.html http://tolstoy.newcastle.edu.au/R/help/01a/0436.html "basehaz" in package survival seems to do a cumulative hazard. extract from the basehaz function: sfit <- survfit(fit) H

Dear R users, I'm currently trying to re-project a geotiff in another coordinate system. For instance, I have a tif image in UTM 19 zone which I would like to reproject into UTM 18. I was wondering if it was possible in R. Furthermore, I looked into 'rgdal' package, but I can't really find out if I'm doing the right thing. So far, here is what I'm doing: library(rgdal)

Hi, I have a dataset (CSV) with some counts of firms located around the globe. Each count is assigned to the longitude and latitude of the specific location. Now I want to plot these counts on a world map using dots (size of dots represent the count). I have been unable to find any info on whether this is possible and if so, how? Can you please help me? Thanks! Mathijs -- View this message

Hello all and thanks for the help. I am analyzing a dataset using MetaMDS and I would like to project some extra samples onto the plot such that the extra samples do not play a role in defining the axes. I have been thinking of different ways of doing this and I was wondering if anyone had a suggestion for an easy way to do this (possibly a way to weight samples in the analysis?) Thanks --

Hi R helpers, I tried to convert a list of LatLong coordinates (in DD format) into UTM zone 11 NAD 27. I first tried this from PBSmapping: library(PBSmapping) LatLong<-cbind(c(56.85359, 56.85478),c(-118.4109, -118.4035)) colnames(LatLong)<-c("X","Y") attr(LatLong, "projection") <- "UTM" attr(LatLong, "zone") <- 11

Hey all, I'm new to the forums, but not new to Wine/Linux. I'm having a bit of an issue. Playing with winetricks, I've gotten the Liteshow Manager install file to start and install up to 70% where it then fails. The application itself stays open, but it lets me know the installation fails. There's no logfile or explanation of why it fails. I'm able to wirelessly connect to

Hello, I have a survivor curve that shows account cancellations during the past 3.5 months.  Fortunately for our business, but unfortunately for my analysis, the survivor curve doesn't yet pass through 50%.  Is there a safe way to extend the survivor curve and estimate at what time we'll hit 50%? We started a new program 3.5 months ago, and I believe that this set of accounts behaves

I am working on creating an R package for doing fire department analysis and am trying to create a function that can display emergency incident densities. The following code sort of does the trick, but I need a display that shows the number of incidents per square mile. I believe the code below shows incidents per square unit (in this case, degrees lat/long). To solve this problem, I believe

Buenos Días, en esta ocasión les escribo para tratar de solucionar una duda que puede resultar muy sencilla, pero de la que no encuentro solucion. necesito hallar la descomposición en valores singulares; tengo la siguiente matriz de datos: DPTO PERSONAS HACINA INTERNET COMPU TV_SUCRIPC Antioquia 4962749,1 1,61 0,39 0,47 0,77 Atlántico 2191360,7 1,97 0,3 0,31 0,42 Bogotá D.C. 7499198 1,67

Estimados me ha surgido una duda con un simple stem and leaf el mismo es el siguiente: > stem.leaf(Datos$ej1.21) 1 | 2: represents 0.12 leaf unit: 0.01 n: 40 LO: 0.72 0.85 3 10 | 9 11 | 4 12 | 4 5 13 | 7 7 14 | 710 (*) 9 15 | 18 16 16 | 3447899 (6) 17 | 045599 (***) 18 18 | 2568 14 19 | 23710 (**) 10 20 | 389 7 21

Dear all, I'm analysing a split-plot experiment, where there are sometimes one or two values missing. I realized that if the data is slightly unbalanced, the effect of the subplot-treatment will also appear and be tested against the mainplot-error term. I replicated this with the Oats dataset from Yates (1935), contained in the nlme package, where Variety is on mainplot, and nitro on

I have data that looks like this: > df1 group id 1 red A 2 red B 3 red C 4 blue D 5 blue E 6 blue F I want a list of the groups containing vectors with the ids. I am avoiding subset(), as it is only recommended for interactive use. Here's what I have so far: df1 <- data.frame(group=c("red", "red", "red", "blue",

Hi, I am out of town and will get back to you on the 13th of July. Leo >>> "r-help at stat.math.ethz.ch" 06/27/03 00:32 >>> Send R-help mailing list submissions to r-help at stat.math.ethz.ch To subscribe or unsubscribe via the World Wide Web, visit https://www.stat.math.ethz.ch/mailman/listinfo/r-help or, via email, send a message with subject or body

I am projecting to migrate applications written with the Dialogic API to Asterisk. There are a few things that I do with Dialogic that I am still not sure are possible with Asterisk : - play a file from an offset (so a user can press a key to rewind 3 sec., pause, etc.) - dynamic volume control (the user can press a key to increase or decrease the volume) Thanks if anyone has a clue... Yves