similar to: proj() function for lm objects

Hi all, I'm fitting a linear model (using lm) to some 2500 data points. The model consists of 4 single terms and two combined terms. I get the following warning message: "Extra arguments projections are just disregarded. in: lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) " Can anybody clarify this ? I don't seem to find any pointer to what this might

Hi all, I am faced with the situation where I want to store/analyze relatively large, organized sets of numerical data, which depend on a number of conditions (biological properties, exposure times, concentrations etc etc). Imagine about a hundred dataframes of a few thousand numerical values, with some annotation in text for some entries. Intuitively, I would like to be able to slice

Hi all, Never really managed to build a for-loop with multiple running variables in an elegant way. Can anybody hint ? See below for an example of what I would like. EXAMPLE a<-c(1,2,3) b<-c("name1","name2","name3") for( number in a, name in b ) { print( number ) ##take a value print( name ) ##and have its name available from a second list } Does R

Hi all, I need to enlarge te fonts used oo R-plots (plots, histograms, ...) in labels and titles etc. I seem to be unable to figure out how to do it. The problem is that the titles of the plots are simply unreadable when I insert them into my LaTeX text, since they are relatively small compared to the entire plot. I am sure it is pretty simple, can anybody give me a hint ? Please reply

Hi all, I'm new to R, and have the following problem: I wish to draw a boxplot of simple data in two columns. The x-axis (taken from first column) is grouped to intervals (using 'cut'). These intervals serve as x-value to the boxplot, and the data from the second column are the y-values. The problem is that I want to give an impression of the trend of the data in the x-range.

Hi all. I have a re-occuring typical problem that I don't know how to solve efficiently. The situation is the following: I have a number of data-sets (A,B,C,...) , consisting of an identifier (e.g. 11,12,13,...,20) and a measurement (e.g. in the range 100-120). I want to compile a large table, with all availabe identifiers in all data-sets in the rows, and a column for every

Hi, Whats the easiest way to set a desired interval for tics on axes in R ? And will the 'grid' command put gridlines on all tics automatically ? I tend to get stuck with graphs with ranges 0-1 with only 2 tics and 2 gridlines in the 0-1 range, while I would like 10 tics (or every 0.1) tnx, Piet please reply to pvremortNOSPAM at vub.ac.be

Hi all, Perhaps this in an inappropriate post, but I've found a bug in proj I'd like to track down a bit further before making a formal bug report. The example below shows the problem, if you change the rownames proj fails. The problem seems to be that there is a mismatch in the rownames in the qr objects constructed by aov and the rownames that proj is expecting them to have.

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 905 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010420/61ada248/attachment.bin

Dear All, I need to convert some Northing-Easting coordinates from the Gauss-Krüger system into latitude-longitude. Any suggestions on how to do it? Regards, Marco Marco Giannitrapani Statistical Consultant Tel: +44151373 5945 Email: Internet: http://www.shell.com [[alternative HTML version deleted]]

Hi R-user, I am new with R and I have a problem with the map and mapproj fonctions : with the following code : >test_proj=mapproject(LONG_d01_vec_t1, LAT_d01_vec_t1, projection="lambert", parameters=c(30,60)) >longitude_vec=test_proj$x >latitude_vec=test_proj$y >longitude_mat=matrix(longitude_vec,ncol=e_sn_d01) >latitude_mat=matrix(latitude_vec,ncol=e_sn_d01)

Hi, Could any one point me to the projection, and parameters if necessary, that would show each country/continent with it's area accurately refelcted on the plot? E.g. aitoff vs. albers vs. bonne vs. cylequalearea vs. guyou - they don't all look the same to mee but some of the documentations suggests they are equal area? Of course this isn't my field, so I am largely guessing and am

Hi, I am using in my workflow gIntersection from sp package. Part of my relevant sessionInfo is: R version 4.0.3 (2020-10-10) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 10 x64 (build 17763) attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] CircStats_0.2-6 boot_1.3-25 MASS_7.3-53 stringr_1.4.0

I'm using R 1.8.1 on Red Hat Linux 9. I've worked out the linear model decomposition by hand. Using "aov" with "Error" to fit a split-plot model, I get very different results depending on whether I use the "data" argument for "aov" or not. When I use the "data" argument, the ANOVA table from "summary" is correct but the

Using R 2.5.1 on Windows XP Professional, and PBSmapping package version 2.51, I have encountered some behavior which puzzles me. I am including the package's listed maintainer on this email but also seek the thoughts of the R-help community. I have a set of EventData, which I want to plot as points, and to color the points according to some criterion. It turns out that some of my

I am trying to set up a Grass project and need to set up the region so that I can view the map. I can look at a map and find the lat/lon, but the map projection is in UTM NAD38 WGS84 and I need to set the eastings and northings. Is there a package that will help me calculate this in R. thanks -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little

Hi List, I have a table that contains several fields that must be initialized when first created, but may not be altered afterwards (when updating the other fields) because they can be changed by external processes. I noticed that the ActiveRecord update commands (even update_attribute) only seem to work by loading the complete record, changing fields that need to be changed and doing a full

Hi R-user, I was trying to make a raster map with WGS84 projection in R, but I could not make it. I found one data set in Google that data is almost the same format as of mine. I wanted to make a raster map of temperature with 1 degree spatial resolution for the global scale. I could make it in GIS software but I do have many variables (to be many raster images) and ultimately I am importing them

Hello, A made a xyplot using the lattice library in R (latest version). The publisher of our paper has requested: 1. all tick marks should point inwards instead of outwards. 2. All lines should be thicker (lines, axes, boxes, etc. Everything). Lines is easy...I used: lwd=1.5 but what about the lines of the axes, and the lines that build up the plot itself?....? Any

Hi all, By default a call to xyplot from the Lattice package when using 2 factors [eg xyplot( dv~iv | XY * AB ) ] yields the following shingle structure: |_A_|_A_|_B_|_B_| |_X_|_Y_|_X_|_Y_| However, I'm wondering if it is possible to merge the upper shingle within levels of that factor, as in: |___A___|___B___| |_X_|_Y_|_X_|_Y_| Mike -- Mike Lawrence Graduate Student, Department of