search for: aquator

Hi, Fitting all possible models (GLM) with 10 predictors will result in loads of (2^10 - 1) models. I want to do that in order to get the importance of variables (having an unbalanced variable design) by summing the up the AIC-weights of models including the same variable, for every variable separately. It's time consuming and annoying to define all possible models by hand. Is there a

I have a specialist piece of scientific software that I am trying to get to run in WINE (0.9.59 running in Ubuntu 8.04). I have installed the software, but when I run it I get a message saying "Microsoft Scripting needs to be installed. Please install the latest version of Internet Explorer." Any ideas what to do?

Hello Instead of testing against 0 i would like to test regression slopes against -1. Any idea if there's an R script (package?) available. Thanks for any hint. Cheers Lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology ?berlandstrasse 133 CH-8600 D?bendorf Switzerland Phone: +41 (0) 71 220

heya, i am fitting linear mixed effect model to a response Y. Y shows an s-shaped distribution when using QQ-plots (some zero values and some very high values). hence, which transformation should i apply that Y follows a normal distribution? any r-function/package available to do this? thanks for any hint, regards, lukas ??? Lukas Indermaur, PhD student eawag / Swiss Federal Institute of

Hi, I would like to evaluate the frequency of the variables within the best selected model by AIC among a set of 12 competing models (I fit them with GLM) with a bootstrap procedure to get unbiased results. So I would ike to do the ranking of the 12-model-set 10'000 times separately and calculate the frequency of variables of the 10'000 best ranked models. I wrote a script doing the model

Hello I want to batch job the calculation of many GLM-models, extract some values and store them in a file. Almost everything in the script below works (read file, extract values and write them to file) except I fail in indexing the GLM with the modelstructure it should run. Running GLM's conventionally is no problem. Conventionally a GLM is calculated as:

Hi everybody, Well, I really hoped that until now i know how to do a selection ..... but obviously i don't. So risking your laughing at me - here you go. I have a table with lots of vegetation data, and one column specifies in which region the data is. There are only 2 regions, region 1 and region 3. I want to select all data that belong to region 1 and after do some xtabs on it. If i do

Hello All - Any recommendations or suggestions for neat ways to visualize data taken from a 3-axis accelerometer? My study species is aquatic, so I would be interested in movement in the 3 dimensions in addition to being able to incorporate the time series as well. Is there a package in R that might be useful for this? Thank you in advance, Jen ~~~~~~~~~~~~~~~~~~~~~ Jennifer Maresh, PhD

Hi, I'm a PhD student and I want to use the function 'daisy' from the package 'cluster' to compute dissimilarities. My variables are of mixed types so I use the argument 'stand' in daisy to define the type of my variables. I have the following error message : Warning message: binary variable(s) 13, 16, 17, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,

Buenos días. Tengo dos mapas que hago con ggplot (abajo) y necesito superponer el 2º al 1º. En el 1º se representa la variable PA, binomial (0, 1), en gris, y en el 2º la variable ExtTG7085, que es numérica y va de 0 a 200. La 2ª es un gradiente en rojo. Lo ideal sería que donde coinciden el rojo de la 2ª se viera más oscuro, pero supongo que eso no es fácil. Podría también ser en otro color, pero

Well.. plans are there from a long time to rewrite SciViews completely and make it platform independent (to work on Linux/Unix and MacOS X, as well as Windows). I have done some work in this direction when time permitted, but I am pretty busy with other work. During the holidays, I will continue to work in this direction. I will try to package a first running version of SciViews compatible with

Hi all, I have two datasets, one that represents a long-term time series and one that represents summary data for the time series. It looks something like this: x<-data.frame(Year=c(2001,2001,2001,2001,2001,2001,2002,2002,2002,2002,2002,2002), Month=c(1,1,1,2,2,2),Q=c(5,5,5,6,6,6,3,3,3,4,4,5)) y<-data.frame(Year=c(2001,2001,2002,2002),Month=c(1,2,1,2),Threshold_Q=c(5,5,4,4)) What I'd

A statistical question, not specific to R. I'm asking for a pointer for a source of definitive descriptions of what types of data are best summarized by the arithmetic, geometric, and harmonic means. As an aquatic ecologist I see regulators apply the geometric mean to geochemical concentrations rather than using the arithmetic mean. I want to know whether the geometric mean of a set of

Hello all, I have a *.csv file that looks like this (actual file is orders of magnitude larger): Site taxa no.ind forest LMA 1 forest LCY 1 forest SCO 1 meadow LMA 2 meadow LCY 1 meadow PNT

Estimate region of highest probabilty density Dear R-community I have data consisting of x and y. To each pair (x,y) a z value (weight) is assigned. With kde2d I can estimate the densities on a regular grid and based on this make a contour plot (not considering the z-values). According to an earlier post in the list I adjusted the kde2d to kde2d.weighted (see code below) to estimate the

Dear all, I am a new user of R. I would like to know how to get fund's returns in percentage (%). For example, I use: R <- ts(read.xls("FundData"), frequency = 12, start = c(1996, 1)) Whith this program, the returns are like 0.0152699. But, I would like to have 1.52%. Please advise me about the function. Thanks! Fabrice

> Background: > OS: Linux Ubuntu Dapper > release: R 2.3.0 > editor: GNU Emacs 21.4.1 > front-end: ESS 5.2.3 > ----------------------------- > Colleagues > I have two spatial datasets (latitude, longitude, fish eggs) and (latitude, longitude, fish larvae) at the same 280 stations (i.e. 280 cases). I want to determine if the 2 datasets are spatially correlated. In other

Buenos días. Hago mis mapas con el código que os copio abajo, pero me gustaría que me dibujase las líneas del mapa después de poner los puntos para que se vean. Puedo hacer los puntos más pequeños, pero entonces quedan separados y yo quiero que rellenen el mapa. Gracias, como siempre, Manuel print(ggplot(legend=FALSE)+geom_path(data=map_data('world'), aes(x=long,

Hi all, I'm trying to identify a particular digit or value within a vector of factors. Specifically, this is environmental data where in some cases the minimum value reported is "<" a particular number (and I want to manipulate only these). For example: x<-c("1","2","3","4","<1") For a dataset that is hundreds or

Hi all, I want to modify a large number of text files (ca 4000) by replacing a value found on a particular line in them with a value from an R object. For a single file I would normally use: con<-file ("foo.txt", open="r+") content<-readLines(con) content[n]<-"test" writeLines(content,con) close(con) For repeating this for several files I can