similar to: How to stack data sets?

Hi All, I'm generating 10 different data sets with 1 and 0 in a matrix form and writing the output in separate files. Now I need to stack all these data sets in one vector and I know that stack only operates on list or data frame however I got these data sets by converting list to a matrix so can't go backwards now. Is there a way i can still use Stack? Please see the program:

Hello All, I'm trying to print specific row and column for Observed_Scores matrix, however, when I execute the command "Observed_Scores[1,1]", I get the message "Error in Observed_Scores[1, 1] : incorrect number of dimensions". Could you please help and let me know where is the mistake? Here is my program: library(ltm) library(psych) # Settting the working directory

I am trying to simulate a dataset using Parallel Latent CTT model and this is what i have done so far: (START) #Importing psych library for all the simulation related functions library(psych) # Settting the working directory path to C:/NCME path="C:/NCME" setwd(path) #Using the function to generate the data GenData <- congeneric.sim(N=500, loads =

Hi, I am a new user of R 2.8.1. I use Tinn-R for code editing. I use a windows 2003 system with 1 GB RAM. I am interested to generate dashboard and reports based on data from MS Access. These reports need to be posted on a weekly basis to the web. The reporting interface should provide facilities for "what if" scenarios. Is it possible to interface R analysis results to good open

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

Hi First off, I apologise if this is covered somewhere and I just don''t get it... I''ve looked far and wide .... As a part of my current project I need to hook up to a particular XMLRPC webservice that won''t be available for me to use until close to launch date... So I am trying to recreate the webservice that mirrors the live one so I can test my application

Version 2.2-0 of the rms package is now available. This is a somewhat major update. One major change is not downward compatible: Instead of specifying predictor=. or predictor=NA to Predict, summary, nomogram, survplot, gendata, you just specify the name of the predictor. For example, to get predictions for the default range of x1 and for just 2 values of x2 you might specify Predict(fit,

Version 2.2-0 of the rms package is now available. This is a somewhat major update. One major change is not downward compatible: Instead of specifying predictor=. or predictor=NA to Predict, summary, nomogram, survplot, gendata, you just specify the name of the predictor. For example, to get predictions for the default range of x1 and for just 2 values of x2 you might specify Predict(fit,

Hi... I really need some help on something I am working on... I am trying to create an XML-RPC webservice that sends the following as a response... As you can see it looks like I can use a Struct with a boolean: member but I don''t seem to be able to have another Struct as a member? Am I barking up the wrong tree trying to do it this way? Should I try something else? I would really

I am fitting 3 parameter model to my response matrix and want to generate item characterstic curve. I want to specify file name to save item characterstic curve by passing it as external parameter to the R batch script. The following is the code I have written for this. *R Script:* library(ltm) cmd_args = commandArgs(); for (arg in cmd_args) cat(" ", arg, "\n",

Dear All, I am trying to simulate data for a spline/piecewise regression model. I am missing something fundamental in my simulation procedure because when I try to fit my simulated data using the Gauss-Newton method in SAS, I am getting some wacky parameter estimates. Can anyone please check my simulation code and tell me what mistake I am making in generating data for spline model? Thank you

If the observations are normally distributed and the 2xk design is balanced, theory requires that the tests for interaction and row effects be independent. In my program, appended below, this would translate to cntT (approx)= cntR*cntI/N if all R routines were functioning correctly. They aren't. sim2=function(size,N,p){ cntR=0 cntC=0 cntI=0 cntT=0 cntP=0 for(i in 1:N){

I would like to simulate a stochastic process by specifying an intensity function, say lambda(t)=0.05 + 3 exp(鈭抰) I have the following code, but it has a number of problems. One issue I have is I do not know where the time variable is. I would appreciate your help with this code. library("survival") n <- 5000 lamc<-0.08 k1<-0.05 k2<-3 ld<-function(t){k1+k2}

Greetings useRs, What is the easiest way to create a design matrix of several factor variables? Function gendata in Design seems to do that for a fitted model, but how to do that only on several factor vectors?? The result should be a df with one row for each distinct combination of levels of factors eg for (M,F) (Y,O) We get M Y M O F Y F O In reality I will have more than 1000 rows so doing

Trying to use contrast to look at differences within an lme lme.fnl.REML <- lme(Max ~ S + Tr + Yr + Tr:Yr, random = ~1 |TID, method = "REML") I have three levels of Tr I'm trying to contrast among different years (R, T97, T98), years = 1997-1999, so I'm interested in contrasts of the interaction term. > anova(lme.fnl.REML) numDF denDF F-value

I am cleaning up the rms package to not export functions not to be called directly by users. rms uses generic functions defined in other packages. For example there is a latex method in the Hmisc package, and rms has a latex method for objects of class "anova.rms" so there are anova.rms and latex.anova.rms functions in rms. I use:

Dear R-users I'm modelling some longitudinal data (1 response variable measured at 6 occasions, 1 baseline, one "treatment" variable) collected in the same subjects using the following model: library(nlme) model.lme <- lme(response ~ V0+ time+ tt + tt:time, random = ~1|subject, correlation = corSymm(form = ~ 1 | subject), na.action=na.omit, data=arm.total2) Then I would lke to

Hi! I try to explain the efffect of (1) forest where i took samples's soils (* Lugar*: categorical variable with three levels), (2) nitrogen addition treatments (*Tra*: categorical variable with two levels) on total carbon concentration's soil samples (*C: *continue* *variable) during four months of sampling (*Time:* categorical and ordered variable with four levels). I fitted the

> On Nov 1, 2017, at 12:51 PM, Tiby Kantrowitz <tlkantro at gmail.com> wrote: > > The prob package has been archived because it depends upon some other > packages which have issues. > > However, such projects as Introduction to Probability and Statistics in R > depend upon it for learning. There are a few other resources that also use > it. > > Does anyone

Hi. First, i would like to thank samba developers for producing such a good product. Second, i have a few questions/remarks : I have recently replaced a w2k file server running in w2k domain (native mode) with samba 2.2.7a on RH 7.3 with the latest kernel, no acl, configured winbind, and ran into the problem described here :