similar to: procrustes (vegan) plot of residual differences

Hello I will appreciate your help with the following. Running a script in R 2.14.1 under windows vista I get the following error message: Error in ifelse(append, "a", "w") : (list) object cannot be coerced to type 'logical' However, the very same script runs perfectly well under Ubuntu. My understanding is that this type of error is associated to

Hello I would appreciate your help on the followig. I want to generate random binary matrices but I need to discard those with all-1 rows. That is, for a 10x10 matrix with five 1's [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 0 0 0 0 0 0 0 0 0 1 [2,] 0 0 0 0 0 0 1 0 0 0 [3,] 0 0 0 0 0 0 1

Dear List, I am writing a custom panel function and xyplot method to plot the results of a procrustes analysis from the vegan package. I am having trouble getting the call to panel.arrows to work as I wish when conditioning. The attached file contains the function definitions for the xyplot method and the custom panel and prepanel functions I am using. This example, using data and functions from

If you submit these lines, you end up with variable "vare.dis". I want to export vare.dis to csv. Stuck I am. library(vegan,logical.return = TRUE) #return=true verifies package is available library(MASS,logical.return=TRUE) #return=true verifies package is available data(varespec) #varespec is an example data file in the vegan package vare.dis <- vegdist(varespec)

Dear all, I realised a correspondence analysis with function cca() of vegan library. Just like in Okansen (2010) in the example of R help: library(vegan) data(varespec) data(varechem) vare.cca<-cca(varespec~ Al + P + K, varechem) With plot.cca() function I represented the species matrix in the next way: plot(vare.cca,display="species") Being similar to: plot((c(-2,2)),(c(-2,2)),

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >

I performed canonical correspondence analysis (cca) with the example data of vegan, but I'm not able to obtain a table like scores() for the constraining variables. I can see them in the summary() mode, but it would be great to have in a separate table. Any suggestion?, thanx Gianandrea require(vegan) data(varespec) data(varechem) vare.cca<-cca(varespec,varechem) scores(vare.cca)

Hi I have just (finally) started to poke around in R and wanted to analyse a stream fish dataset with 28 sites and 18 species. When trying to follow the Vegan manual to run nmds from distance measures calculated by the vegdist function it turns out that I have two sites (streams) with the exactly the same four species (I have used pres-abs data in this case). When I try to run isoMDS I get an

Let's say I have a program that returns variables whose names may be any string within the vector NAMES=c("varA","varB","varC","varD","varE","varF"..."varZ"), but I do not ever know which ones have actually been created. So in one example output, "varA", "varC", and "varD" could exist, but

Dear R-list, I have a data set (in the following example called "a") which have: one "subject indicator" variable (called "id") three dependent variables (varD, varE, var F) three independent variables (varA, varB, varC) I want to fit 9 lme models, one per posible combination (DA, DB, DC, EA, EB, EC, FA, FB, FC). In stead of writting the 9 lme models, I want to

HI Guys, I know that this forum is not for homework but I am trying to interpret R output code. I was just wondering if someone might be able to help. I have been given the following. For (X1,X2) distributed bivariate normal with parameters mu1 = 5.8 mu2 = 5.3 sd1 = sd2 = 0.2 and p = 0.6 The r-code and inpit/output are as follows input m <- 5.3 + 0.6*(6.3 - 5.8) s <-

I don?t see how neither a) or b) applies to this question nor the technical merit of the remark about mixture models. Do you have a suggestion for a more appropriate forum for this issue/question? (stackoverflow basically sent me here). Kind regards > On 29. Jun 2017, at 16:58, Bert Gunter <bgunter.4567 at gmail.com> wrote: > > Offlist, because this is (a) an opinion and (b)

Would package "teigen" help? Ranjan On Thu, 29 Jun 2017 14:41:34 +0200 vare vare via R-help <r-help at r-project.org> wrote: > Hello! > > I am new to R (before used python exclusively and would actually call the R solution for this issue inside a python notebook, hope that doesn?t disqualify me right of the batch). > > Right now I am looking for a piece of

Offlist, because this is (a) an opinion and (b) about statistics and therefore offtopic. I don't know whether any such package exists, but I would predict that this is likely to be overdetermined (too many parameters) and therefore unlikely to be a successful strategy. Fitting a mixture of Gaussians is already difficult enough. Feel free to ignore, of course, and no need to reply. Cheers,

Hi... I have the next problem.. I have a FXO card with i can make calls but i cant recive calls. At the consol, i get the next error: -- Zap/2-1 is ringing -- Zap/2-1 is ringing -- Zap/2-1 answered Zap/1-1 -- Attempting native bridge of Zap/1-1 and Zap/2-1 WARNING[262160]: File chan_zap.c, Line 2857 (zt_handle_event): Ring/Off-hook in strange state 6 on channel

Add this: --delete-excluded See the manual. Rainer -----Urspr?ngliche Nachricht----- Von: rsync-bounces+rsync=diplan.de@lists.samba.org [mailto:rsync-bounces+rsync=diplan.de@lists.samba.org]Im Auftrag von Mandrake Pe Gesendet: Mittwoch, 29. Oktober 2003 13:33 An: rsync@lists.samba.org Betreff: trying to figure out how the --delete option works. I use this comand to back up my entirely drive

Hi, I would like to know whether it is possible to obtain a value of significance for random effects when aplying the lme or related functions. The default output in R is just a variance and standard deviation measurement. I feel it would be possible to obtain the significance of these random effects by comparing models with and without these effects. However, I'm not used to perform

Hi, I want to calculate the t- and p-values for a linear model using the Newey West estimator. I tried this Code and it usually worked just fine: > oberlm <- lm(DYH ~ BIP + Infl + EOil, data=HU_H) > coeftest(oberlm, NeweyWest(oberlm, lag=2)) t test of coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.1509950 0.0743832 2.0300 0.179486 BIP

Hii rstudio members I am learning rstudio, For my manuscript I am trying to plot CCA using species and environmental data. But I am getting error like Error in cca.default(sptrans, envtrans) : all row sums must be >0 in the community data matrix *My code is like * library(vegan) library(ggplot2) library(dplyr) rassspec<-read.csv("C:/Users/hp/Desktop/R_data/rassspec.csv",

Hi, trying to do a robudt regression of a two-way linear model, I keep getting the following error: > lqs(obs ~ y + s -1,method="lms", contrasts=list(s=("contr.sum"))) Error: lqs failed: all the samples were singular Robust regression with M-estimators works (also regular least square fits, of course): rlm.formula(formula = obs ~ y + s - 1, method = "M",