similar to: iterating over groups of columns

I have a 45 x 16 data frame consisting of dissimilarities among 10 colors, giving in each column the 45 = 10*9/2 pairwise judgments for one of 16 subjects. The rownames identify each pair of colors, e.g, "AC" = ("A","C"), and the pairs are ordered by columns in the lower triangle of each distance matrix. > helm.raw <-

Hola!! tengo un data.frame donde cada fila corresponde a un año y cada columna a un mes (De enero a diciembre) > head(valT) V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 1941 18.0 16.3 15.2 10.1 8.1 8.3 8.8 9.2 7.9 12.2 11.9 14.6 1942 17.2 15.9 13.6 11.6 8.7 6.2 6.4 7.2 9.7 12.0 14.1 16.7 1943 17.6 17.3 13.5 12.5 10.5 7.0 8.2 7.9 -999.9 -999.9

Dear all, I am attempting to explain patterns of arthropod family richness (count data) using a regression model. It seems to be able to do a pretty good job as an explanatory model (i.e. demonstrating relationships between dependent and independent variables), but it has systematic problems as a predictive model: It is biased high at low observed values of family richness and biased low at

R Users, I am new to R and have tried to figure out how to automate this process instead of using excel. I have read in this dataframe into r with read.table. I need to reshape the data from the first table into the format of the second table. TractID StandID Species CruiseDate DBHClass TreesPerAcre Carbon Stand 1 Loblolly Pine 5/20/2010 10 1.2 Carbon Stand 1 Loblolly Pine

Dear all, I'm a trainee statistician in a company and we'd like to understand svm mechanism, at first with simple examples. I use e1071 package and I have several questions. I'm working with data extracted from cats data (from R). My dataset corresponds to a completely separable case with a binary response variable ( Sex with 2 levels: F and M), two explanatory variables (Bwt

Hi all! I tried to perform Shapiro-Wilk test for my sample of 243 values. > Us [1] -10.4 -13.1 -12.2 38.1 -18.8 -13.3 -11.7 29.3 49.7 6.8 12.7 16.3 [13] 5.8 -0.7 -29.4 4.1 38.8 -1.4 8.8 15.6 32.9 -5.3 19.1 35.8 [25] 4.0 -1.5 0.6 -4.2 -10.0 -4.0 1.1 48.9 -21.0 -5.3 5.8 -10.8 [37] 21.9 8.2 -3.2 -3.9 -2.3 12.6 -4.7 -8.0 11.8 27.4 -9.5 -20.8 [49]

Hello, I have a matrix (class matrix) composed of GridCell (row and column). The matrix value is the beta diversity index value between two grids. Now I would like to get the average value of each GridCell. Please kindly advise how to make the calculation. Thank you. Elaine The matrix looks like (cited from Michael Friendly) I would like to get the average value of each color. Obs

Reading data with variable column widths. Here are several lines of a txt data set I would like to read. The number of variables is fixed at 13 . The problem is how to read the first variable when it can contain blank space-- for example " Alabama (Seasonally Adjusted)" , "St. Clair", etc. Alabama (Seasonally Adjusted) 2,168,870 2,162,604 2,122,787 1,954,895 1,956,026

Hello, I'm using the "sem" package to do a confirmatory factor analysis on data collected with a questionnaire. In the model, there is a unique factor G and 23 items. I would like to calculate the standardized residual variance of the observed variables. "Sem" only gives the residual variance with the "summary" function, or the standardized loadings with the

On Thu, Apr 04, 2019 at 11:52:46AM -0400, Michael S. Tsirkin wrote: > I simply love it that you have analysed the individual impact of > each patch! Great job! Thanks! I followed Stefan's suggestions! > > For comparison's sake, it could be IMHO benefitial to add a column > with virtio-net+vhost-net performance. > > This will both give us an idea about whether the

On Thu, Apr 04, 2019 at 11:52:46AM -0400, Michael S. Tsirkin wrote: > I simply love it that you have analysed the individual impact of > each patch! Great job! Thanks! I followed Stefan's suggestions! > > For comparison's sake, it could be IMHO benefitial to add a column > with virtio-net+vhost-net performance. > > This will both give us an idea about whether the

http://www.nabble.com/file/p20332968/spe.dat spe.dat Hi, i have a problem in sorting out some command in R and i am really hoping some expert can help me out please! i have the spe.dat file which i upload here, and when u read into R u got something like this: A5 <- read.table('spe.dat' ,header=TRUE) A5 EXPEND ECAB MET GROW YOUNG OLD WEST STATE 1 256 85.5 19.7 6.9 29.6

Dear list, this works fine: x <- split(iris, iris$Species) x1 <- lapply(x, function(L) transform(L, g = L[,1:4] * 3)) but I would like to multiply each Species with another factor: setosa by 2, versicolor by 3 and virginica by 4. I've tried mapply but without success. Any thoughts? Thanks for any idea! Patrick

I have a set of modeled climate data recorded at irregular intervals. The format of the data is such that there are monthly measurements for the years 2000, 2020, 2050, 2080, etc. Therefore I have 12 regular records, a skip of some number of years, then 12 more monthly records, another skip, and so on.... I created a dataframe from .txt with the read.table() command. For starters I need

I wrote a simple log likelihood (for the ordinary least squares (OLS) model), in two ways. The first works out the likelihood. The second merely calls the first, but after transforming the variance parameter, so as to allow an unconstrained maximisation. So the second suffers a slight cost for one exp() and then it pays the cost of calling the first. I did performance measurement. One would

Hello all, I'm trying to use confint from the MASS package to compute confidence intervals for an nls object. When I plot the results, however, they don't make sense - lines cross over the fitted model or just don't match the data. Code is : Thanks for help dat<-data.frame(a,b) with(dat, plot(a,b)) model<-(nls(b~(1/exp(a*x))*n, data=dat, start=list(x=.001,n=20),

On Thu, Apr 04, 2019 at 06:47:15PM +0200, Stefano Garzarella wrote: > On Thu, Apr 04, 2019 at 11:52:46AM -0400, Michael S. Tsirkin wrote: > > I simply love it that you have analysed the individual impact of > > each patch! Great job! > > Thanks! I followed Stefan's suggestions! > > > > > For comparison's sake, it could be IMHO benefitial to add a

Dear all, I am trying to fit a model with crossed random effects using lme. In this experiment, I have been measuring oxygen consumption (mlmin) in bird nestlings, originating from three different treatments (treat), in a respirometer with 7 different channels (ch). I have also measured body mass (mass) for these birds. id nest treat year mlmin mass ch hack 1EP51711 17

I am aware of the inherent risks of having plots with more than two axes, but I am trying to produce the graphs that I have been tasked with. That being said I am having a hard time figuring out how to have two axes onto a boxplot. below is the sample code. I would like BC on the plot produced with this code to be on a second axis with all of the others being on the first axis. This will be

Dear all I am intrested in making iterative estimation (thro' loop statements) of, say, linear regression model. For this purpose, I have written the following programme and that I have made use of a sample data (viz., exp.txt): ? Programme: ? # Linear regression modelling with sample data (try5.txt) # Repeated estimation through loop statement