Displaying 20 results from an estimated 100 matches similar to: "write.table() error"
2013 Apr 13
2
Comparison of Date format
Hi,
?In the example you provided, it looks like the dates in Date2 happens first.? So, I changed it a bit.?
DataA<- read.table(text="
ID,Status,Date1,Date2 ??? ??? ??????
1,A,3-Feb-01,15-May-01 ??? ???
1,B,15-May-01,16-May-01 ??? ???
1,A,16-May-01,3-Sep-01 ??? ??? ??? ??? ???
1,B,3-Sep-01,13-Sep-01 ??? ??? ??? ??? ???
1,C,13-Sep-01,26-Feb-04 ??? ??? ??? ??? ???
2008 Sep 17
2
adding rows to table
Greetings everyone,
I'm trying to add a specific table or a specific number of rows (e.g.44) to a table with no success.
This is my basic table
> head(dataA)
year plot spp prop.B DCA1 DCA2 DCA3 DCA4
1 2000 1 a1 0.031079 -0.0776 -0.0009 0.0259 -0.0457
2 2000 1 a2 0.968921 -0.0448 0.1479 -0.1343 0.1670
3 2000 2 a1 0.029218
2009 Jun 09
2
is it possible to combine multiple barplots?
i am working with two sets of likert scale type (4 distinct values) data:
dataA <- rep(1:4, c(3,2,2,4))
dataB <- rep(1:4, c(5,4,3,2))
i can now (bar)plot both of these separately and compare the distributions.
plot(table(dataA), type='h')
plot(table(dataB), type='h')
is there a way to plot both of them in one plot, so that the bars for
value "1" (dataA: 3,
2010 Apr 30
1
gluster-volgen - syntax for mirroring/distributing across 6 nodes
NOTE: posted this to gluster-devel when I meant to post it to gluster-users
01 | 02 mirrored --|
03 | 04 mirrored --| distributed
05 | 06 mirrored --|
1) Would this command work for that?
glusterfs-volgen --name repstore1 --raid 1 clustr-01:/mnt/data01
clustr-02:/mnt/data01 --raid 1 clustr-03:/mnt/data01
clustr-04:/mnt/data01 --raid 1 clustr-05:/mnt/data01
clustr-06:/mnt/data01
So the
2013 Jan 18
2
A smart way to use "$" in data frame
Hello all,
I have a data frame dataa:
newdate newstate newid newbalance newaccounts
1 31DEC2001 AR 1 1170 61
2 31DEC2001 VA 2 4565 54
3 31DEC2001 WA 3 2726 35
4 31DEC2001 AR 3 2700 35
The following gives me the balance of state AR:
2009 Jul 20
1
a complicated merging task
I would like to merge two dataframes, but i have a condition that needs to
used for the merge as well.
the rows (observations) in each dataframe are identified by each person's ID
and by the date of the observation.
Basically I would like it to be merged based on both ID (exact match) and
date (a condition where one dataframe's date must be after the other
dataframe's date).
2009 Jun 18
1
simple question on glm
Hi,
I am trying to use glm to fit my data, wondering if there is a easy way to
fit a glm without typing all the explanatory variable names. For example, if
I have 100 explanatory variables x1, x2, ..., x100 and response variable is
y, I don't want to do something like
glm1 <- glm(y ~ x1 + x2 + ... + x100, family = gaussian, data = dataA)
since it would be a lot of typing.
Many thanks,
2005 Sep 22
1
R2WinBUGS: Data loading error
Hi R-Help!
I am trying to use R2WinBUGS but I get the following error message in WinBUGS
(and there must be something wrong with my R statement as I tried it directly in
WinBUGS and it worked):
display(log)
check(C:/Documents and Settings/Daikon/Roche/pop_model.txt)
model is syntactically correct
data(C:/Documents and Settings/Daikon/Roche/data.txt)
expected key word structure
compile(7)
...(and
2012 May 06
3
PLot a matrix
Hi,
I want to plot this matrix (I attach the data), it is suposed that each
column is a different time series.
If I do
g<-read.table("dataADF.txt", header=F)
and
plot(g[,1],type="l")
it plots the first column plot if I want in a unique graph each colums of
dataA, all in one. How should I proceed?There is a direct pre-defined code?
And If I wanted a plot by each
2006 Sep 13
2
recursive methods for concatenating sets of files
Hello,
I would like to read sets of files within a folder, perhaps using recursive
methods.
Right now, I rename the files before import.
It would be even better to do this without renaming files, without providing
explicit filenames, perhaps by importing files based on chronology,
and translating each filename into a header?
Please excuse my ignorance, and help cure my clunky programming
2012 Nov 06
1
Multinomial MCMCglmm
Thanks for your answers Stephen and Ben,
I hope I am posting on the correct list now.
I managed so far to run the multinomial model with random effect with the
following command:
MCMCglmm(fixed=cbind(Apsy,Mygl,Crle,Crru,Miag,empty) ~
habitat:trait,random=~idh(trait):mesh,family="multinomial12",
data=dataA,rcov=~trait:units)
(where multiple responses are different species,
Habitat
2009 Jul 08
3
matching each row
I have two dataframes, the first column of each dataframe is a unique id
number (the rest of the columns are data variables).
I would like to figure out how many times each id number appears in each
dataframe.
So far I can use:
length( match (dataframeA$unique.id[1], dataframeB$unique.id) )
but this only works on each row of dataframe A one-at-a-time.
I would like to do this for all of
2011 Nov 30
0
formula for calculating the survival probability for nomogram
Hi, I used Dr. Harrell's rms package to make a nomogram.
Below is my code for nomogram and calculate total points and probability *in
original data set* used for building the nomogram. *My question is how I get
the formula for calculating the survival probability for this nomogram. Then
I can use this formula to do validation by using other data set. *
f1 <- cph(Surv(retime,dfs) ~
2010 May 04
1
Posix warning : Access to ... is crossing device
I have a distributed/replicated setup with Glusterfs 3.0.2, that I'm
testing on 4 servers, each with access to /mnt/gluster (which consists
of all directories /mnt/data01 - data24) on each server. I'm using
configs I built from volgen, but every time I access a file (via an
'ls -l') for the first time, I get all of these messages in my logs on
each server:
[2010-05-04 10:50:30] W
2011 Oct 10
1
[LLVMdev] Adding fixups and relocations late in code generation
Jim,
Both the branch and the branch target are in the same function. This is known delta and should not matter where it is relocated because the delta will remain the same once it is a .o or later.
I just want to know how to tell the compiler that this expression is target offset - branch offset.
More generally, I want to know the rules for setting up expressions in general. I have no idea if
2002 Apr 15
8
Problem
Hello! If I have a matrix as 1 2
2 3
and I want to change the value 2 in 0, what can I do?
Thank you
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)
2012 Aug 13
4
dimnames in an array(I'll be grateful if this message will be passed to all list users)
Hello,
I'm hoping someone with a wide experience with R may be able to see
what the program is trying to tell me.
I've got an array:
y1=rnorm(41,0.2)
y2=rnorm(41,0.2)
y3=rbind(y1,y2)
data11<-array(0,c(41,2,2))
data11[,1,]=y3
data11[,2,]=y3
rownames(data11)<-rownames(data11, do.NULL = FALSE, prefix = "Obs.")
colnames=c("V","R")
2012 Aug 08
1
dimnames in array
Hello,
I'm working with an array; I'm trying to make it so that an array of
dim(42,2,2) has names whose length corresponds to that of the array,
and am hoping someone with experience with this can see what I'm not
doing correctly:
data11 = array(0,c(41,2,2))
y = lsoda(x0,times,fhn$fn.ode,pars)#This is make.fhn() from colloc
infer package#
y = y[,2:3]
2008 Mar 03
1
Formating a zoo dataset .
Suppose I have following dataset :
> head(data1)
Date Return
1 03/31/00 0.14230650
2 04/28/00 -0.03276228
3 05/31/00 -0.06527890
4 06/30/00 -0.04999873
5 07/31/00 -0.01447902
6 08/31/00 0.22265729
Now I convert it to zoo object :
> data11 = zoo(data1[,2], as.Date(data1[,1], format="%m/%d/%y"))
> head(data11)
2000-03-31 2000-04-28 2000-05-31
2011 Nov 30
1
Nomogram with stratified cph in rms package, how to get failure probability
Hello,
I am using Dr. Harrell's rms package to make a nomogram. I was able to make
a beautiful one. However, I want to change 5-year survival probability to
5-year failure probability.
I couldn?t get hazard rate from Hazard(f1) because I used cph for the model.
Here is my code:
library(rms)
f1 <- cph(Surv(retime,dfs) ~
age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11,