search for: cses

...version 3.1-24 of nlme in R 1.4.1 under Windows 2000, and both S-PLUS 2000 and 6.0, again under Windows 2000. I've noticed discrepancies in a couple of instances. Here's one, using data from Bryk and Raudenbush's Hierarchical Linear models: From R: > attach(Bryk) > cses <- meanses <- numeric(length(ses)) # initialize > for (sc in unique(school)) { + meanses[school==sc] <- mean(ses[school==sc]) + cses[school==sc]<-ses[school==sc]- meanses[school==sc] + } > Bryk$cses <- cses > Bryk$meanses <- mean...

...ame machine but R is much slower - like an order of magnitude - than S-PLUS on this problem. It looks curious to me when I run system.time() on this that timing vector shows almost no user time but a huge amount of system time. R> system.time(assign("fm2", + lme(MathAch ~ MEANSES*CSES + Sector + CSES:Sector, ttt, ~ CSES | School, + RML = TRUE ))) iteration = 0 Step: [1] 0 0 0 Parameter: [1] 2.61400 0.16608 3.88195 Function Value [1] 44964 Gradient: [1] -0.014379 -0.143176 -0.646558 iteration = 1 Step: [1] 0.002888 0.028756 0.129857 Parameter: [1] 2.61688 0...

Hi, Can someone please offer a procedure for going from CIs produced by intervals.lme() to fixed-effects response CIs. Here's a simple example: library(mlmRev) library(nlme) hsb.lme <- lme(mAch ~ minrty * sector, random = ~ 1 | cses, Hsb82) (intervals(hsb.lme)) (hsb.new <- data.frame minrty = rep(c('No', 'Yes'), 2), sector = rep(c('Public', 'Catholic'), each = 2))) cbind(hsb.new, predict(hsb.lme, hsb.new, level = 0)) Is the following correct (I know from the previous command that t...

Hi all, I''ve got the following two objects and the habtm relationship isn''t putting records in the DB: class Cse < ActiveRecord::Base has_and_belongs_to_many :conditions end class Condition < ActiveRecord::Base has_and_belongs_to_many :cses end I''m creating a Cse object and filling the collection of Conditions, then calling save on the Cse object. Nothing is getting saved in the conditions_cses table. Is there anyway to do fill the collection of conditions prior to saving the cse object? -- Posted via http://www.ruby-...

hi, I want to do hierarchical clustering with Jaccord index. I tried to do with vegan package for finding index and hierarchical clustering with hclust function. While doing clustering it is showing an error message as "invalid distance method". I would be grateful if anyone tells how to rectify the error. Thanks in advance,   kind regards, Ms.Karunambigai M PhD Scholar Dept. of

...ues display as simply numeric (i.e., 10,100,1000, ... on the X axis and proportions on the Y axis). Suggestions much appreciated. R 2.8.1 on mac OSX 10.4.11 Don McKenzie, Research Ecologist Pacific WIldland Fire Sciences Lab US Forest Service Affiliate Professor College of Forest Resources CSES Climate Impacts Group University of Washington desk: 206-732-7824 cell: 206-321-5966 dmck at u.washington.edu donaldmckenzie at fs.fed.us

Hello everyone, when I run the "unique" command on my data frame, it deletes the majority of duplicate rows, but not all of them. Here is a sample of my data. How do I get it to delete all the rows? 6 -115.38 32.894 195 162.94 D 8419 D 7 -115.432 32.864 115 208.91 D 8419 D 8 -115.447 32.773 1170 264.57 D 8419 D 9 -115.447 32.773 1170 264.57 D 8419 D 10 -115.447 32.773 1170

Does anyone know of any R code for computing partial cross-correlation? I have examples of cross correlation functions (ccfs) that are not smooth but rather consist of a peak of several high values in consecutive lags, with sharp drops on either side. This indicates that y(t) is a function of some average of x(t-tau) at the set of lags tau over which the ccf is high. I could sort out these

Hello, I have estimated an ARIMA model and I would like to make simulations from this estimated model 1,5 and 10 steps ahead. Does anybody know how to do this? Thank You Felipe Parra [[alternative HTML version deleted]]

Hi, I am trying to make a scatter plot with 4 different categories using the jitter function. My code returns a variable length error and will not plot because my four categories have different numbers of samples. When I delete samples from my spreadsheet so that each of the categories has the same number of variables, it plots just fine. Is there any way to get around this and plot all of my

...hank you for your time. -- ___________________________________ Don McKenzie, Research Ecologist Pacific Wildland Fire Sciences Lab USDA Forest Service 400 N 34th St. #201 Seattle, WA 98103, USA (206) 732-7824 donaldmckenzie at fs.fed.us Affiliate Assistant Professor College of Forest Resources CSES Climate Impacts Group University of Washington dmck at u.washington.edu

Hello everyone, When I run a for loop I noticed that the looping variable gets assigned the last value in the loop: For example, if I let h=200, and I run for(i in 1:h), i gets the value 200. What's wrong here? -- -Tony [[alternative HTML version deleted]]

Hi, The first column in as.matrix(x) has extra spaces (" "), which I don't want. Is there a way not to generate those spaces? Regards, Peng > x=data.frame(x=1:10,y=letters[11:20]) > as.matrix(x) x y [1,] " 1" "k" [2,] " 2" "l" [3,] " 3" "m" [4,] " 4" "n" [5,] " 5"

I am having a hard time interpreting the phase and coherency plots. x is two timeseries that occur at the same time i.e. a b 1 11.2 12.3 16 11.3 12.4 31 11.4 12.5 46 11.5 12.6 ...etc even though my example is does not show this they are oscillating at more or less the same frequency just shifted by t=x (imagine two sine waves offset with the 2nd sine

I'm trying to do this in several ways but havent had any result. Im asked to install python, or perl.... etc. Can anybody suggest a direct, easy and understandable way? Every help would be appreciated. Thx. -- View this message in context: http://r.789695.n4.nabble.com/Export-R-dataframes-to-excel-tp3330399p3330399.html Sent from the R help mailing list archive at Nabble.com.

This is has been bugging me for a long time. Nobody around me seems to have this problem. I hope someone on the forum could help me. When I generate a R graph and want to bring the image into Word. I cannot copy and paste it like many of my classmates can do. The "Select" in Edit menu have all options grayed out. So what I ended up doing is to save it as pdf file and bring it

Dear R experts, i have a vector z , i have to do something after z is sorted. how can i find the original index, i.e., before sorting, of a certain element in the sorted vector . thanks in advance regards ----------------------------------------------- Peng Jiang ?? ,Ph.D. Candidate Antai College of Economics & Management ???????? Department of Mathematics ??? Shanghai Jiaotong

Dear all, I've tried a multiple regression, and now I want to try a cross-validation. I obtain this error (it must be sth related to df) that I don't understand, any help would be appreciated. cv.lm(df= dat, lm2.52f, m=3) Error en `[.data.frame`(df, , ynam) : undefined columns selected lm2.52f is my lm object, dat is a dataframe where the variables involved in .lm are I tried CVlm

On Fri, Feb 27, 2015 at 2:13 PM Ahmed Bougacha <ahmed.bougacha at gmail.com> wrote: > On Fri, Feb 27, 2015 at 1:42 PM, Eric Christopher <echristo at gmail.com> > wrote: > > > > > > On Fri, Feb 27, 2015 at 1:38 PM Ahmed Bougacha <ahmed.bougacha at gmail.com > > > > wrote: > >> > >> On Thu, Feb 26, 2015 at 2:33 AM, Kristof Beyls

Dear R-Help Team, I have been trying to sort (all columns of) a matrix: > a<-matrix(a<-c(1,3,4,6,6,4,6,56,4,64,86,39,4,2),length(a),2) > a [,1] [,2] [1,] 1 1 [2,] 3 3 [3,] 4 4 [4,] 6 6 [5,] 6 6 [6,] 4 4 [7,] 6 6 [8,] 56 56 [9,] 4 4 [10,] 64 64 [11,] 86 86 [12,] 39 39 [13,] 4 4 [14,] 2 2 > I