similar to: non-metric multidimensional scaling

I measured nitrate concentration and primary production (PP) biweekly for 23 months in one headwater stream. I would like to use linear regression to determine if PP is related to nitrate concentration. My dataframe is called "data" and consists of the vectors Rdate, PP, and nitrate. Rdate is the observation date in class "date" and PP is primary production. I first

#is there a way to get NA in the table of descriptive statistics instead of the function stopping Thank you in advance #data x.f <- structure(list(Site = structure(c(9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L), .Label = c("BC", "HC", "RM119", "RM148", "RM179", "RM185",

I am trying to calculate monthly means by year of phosphates and nitrates from a multi year data set. Can anybody help me out with the most effective way to do this? My data looks like this: Collection_Date                   Test.Name              Value 2000-01-24 17:00:00            Phosphate               0.108 2000-01-24 17:00:00            Nitrate                     0.037 2001-11-12

Hello, my name is Toby, and i live in Sydney Australia. i would like to buy a 25Kg bag of strontium nitrate and have it delivered to an address in sydney. i understand that long distance postage can sometimes be very costly, so can you please give a a quote for the 25Kg bag, plus the postage cost to my address in sydney. thankyou for your time

Dear all, When analyzing data from a climate change experiment using linear mixed-effects models, I recently came across a situation where: - the summary(model) showed a significant difference between the levels of a two-level factor, - while the anova(model) showed no significance for that factor (see below). My question now is: Is the anova.lme() approach correct for that model? And why does

Dear list users, I have trouble with covariance analysis. I measured nitrate concentrations in the soil (NO3) and the percentage of legumes (LEG, continuous), affected by 2 different CO2 concentrations (CO2, discrete). I suspect that CO2 has an effect on LEG and NO3, but also that LEG has an effect on NO3, so this is the formula I wrote to test this: NO3 ~ CO2 + LEG + CO2:LEG Will LEG be

Hello, I tried to make a betadisper plot; however, it is quite messy at the moment with lines and symbols. I made two plots, one focusing on sites and the other on treatments. This is the code that I used: plot(betadisper(vegdist(y.nth,method="euclidean"),site)) plot(betadisper(vegdist(y.nth,method="euclidean"),treatment)) I have a few questions pertaining to how I could

Any suggestions on the best way to configure a puppet master replica/puppet HA/failover in different data center and subnet? Thought about masterless puppet, but seems like it is only good for some number of puppet agents. Thanks, Myra -- You received this message because you are subscribed to the Google Groups "Puppet Users" group. To view this discussion on the web visit

I have a question concerning the Wilcoxon signed-rank test, and specifically, which R subroutine I should use for my particular dataset. There are three different commands in R (that I'm aware of) that calculate the Wilcoxon signed-rank test; wilcox.test, wilcox.exact, and wilcoxsign_test. When I run the three commands on the same dataset, I get different p-values. I'm hoping that

Hello, Firstly, real new to R here. I have a function intended to evaluate the values in columns spread over many tables. I have an argument in the function that allows the user to input what sequence of tables they want to draw data from. The function seems to work fine, but when the user inputs a single number (over 9) instead of a sequence using the : operator, I find an error message: the

Dear list, I've been learning how to make a 2x2 paneled dotplot in lattice without any previous experience using lattice. my code thusfar is: nut<-read.table("/Users/colinwahl/Desktop/nutsimp_noerror.csv", T, sep= ",") attach(nut) nut1<-data.frame(Nitrate, Total_Nitrogen, Phosphate, Total_Phosphorus) nut1<-as.matrix(nut1) rownames(nut1)<-group

Hi, how are you? My name is Irucka Embry and I would like assistance on properly subsetting some data obtained from a .csv file. > file <-"info_n.csv" > INFO<- getMetaDataFromFile(file) INFO is 92 observations of 14 variables > INFO<- subset(INFO,site.no==02169570) INFO is 0 observations of 14 variables When the site.no is used I should have 1 observation of 14

Hi, I am trying to run a 1D nutrient-phytoplankton-zooplankton model in R using the package 'deSolve'. The code is shown below: DEPTH = seq(2.5, 147.5, 5) NPZ = function(t, state, params){ with(as.list(params), { P <- state[1:NB] Z <- state[(NB + 1): (2*NB)] N <- state[(2*NB + 1): (3*NB)] F.I = function(z, hr){ I0 = function(hr){

Does anyone have the code to calculate the drought severity index? -- Muhammad Rahiz | Doctoral Student in Regional Climate Modeling Climate Research Laboratory, School of Geography & the Environment Oxford University Centre for the Environment, University of Oxford South Parks Road, Oxford, OX1 3QY, United Kingdom Tel: +44 (0)1865-285194 Mobile: +44 (0)7854-625974 Email: muhammad.rahiz

hi i try to use arima and holtwinter to predict drought from 1895-2006 but i cannot read whole period of time and i try to do the exponent fitting, but it comes out as the coordinate x-y error i send the source code and data to take a look if anyone can help me, i am really new in R thank u so much jam

Hi, Could someone tell me if this is the correct model syntax for the following dataset: lme(height~treatment+genotype+treatment*genotype,drought,random=~genotyp e) The dataset has two factors: one fixed - treatment, and one random - genotype. I need to test the effect of both factors to identify their significance. There are multiple (but not equal) replicates at each level of genotype (the

Hello I am trying to fit a REML to some soil mineral data which has been collected over the time period 1999 - 2004. I want to know if the 19 different treatments imposed, differ in terms of their soil mineral content. A tree model of the data has shown differences between the treatments can be attributed to the Magnesium, Potassium and organic matter content of the soil, with Magnesium being the

Dear R-Help friends, I am unable to get the latest version of R (2.4.1) to compile on my solaris 10 system - has anybody else experienced this problem and are you able to offer me any advice? I appreciate your time, many thanks, Jenny Barnes Here are my CURRENT specifications: platform sparc-sun-solaris2.10 arch sparc os solaris2.10

Dear R-help community, I have 2 different arrays of precipitation data each of the same dimensions of [longitude, latitude, time] dim=[30,32,43], called array1 and array2. I need to correlate them. This is the code I used to get one overall correlation value for the whole of the area of interest: > result <- cor(array1,array2,use="complete.obs") > result This give me a

Dear R-help, Which random number generator function would you recommend for simply picking 15 random numbers from the sequence 0-42? I want to use replacement (so that the same number could potentially be picked more than once). I have read the R-help archives and the statistics and computing book on modern Applied statistics with S but the advice seems to be for much form complicated