similar to: selecting maximum values

I calculate an anova test in the following way: expdata<-read.table("/home/dorien/UA/meta-music/optimuse/optimuse1-build-desktop/results/results_processedCP", header=TRUE)

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;

Hi, I would be very grateful for some advice. I have read the help pages for Date, strptime, etc. All examples seem to use some version of day month year as date format. However I have Weekly composite data so ONLY want to input the dates as Week.Year (eg 35.2011). strptime seems to show this is possible using %W for week (UK convention) and %Y for year. My data is in a df called chlorophyll and

Hi, I have a problem. I have a data frame looking like: ID val A? .3 B? 1.2 C? 3.4 D? 2.2 E? 2.0 I need to CREATE the following TABLE: CASE?? DIFF A-A??? 0 A-B??? -0.9 A-C??? -3.1 A-D??? -1.9 A-E??? -1.7 B-A??? ... B-B??? ... B-C B-D B-E C-A C-B C-C C-D C-E D-A D-B D-C D-D D-E E-A E-B E-C E-D E-E WHERE CASE IS THE COUPLE OF ELEMENTS CONSIDEREDM AND DIFF IS THE computed DIFFERENCE between

I am having difficulty understanding the syntax of the nlme() function for nonlinear mixed models. The data frame is called Marouane.chlorophyll. The model involves a dependent variable (Absorb) and an independent variable (Ch.surf), which are both numeric variables in the data frame. The data are hierarchically grouped as Espece, Plante nested within Espece, and Feuille nested within Plante

Hi everybody, I'm trying to modify the boxplot just to set the upper whisker to the 90 percentile value, but I still couldn't find the solution. Can anyone help me? Thanks Alessandro Carletti

Hello all, I can't seem to figure how to use a greek character in expression() in plot() labels without adding a space. So for example below when plotting this out x<-1:10 plot(x,x^2, xlab=expression(Chlorophyll~italic(a)~mu~g~cm^-2)) the axis label read as μ g cm^-2 because I have space there with a tilda. But if I remove the tilda then my units are mug cm^-2. Can anyone recommend a

On 21/03/2018 11:44 AM, Stefano Sofia wrote: > Dear list users, > let me ask you this trivial question. I worked on that for a long time, by now. > Suppose to have a data frame with NAs and to sum some columns with rowSums: > > df <- data.frame(A = runif(10), B = runif(10), C = rnorm(10)) > df[1, ] <- NA > rowSums(df[ , which(names(df) %in%

Dear list users, let me ask you this trivial question. I worked on that for a long time, by now. Suppose to have a data frame with NAs and to sum some columns with rowSums: df <- data.frame(A = runif(10), B = runif(10), C = rnorm(10)) df[1, ] <- NA rowSums(df[ , which(names(df) %in% c("A","B"))], na.rm=T) If all the elements of the selected columns are NA, rowSums

Hi everybody, I'm performing a cluster analysis (pkg "cluster") on a dataset which includes 15 variables: is there a way to know how much each variable weighs on the final clustering output? Thanks Alessandro

Hi, I'm trying to convert a vector containing dates in character format ("dd/mm/yy"): mdy.date (from date package) seems to be able to do that, but it returns to me a vector containing julian dates... but negative! for example: 16/12/03 is converted into -20470 it is because R recognizes year ../03 as 1903, instead of 2003, but how can I do to solve this problem? (of course, I

Dear R users, I need to produce rainfall maps using R. I know that this is possible, I looked though the web, I found the example below reported (the author is Andrew Tredennick). I would ask you if this is the most performing way to make rainfall maps; if yes would someone be able to give me an example of how file.asc and pointfile.csv should be? If no would somebody please show me another way

Hello, Several machine cannot joint the cluster after a reboot. Here is a example of the error : caserta kernel: CMAN: Join request from como.dmz.lexum.pri rejected, node ID 3 already in use by ancona.dmz.lexum.pri All cluster.conf are the same (checked with md5sum) I googled this but did find nothing Any ideas? tx

Dear list users, I have to aggregate wind direction data (wd) using a function that requires also a second input variable, wind speed (ws). This is the function that I need to use: my_fun <- function(wd1, ws1){ u_component <- -ws1*sin(2*pi*wd1/360) v_component <- -ws1*cos(2*pi*wd1/360) mean_u <- mean(u_component, na.rm=T) mean_v <- mean(v_component, na.rm=T) mean_wd

Hi R-users I'm a relative newbie and uses R mostly for graphical purpose. I have a layout problem when plotting a scatterplot and a barplot using par(new=TRUE). The baseline of the x-axis is not positioned equal for the two plotting functions (see picture) and I have been unable to find out how this is changed. http://www.nabble.com/file/p18025066/pic.jpeg I have added the script if this

Should not the result be NULL if you have removed the NA with na.rm=TRUE ? B. > On Mar 21, 2018, at 11:44 AM, Stefano Sofia <stefano.sofia at regione.marche.it> wrote: > > Dear list users, > let me ask you this trivial question. I worked on that for a long time, by now. > Suppose to have a data frame with NAs and to sum some columns with rowSums: > > df <-

Hi, You might get more help from the R-sig-geo list, which is devoted to spatial topics. However. The *.asc file is an ArcGIS raster export format. You should use whatever the appropriate import commands are for your own gridded rainfall data. If you have a different format, you might or might not be able to import it directly with raster. ?raster will tell you more about the kinds of formats

What do you mean by "should not"? NULL means "missing object" in R. The result of the sum function is always expected to be numeric... so NA_real or NA_integer could make sense as possible return values. But you cannot compute on NULL so no, that doesn't work. See the note under the "Value" section of ?sum as to why zero is returned when all inputs are removed.

Hi everybody, I'd like to know if there's an easy way for extracting outliers record from a dataset, in order to perform further analysis on them. Thanks Alessandro

Fwiw the engine behind geom_raster needs explicit observation-per-row form for input (with no structural normalization), so conversion to points is perfectly proper here, albeit confusing in context. (It's closer to what graphics devices actually use ultimately, but the expansion is laid out very early in ggplot2 because there's no standard for intermediate forms.) Cheers, Mike On Wed,