Re: [R] SEM specify RAM model

Dear George, I am not very used to sem package, but it seems to me that you can use a file as an argument of specifyModel() or specifyEquations() functions, instead of input each line individually on the shell. see the first argument on ?specifyEquations Maybe you could write such matrix insid

Re: [R] RODBC and PosgreSQL problems

Dear All, I just wanted to follow up my question with an answer, which I owe to Robbie Bingler at UVA's IATH. The code chunk that bombed is here: sqlQuery(DRCch,paste(" + SELECT * + FROM tblCeramicWare + ")) [1] "42P01 7 ERROR: relation \"tblceramicware\" does not exis

Re: [R] Code for generating states and observations for HMM

Dear Dr Rolf, Thank you for your response and suggestion. But i still find it very difficult to relate my own problem to the sim.hmm() function you suggested. Because the states and the observations should come from Poisson and Negative Binomial distributions where state "0" to come from Poisso

[R] I need help computing PRESS statistics (qpcR package) of a nlsLM model (minpack.lm package)

Dear all, I have used nlsLM in minpack.lm package for fitting a model using 22358 observations. Now I would like to compute PRESS residuals and have been trying it using qpcR package but it does not seem to work (I have fitted my data with a lm model and used qpcR PRESS and it worked fine).

Re: [R] Code for generating states and observations for HMM

Sorry, I haven't the time at the moment to delve into this in detail, but at a quick reading, your question makes no sense. If you are trying to simulate a hidden Markov chain, the states have a distribution as determined by the transition probability matrix and the initial state probability

Re: [R] Confusing behavior when using gsub to insert unicode character (minimal working example provided)

Yep. You are right. That is better. -tgs On Thu, May 29, 2014 at 5:23 PM, Ista Zahn wrote: > 10Hi Thomas, > > On Thu, May 29, 2014 at 9:15 AM, Thomas Stewart > wrote: > > Thanks to to Ista Zahn, I was able to find a work around solution. The > key > > seems to be that string1 needs to be en

[R] Collapsing data.frame to its or xts

Hello, I have a data.frame of a time series sampled every 15 minutes. Also the data are in reverse order (recent to older). It looks like the following: Date 0:00 0:15 0:30 0:45 1:00 1:15 1:30 1:45 2:00 2:15 2:30 2:45 3:00 3:15 3:30 3:45 4:00 4:15 4:30 4:45 5:00 5:15 5:30 5:45 6:00 6:15 6:30 6:

[R] Preserving topology when simplifying Spatial Polygons

I'm trying to simplify a group of adjacent polygons without gaps and line overlaps forming between them. Ideally I'd like results similar to what mapshaper produces. I've tried using gSimplify() (from rGEOS), thinnedSpatialPoly() (from maptools), and dp() (from shapefiles) but haven't been able to

[R] SEM specify RAM model

Dear all, I am trying to create multiple SEM models in a loop - so that if I change the data they can easily be re-created rather than re-creating each RAM model line by line - therein lies my problem. I have a dataframe set out as follows1 Duration -> RC1 cov1 -0.2563002

[R] StarCluster + R

Hello Alex, I am wondering if you reached a solution to your question asked in 2011 about parallelizing R in conjuction with starcluster infrostructure? if there are some ready to use scripts to run a job inside StarCluster. Thank you, Best Regards, Heba El Kholy [[alternative HTML

Re: [R] tab auto-fill and arrow keys are not working in R on a linux cluster

Thank you, Jeff!  I'll try that. Tao On Tuesday, June 3, 2014 5:22 PM, Jeff Newmiller wrote: FAQ 7.20 --- Jeff Newmiller                        The    .      .  Go Live... DCN:        Basics: ##.#.      ##.#.

Re: [R] recode the ID with sequential order of a dataset

Hi York, Using "dat" (assuming that the data is ordered by "ID") dat$NEWID <- cumsum(c(TRUE,with(dat, ID[-1]!= ID[-length(ID)]))) A.K. Hi arun, Thanks for your reply. I think you misunderstood my question, for the new ID, I would like to have the following pattern, ID TIME NEWID 1254 0  

Re: [R] Selecting first and last row for each date

Hi, May be this helps: dat <- read.table(text="DATE Price 01.01.2010 2 01.01.2010 3 01.01.2010 2 01.01.2010 7 02.01.2010 3 02.01.2010 9 02.01.2010 0 03.01.2010 2 03.01.2010 4 03.01.2010 3 03.01.2010 6 03.01.2010 8",sep="",header=TRUE,strin

Re: [R] tkbind (callback)

And the path may vary a bit between installations within Windows, but 'somewhere in the R folder' is close enough to find it. -Original Message- From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk] Sent: June 4, 2014 2:53 PM To: Greg Snow; Fowler, Mark Cc: R help Subject: Re: [R] tkbind

Re: [R] How to extract data for the week ending on Friday from date column.

Hi, Using the example data, library(zoo) z1 <- read.zoo("dfRaw_20140509.csv",sep="",index.column=2,format="%Y%m%d") z1 #   P_num    P_ID Totalvalue #2013-05-10 25193 1230238  1.203 #2013-05-24 25190 1230238  1.201 #2013-07-05 25191 1230238  1.208 #2013-08-02 25194 1230238 

Re: [R] tkbind (callback)

On 04/06/2014 18:30, Greg Snow wrote: Installed with R is the documentation for tcl/tk (the original language), on my computer it is in the Doc folder/dir of the TCL folder/dir under the R folder. In that documentation on listbox, near That will only be the case on Windows. On other OSes you

Re: [R] tkbind (callback)

Installed with R is the documentation for tcl/tk (the original language), on my computer it is in the Doc folder/dir of the TCL folder/dir under the R folder. In that documentation on listbox, near the bottom in the "Default Bindings" section it mentions the virtual event "<>", so I tried that and

Re: [R] tkbind (callback)

Thanks for responding, Greg. That works. Plus I wasn't aware of a ListboxSelect option, which is much more convenient. Trivial question if you know off the top of your head - why <> rather than < ListboxSelect>? I never saw any examples of that type of enclosure. I cced the list. Ran into a lo

[R] question about how the mantel test is run in the ade4 package

Hi, I'm running a mantel test of the location of individuals in space against the behavioural score. I wanted to see whether those with similar behavioural scores are closer together in space relative to those that are further apart. This is the code I have written: nest.dists<-dist(cbind(aa$lon

Re: [R] plot in package psych with function error.bars.by

Hi Jim, Thank you for your answers. I also don't know where does the error like you mentioned come from. But for me, R works normally with the codes following: require(psych) keys.list=list(Agree=c(-1,2:5),Conscientious=c(6:8,-9,-10),Extraversion=c(-11,-12,13:15),Neuroticism=c(16:20),Openness =

Re: [R] Loop Autoregression

Hi I am not sure if Srates is a function I found in R-Forge. If yes, NSSParameter has to be time series as there is command time(x) in the function which, in case of data frame, gives this error. > time(sadr) Error in hasTsp(x) : invalid time series parameters specified Regards Petr > -Or

Re: [R] Error clmm(){ordinal}

Aurore, I don't know if car::Anova is able/should be able to produce anova tables for clmm objects; I usually use drop1() (and sometimes add1) to test terms in CLMMS: > library(ordinal) > fm1 <- clmm(rating ~ temp + contact + (1|judge), data=wine) > drop1(fm1, test="Chi") Single term deletions M

Re: [R] Error clmm(){ordinal}

Dear Rune, Thanks a lot for your comment! I tried with "Location" as a fixed-effect and clmm() ran! It take quiet a long time but I have approximatively 67000 observations in my datatable so I guess it's quiet normal! I want to run function Anova{car} on Modele but it gave my another error : >Ano

Re: [R] Loop Autoregression

Hi Jonas Actually I am a poor chemist and so I do not know what you are talkin' about. On stackoverflow I found that the error is probably due to wrong format of NSSParameter which probably is data frame and it shall be xts. Maybe others can help you more but I believe they will expect more inf