similar to: Interpretation of randomForest results

I have a multilevel dataframe (df): ID Date Segment Slice Tract Lesion 1 CSPP005 12/4/2007 1 1 LCST 0 2 CSPP005 12/4/2007 1 1 LPC 2 3 CSPP005 12/4/2007 1 1 RPC 3 4 CSPP005 12/4/2007 1 1 RCST 1 5 CSPP005 12/4/2007 1 1 LGM 0 6 CSPP005 12/4/2007 1 1

Dear R-helpers, I would appreicate if someone can help me on the transfer parameter in ARIMAX and also see what I am doing is correct. I am using ARIMAX with 2 Exogeneous Variables and 10 years data are as follows: DepVar Period, depVar, IndepVar1 Period, indepVar1, IndepVar2 Period, indepVar2 Jan 1998,708,Jan 1998,495,Jan 1998,245.490 Feb 1998,670,Feb 1998,421.25,Feb 1998,288.170 Mar

I am trying to call a C++ k-means function within R and I am struggling. I know that the below code is used to call a C++ function for gbm but how do I do it for k-means? gbm.obj <- .Call("gbm", Y=as.double(y), Offset=as.double(offset), X=as.double(x), X.order=as.integer(x.order),

I have a dataset of 10000 records which I want to use to compare two prediction models. I split the records into test dataset (size = ntest) and training dataset (size = ntrain). Then I run the two models. Now I want to shuffle the data and rerun the models. I want many shuffles. I know that the following command sample ((1:10000), ntrain) can pick ntrain numbers from 1 to 10000. Then I just

Dear list, I used glm to analyze the effect of treatment, measurement-methods etc. on lesion size in a time course experiments. Could anyone points me out how to extract detailed comparisons of interest? For example: Measurement.T2.M6 vs. Measurement.T2.Baseline, Or even Measurement.T2.M6.Treatment vs. Measurement.T2.M6.placebo (I did not include Treatment as a factor in the query below yet).

Hi I have data from 12 subjects. The measurement is log(expression) of a particular gene and can be assumed to be normally distributed. The 12 subjects are divided into the following groups: Infected, Vaccinated, Lesions - 3 measurements Infected, Vaccintaed, No Lesions - 2 measurements Infected, Not Vaccinated, Lesions - 4 measurements Uninfected, Not Vaccinated, No Lesions - 3 measurements

Dear All, I am far from being a guru about parallel programming. Most of the time, I rely or randomForest for data mining large datasets. I would like to give a try also to the gradient boosted methods in GBM, but I have a need for parallelization. I normally rely on gbm.fit for speed reasons, and I usually call it this way gbm_model <- gbm.fit(trainRF,prices_train, offset = NULL, misc =

Hi R User, I was trying to find a final model in the following example by using the Boosted regression trees (GBM). The program gives the fitted values but I wanted to calculate the fitted value by hand to understand in depth. Would you give moe some hints on what is the final model for this example? Thanks KG ------- The following script I used #----------------------- library(dismo)

Dear R-helpers. Would appreciate if someone can explain the warning messages below, after auto.arima. I couldn't find any clue in the archived help. Also, how do I retrieve the AICs of each tried model in auto.arima? The purposes are (1) to output to a text file, and (2) to find the 2nd best model by finding 2nd lowest AIC instead of eyeballing thru the value at the console

When I mount samba partitions using smbmount on a Linux box under /mnt/linnet_g and run this command: find /mnt/linnet_g/ -type f|xargs wc I get a number of errors which look like this: wc: /mnt/linnet_g/analysis.lesion/b19o61/b61o19.9/b61o18.9-042303.332.cbin: Input/output error What does this mean? The error is repeatable in a I am using samba-3.0.14a-6 on a Debian

Dear R users & Experts, This is just a curiousity, I was wondering why the dominant eigenvetor and eigenvalue of the following matrix is given as the third. I guess this could complicate automatic selection procedures. 0 0 0 0 0 5 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 Please

Hi all, Just migrated users to a samba 4 server on built on Solaris 11 samba pkg. Some users are mapping all required drives (G:, H:, and I:) but some are not able to map them. The affected users user log files indicate that the %g variable is being translated to '-1' instead of the users Unix group. I would welcome any feedback or suggestions on how to resolve this issue. >From my

Hi all, Samba 2.2.8a is currently running on Solaris 8 as a member of a domain. For testing purposes I have installed Samba 3.0.10 on a test server (Solaris 9) a a member of a Workgroup. In a couple of weeks we will be upgrading our production server to Solaris 9 and are planning on doing a fresh install of Samba 3.0.10 at the same time. I will be saving the following files to be restored after

I am in the process of planning for a samba upgrade. The current version of Samba (2.2.8a) runs on a Sunfire 280R with Solaris 8. The intended new platform is a sunfire V240 with Solaris 9 as the OS. The new version of Samba will be Samba-3.0.10. Please note that the Samba server will be a member of a domain. I've read Samba-3 by Example, Chapter 8 which describes to the process of updating

another issues with cisco 7975 I have phone registered on asterisk have 2 different issues on different versions of firmware, on 9-4-2-1S I have not working 3way conference, when I trying to connect second call, phone says ?unable to set up conference? and sending some cisco xml data to asterisk which cannot be handled, thats the problem, I know on firmware 8-5-4 3way conference works just

I have a system that is running into problems with high load, however, even when I shut down everything that could be loading the system, the load continues to rise. I can't see anything that is causing the load. Any ideas? 08:36:06 up 1 day, 10:57, 3 users, load average: 19.98, 19.30, 17.29 90 processes: 75 sleeping, 2 running, 13 zombie, 0 stopped CPU states: cpu user nice

I have several datasets that come from different studies (fv02 and fv03), they represent different levels (patients and lesions), and they have different patient populations (itt, mitt, mitt3). I wanted to write some code that would pass my three requirements into a function I wrote, produce the output, but not have to require me to also pass a unique plot title or output filename for each

I can't figure out why this function wont work. #Custom Function used fallInBand <- function(x){ #Assume the bands fall go up form 0 to 100 in intervals of 5 then from 100 to 1000 in intervals of 100. #returns the location (band number) if (is.na(x) == FALSE) { if(x < 100) {#bands of 5 if((x %% 5) >=0){#falls within the band Result = as.integer(x/5)+1

Hi Everyone, I made some small changes in my dovecot setup to switch it from looking up users and passwords from a mix of ldap (i.e. freeipa) and password files. One of the changes was to switch from using one id for all authentication to using individual ids) It's working fine with Evolution. I have one account authenticating with GSSAPI, which is my userid for logging into my desktop and

If you used MySQL FULLTEXT indexes and just upgraded to Rails 1.1, watch out for this one. Rails 1.1 now defaults to Ruby as your schema format. So instead of development-structure.sql, Rake will use schema.rb. That''s portable and good stuff overall but there is a slight compatibility issue. Ruby schema doesn''t know how to do FULLTEXT index syntax for MySQL. To work around it,