similar to: problem of data manipulation

Hi, My problem maybe a little bit complicated, so forgive me if the following words are too much. #date set a0<-matrix(c(1.1,1.3,1.1,1.3,1.3, 2.0,1.8,2.0,1.8,1.8, "12/01/2008","05/20/2007","12/06/2008","05/10/2007","05/06/2007", "N","N","A","C","A", 1,2,3,4,5),ncol=5,byrow=FALSE)

Hello. I am probably wrong... I am wondering if it could have a mistake in the code of the ansari.test function. For me, it seems that the function do not recover the p value at the correct side of the normal law N(0, 1) when it use the normal approximation (presence of ties) in a one tailed test. Exemple : quanti<-c(197, 205, 228, 234, 237, 195, 233, 226, 244, 227, 259, 185, 198, 253, 207,

Hi, I have Forecast Class and Observed Class in a data matrix as below. > Sample1 FCT OBS 1 1 5 2 2 4 3 3- 3+ 4 3 3 5 3+ 3- 6 4 2 7 5 1 I want to find the difference between Observed and Forecast Classes. How can I get this done? I tried to following to convert the 1 through 5 classes, to 1 through 7 for both OBS and FCT column. > Sample1$OBS2 <- Sample1$OBS

Let's say I have a within-subject experiment with 2 observables, obs1 and ob2 and 2 independent factors, fac1 and fac2. I can do summary( aov( obs1~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) summary( aov( obs2~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) to test the 2 observables separately. > summary( fit<-manova( cbind(obs1,obs2)~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) gives

hi, all I inherited a data set with format like the following: para.A para.B para.C observation.1 observation.2 observation.3 .... Essentially, the same experiment (with same parameters) are conducted 3 times, with corresponding observations. What are the efficient ways of re-arranging this data set in a one row per observation such that: para.A para.B para.C observation.1 para.A

Hi, newbie question. I have a data-frame with 3 named columns: Name, Obs1, Obs2. The Name column members are made of alphanumeric characters: T1, T2, T3 etc. I would like to acess only that subset of the data-frame with Name == T44. X <- dataframe[dataframe$Name=='T44'] does not work. Any ideas on how to do this? I'm sure I'm missing a simple concept here. Thanks, Anjan --

Hello, I want to do fisher test for the rows in data file which has value less than 5 otherwise chi square test .The p values from both test should be stored in one resulted file. but there is some problem with bold if statement. I don't know how implement this line properly. x = cbind(obs1,obs2,exp1,exp2) a = matrix(c(0,0,0,0), ncol=2, byrow =TRUE) #matrix with initialized values

Hi, To simplify, suppose I have 2 observations each day for three days. I would like to define the correlation structure of these 6 observations as follows: the correlation of 2 observations on the same day is, say, alpha, the correlation for 2 observations one day apart is rho and the correlation for 2 observations 2 days apart is rho^2. I.e. I would like to have an AR1 correlation + a

Dear R Users, I have some data that is very similar in form to a 3D image - ie univariate data on a regular 3D grid. I keep this as a 3D numeric array in R with attributes describing the sampling points along the 3 dimensions. I would like to interpolate this onto a new regular 3D grid that I specify (eg by supplying 3 vectors corresponding to the new grid locations on each of the 3

Hi, has anybody managed to compile samba-2.0.6 with the sfio-library? I get the following error: Using FLAGS = -O -I./sfio/include -Iinclude -I./include -I./ubiqx -I./smbwrapper -D_LARGEFILE_SOURCE -D_FILE_OFFSET_BITS=64 -DSMBLOGFILE="/var/samba/log.smb" -DNMBLOGFILE="/var/samba/log.nmb" -DCONFIGFILE="/usr/local/samba/lib/smb.conf"

Full_Name: Martha Nason Version: 2.0.1 OS: Windows XP Submission from: (NULL) (137.187.154.154) I am running simulations using fisher's test on 2 x c tables and a very small p.value from fisher's test (<2.2e-16) is returned as a negative number. Code follows. > set.seed(0) > nreps.outer <-7 > pvalue.fisher <- rep(NA,nreps.outer) > > population1 <- c(

Hi! Well, I have a Linux machine (RedHat 9.0) with samba on it and a Win ME on the other machine and I want to share a partition for music, read-only. (Won't be hard, is mounted read-only) and a transfer share where I can write to, too. My config file looks as follwing: # Samba config file created using SWAT # from 192.168.0.1 (192.168.0.1) # Date: 2003/08/04 23:05:42 # Global parameters

I apologize in advance for the simplicity of this question. I use R 2-3 times a year, and I seem to forget more in the intervening months than I learn during my days of panicked reading.... I HAVE tried looking at the help resources; I'm just not very good at understanding them. I have a dataframe with 18 observations of 5 different things, and a second dataframe with and estimate for those

Thank you very much for your reply. Your code work well with this example. I modified a little to fit my real data, I got an error massage. Error in split.default(x = seq_len(nrow(x)), f = f, drop = drop, ...) : Group length is 0 but data length > 0 On Thu, Jan 31, 2013 at 12:21 PM, arun kirshna [via R] < ml-node+s789695n4657196h87@n4.nabble.com> wrote: > Hi, > Try this: >

Hi to all! I have a problem with a nvidia geforce 520mx [NVd0 generation card (0x0d9110a1)] i have on my notebook (samsung 3 series). In practice i'm not able to use it with bumblebee and bbswitch. The dmesg message is: [ 13.507435] nouveau 0000:01:00.0: power state changed by ACPI to D0 [ 13.507440] nouveau 0000:01:00.0: power state changed by ACPI to D0 [ 13.507448] nouveau

Hello ! I am attempting to switch from being a long time SAS user to R, and would really appreciate a bit of help ! The first thing I do in getting a large dataset (thousands of obervations and hundreds of variables) is to run a SAS command PROC CONTENTS VARNUM command - this provides me a table with the name of each variable, its type and length; then I run a PROC MEANS - for numerical

hello! my question conserns with use of "panel" package (written by R.C.Gentlman) (unfortunately the manual and help sites are very short) 1. is it possible to do analysis just without a(ny) covariate? i suggest do it by introducing a covariate with level=0 in all obervations, this because of Q(z)=Q_o exp(beta*z), but it seemingly doesn't work 2. in the option gamma in the

Colleagues, I am trying to de-duplicate a large (long) database (approx 1mil records) of diagnostic tests. Individuals in the database can have up-to 25 observations, but most will have only one. IDs for de-duplication (names, sex, lab number...) are patchy. In a first step, I am using Andreas Borg's excellent record linkage package (), that leaves me with a list of 'pairs' looking

Hi, dat1<- read.table(text=" id??????????????? t???????????????????? scores 2???????????????? 0??????????????????????? 1.2 2???????????????? 2???????????????????????? 2.3 2???????????????? 3??????????????????????? 3.6 2???????????????? 4??????????????????????? 5.6 2???????????????? 6??????????????????????? 7.8 3???????????????? 0??????????????????????? 1.6 3????????????????

Hi there, I am trying to interface c++ code in R and make a package. With R CMD SHLIB the dll was created, but when I try R CMD check, I am getting 'undefined reference to..' linkage error messages. The relevant c++ source from conf-infomap.cpp: #include "conf-infomap.h" #include "R.h" // R functions #include "Rinternals.h" #include "Rmath.h" //