Thanks for the help! But when I tried it, it does not work the same way I
want. :(
after combining the two matrices, they look like this:
V1V2 X TEL.AML1.C41
Hyperdip.50.C23
1 TEL.AML1.C41 1TEL.AML1.C41 1.000 0.
I have used tm package to import a set of text documents using the
following command:
text - Corpus(DirSource(.),readerControl = list(language =ansi))
I would like to extract only a certain portion of the text in each document
using certain keywords. For example, I would like to include all the
On 26.06.2011 21:26, statfan wrote:
I am trying generate a sample for a truncated multivariate normal
distribution via the rtmvnorm function in the {tmvtnorm} package.
Why does the following produce NaNs?
rtmvnorm(1, mean = rep(0, 2), matrix(c(0.06906084, -0.07463565, -0.07463565,
On 26.06.2011 13:36, scaramouch wrote:
Hi everyone.
My problem: I've 4 distributions: A, B, C and D:
#A
mA=16.6
sA=3.0
#B
mB=18.9
sB=3.2
#C
mC=20.3
sC=2.1
#D
mD=24
sD=0.8
###Graphiques ensembles
plot(function(x) dnorm(x,mA,sA),0,40,col=orange,ylim=c(0,0.5))
plot(function(x)
Dear Firas,
Regression on the subdistribution hazard can be performed by fitting a
weighted Cox model. See Geskus, Biometrics 67, p. 39-49, 2011. Hence, cph
(and coxph) can be used directly; there is no need to use the crr function
in cmprsk. The result from the weighted cph fit should allow you
On Sun, Jun 26, 2011 at 06:34:28PM -0700, Ungku Akashah wrote:
hello. I need some help about this R software. I've been searching
for volcano plot(statistic) script for long, but still not found.
May i request the script for volcano plot. If able, pls include any
tips about volcano plot.
Hmmm. Caution would seem more appropriate than fear.
My point was that if you want the results as numbers rather than characters,
it seems perverse to convert them to characters with print; you should use a
suitable mechanism to return them as numbers. If you're nervous of returning
values via
Hi
r-help-boun...@r-project.org napsal dne 24.06.2011 16:51:27:
Thanks for your reply, but that is not quite what I am looking for...I
do not
want to add up all the values in the row, I want to get the number of
columns in each row that meet the criteria and then get that as a
In your dummy variable trick, I see how subjectwise predictions are
obtained. But to get the Group mean predictions, isn't the s term
with the 0 dummy by variable the same as just omitting the s term
with Subject altogether? I would think so but want to check.
- Yes it's the same as omitting
Good morning,
I was using the packages snow 0.3-3 and snowfall 1.84 with R.10.1 : it is
working great ! (Windows 7)
Since I have tried to update my R version (Windows 7). I cannot make the
packages snow and snowfall works for parallel processing. I have tried the
following version
Please try R-2.13.0 patched with the most recent snow. We found a bug in
R that has been fixed in the meantime.
Best,
Uwe Ligges
On 27.06.2011 10:46, Martin Valere wrote:
Good morning,
I was using the packages snow 0.3-3 and snowfall 1.84 with R.10.1 : it is
working great ! (Windows 7)
Hallo!
Can anybody help with the comparison of ROC curves? I want to compare two
uncorrelated ROC curves, i.e. the data sets of the two ROC curves are from
different samples.
One way to compare ROC curves seems to be the comparison of the AUCs. Although
I
found the package pROC with he
Dear All,
I am trying to use the glmmBUGS package with OpenBUGS under Linux.
Here is a short list of code I modified from the diseasemapping
package document:
library(diseasemapping)
data(popdata)
data(casedata)
model = getRates(casedata, popdata,
On 06/26/2011 07:26 PM, wjcao wrote:
Dear R members
Could someone tell me how to simulate Likert-type data using the rnorm
function.
Let's say, 200*15 random numbers in a variable that goes from 1 to 4 in
steps of
1 (i.e., 1, 2, 3, 4) belonging to a normal distribution?
random.data—
Addition of a cluster() term fits a Generalized Estimating Equations
(GEE) type of model, addition of frailty() fits a random effects model
(Mixed Effect or ME). In glm analysis (linear regression, logistic
regression, etc) the arguments about the advantages/disadvantages of GEE
ve ME would
Thanks! I did not realize you can access variables by name like this.
Serguei
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide
Hi,
I'm mighty new to R. I'm using it on Windows. I'm trying to cluster using a
distance matrix I created from the data on my own and called it D10.dist. I
loaded the cluster package. Then tried the following command...
agnes(E:D10.dist, diss = TRUE, metric = euclidean, stand = FALSE,
method =
Hi dear all,
I am looking for a R algoritm which takes logistic relative risk function
versus exponential relative risk function in cox proportional hazard
model.
Let me know how I can do it?
Thank you very much!
Rasouli
Rasouli
--
This message has been scanned
Is there any easy way I can upgrade R without losing all of my packages, when
I'm using a computer (Mac) on an organisation's network which blocks the
package installer?
Thank you,
Thomas
--
View this message in context:
http://r.789695.n4.nabble.com/Upgrading-R-tp3627205p3627205.html
Sent
I will try it out..Thanks a lot
On Mon, Jun 27, 2011 at 10:08 AM, rgeskus [via R]
ml-node+3627184-664831548-247...@n4.nabble.com wrote:
Dear Firas,
Regression on the subdistribution hazard can be performed by fitting a
weighted Cox model. See Geskus, Biometrics 67, p. 39-49, 2011. Hence, cph
Dear sir,
If I have a vector as:
y - c(1:4)
and a matrix as:
x - matrix (5:12, nrow=4, ncol=2)
Then I create a data frame as:
t - data.frame(y, x)
R will generate a data frame of 3 columns
names(t)
y X1 X2
is.vector(t$X1) returns TRUE
is.vector(t$X2) returns TRUE
R splits the columns of a
Hi, guys:
I try to plot a time series against another to see the relationship between
them. However, both plot.ts and ts.plot sort the independent series values
first then do plotting, i.e. it changes time order. Could anyone help with
this please?
Appreciate
C MEN
--
View this message in
Thanks a lot. That solves the problem :)
best,
salih
On Sat, Jun 25, 2011 at 9:13 AM, Dimitris Rizopoulos
d.rizopou...@erasmusmc.nl wrote:
yet another solution is:
v - c(1, 2, 6, 5)
w - matrix(c(1, 4, 3, 5, 8, 10), ncol = 2, byrow = TRUE)
check - outer(w[, 1], v, -) * outer(w[, 2], v, -)
I wondered if there were standard practices in CRAN for delivery of R source
implementing functions in R packages. I has encountered a couple of packages
where the gzipped version of source contains very little, primarily the Help
files describing the functions in the package. In some cases I
Thank you Uwe !
R-2.13.0 patched did indeed solve my issue and the package snow is
working again fine ...
Best,
Valère
-Ursprüngliche Nachricht-
Von: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de]
Gesendet: Montag, 27. Juni 2011 11:47
An: Martin Valere
Cc:
hi!
I apologize in advance if this is a newbie dumm question, but I really can't
figure it ou.
I have lists of sumeric and character data on some URLs, which look like
this:
photo id=5876248819 owner=13716719@N04 secret=faf9bb7f52
server=5264 farm=6 title=our rose garden ispublic=1 isfriend=0
Hi
First a handy point : When supplying sample data it is a good idea to use
dput(). See ?dput for an explanation. It makes it much easier to see what the
data looks like and to work with it. Sample data pasted into a e-mail can get
badly mangled. Yours was not bad but still need a bit of
On Mon, Jun 27, 2011 at 1:24 PM, Galkowski, Jan jgalk...@akamai.com wrote:
I wondered if there were standard practices in CRAN for delivery of R
source implementing functions in R packages. I has encountered a couple of
packages where the gzipped version of source contains very little,
On Jun 27, 2011, at 7:24 AM, Galkowski, Jan wrote:
I wondered if there were standard practices in CRAN for delivery of R source
implementing functions in R packages. I has encountered a couple of packages
where the gzipped version of source contains very little, primarily the Help
files
It's often difficult to figure out what a problem data set looks like when
someone pastes a mess of output from R into an e-mail. The dput command, used
judiciously—I would hope no one would send a 1 M data set—can make life much
easier.
I wonder if it would be useful to add a hint in the
hi,
because of my ignorance I can't figure out why would you need such a structure.
However the following trick works:
a - matrix(1:6,3)
b - data.frame(1:3)
b$a - a
b
X1.3 a.1 a.2
11 1 4
22 2 5
33 3 6
names(b)
[1] X1.3 a
class(b$a)
[1] matrix
is.matrix(b$a)
[1]
On 2011-06-27 03:02, 新鼎-智慧製造事業處-鄭紹文 wrote:
Dear sir,
If I have a vector as:
y- c(1:4)
and a matrix as:
x- matrix (5:12, nrow=4, ncol=2)
Then I create a data frame as:
t- data.frame(y, x)
R will generate a data frame of 3 columns
names(t)
y X1 X2
is.vector(t$X1) returns TRUE
Dear UseRs,
I built varying coefficient models (in mgcv) for two groups separately, with
one explanatory and one moderator variable (see the example below).
# ---
# Example:
# --
# generate moderator variable (can the same for both groups)
modvar - c(1:1000)
# generate group1
I also tried...
D1 D2 D3 D4 D5 D6 D7 D8 D9
D2 0.608392
D3 0.4974510.537662
D4 0.6345480.3933430.537426
D5 0.5587850.5433990.6322210.726633
D6 0.6594830.701778
I actually thought it was in the posting guide, but I could not find it
with a text search. Did I overlook it?
I'm not sure that adding the dput() function at the bottom of every
e-mail will help, since you don't see it when you write an e-mail, just
when you receive one. Though it cannot
On Mon, Jun 27, 2011 at 3:43 AM, Karen R. Khar karen.k...@gmail.com wrote:
Hi,
I'm mighty new to R. I'm using it on Windows. I'm trying to cluster using a
distance matrix I created from the data on my own and called it D10.dist. I
loaded the cluster package. Then tried the following
On Jun 27, 2011, at 9:44 AM, Ivan Calandra wrote:
I actually thought it was in the posting guide, but I could not find
it with a text search. Did I overlook it?
The Posting Guide has an equivalent method it suggests. It has the
trivial advantantage of being truly cut-and-pastable but the
On Mon, Jun 27, 2011 at 1:24 PM, Galkowski, Jan jgalk...@akamai.com wrote:
I wondered if there were standard practices in CRAN for delivery of R source
implementing functions in R packages. I has encountered a couple of packages
where the gzipped version of source contains very little,
Hi all,
I'm currently working on updating an R package to run for loops in parallel
to speed up computation time. I'm using the 'foreach' package with a
foreach loop. When I run my code inside the loop, I get an error message
that a number of the functions aren't recognized (even though the
xqzhang85 at yahoo.com writes:
Hi, I can not estimate the zero inflated mixed model parameters successfully.
If it is possible, would you
please help me write the code? In my case, I consider the location as a random
effect. My data are:
location y x1 x2 x3
1
1
1
2
2
2
3
3
3
Hey all,
I learned that using the equals sign = to assign objects is generally OK, but
will not work in some cases.
As I always use - for assignments, I have not encoutered any problems.
Could somebody provide an example or explanation, why getting used to = is
not a good idea?
Or is it?
I'm not an expert but some third party libraries won't recognize = as an
assignment for some objects (I learned this empirically) when the
assignment is made within a function. Using - protects you from these
unknows, as you have discovered on your own.
Regards
David
On 6/27/2011 8:34 AM,
Hey again,
I didn't wat questions to get mangled up, so here's my second email.
In matlab, there is the simple possibility to access colums x till last of a
matrix using
mydata(1:3, 5:end).
In R, I so far use
mydata[1:3, 5:ncol(mydata)]
Is there a faster way? (in terms of typing)
Thanks
Hi, I want to fit a GEE with a negative binomial distribution. I have uesd
already a poisson glm and then neg binommial to deal with alot of
dispersion. In my neg binomial residuals i have some patterns so i have
implemented a GEE, but only with a poisson family as i couldnt with neg
binomial.
Fine. Attached. It's waved.
All it has is *.Rd files. Apparently the functions are collected in
functionINIT.R. But 00Index and DESCRIPTION are not helpful.
- j
-Original Message-
From: b.rowling...@googlemail.com [mailto:b.rowling...@googlemail.com] On
Behalf Of Barry Rowlingson
Apologies David,
Every once in a while I hit the wrong reply button. The post was meant to be a
Reply All. I'll post this one to the list.
--- On Mon, 6/27/11, David Winsemius dwinsem...@comcast.net wrote:
From: David Winsemius dwinsem...@comcast.net
Subject: Re: [R] Add use of dput() to
On 27.06.2011 17:14, Galkowski, Jan wrote:
Fine. Attached. It's waved.
All it has is *.Rd files. Apparently the functions are collected in functionINIT.R. But
00Index and DESCRIPTION are not helpful.
Are you joking or just pulling us?
All these files are interpreted by R and essential for
On Mon, Jun 27, 2011 at 11:14 AM, Galkowski, Jan jgalk...@akamai.com wrote:
Fine. Attached. It's waved.
All it has is *.Rd files. Apparently the functions are collected in
functionINIT.R. But 00Index and DESCRIPTION are not helpful.
All the R code is in functionINIT.R. Are you objecting
On Mon, 2011-06-27 at 11:14 -0400, Galkowski, Jan wrote:
Fine. Attached. It's waved.
All it has is *.Rd files. Apparently the functions are collected in
functionINIT.R. But 00Index and DESCRIPTION are not helpful.
- j
The Rd files are the help or manual pages for the functions defined
On 06/27/2011 10:34 AM, Berry Boessenkool wrote:
Hey all,
I learned that using the equals sign = to assign objects is generally OK, but
will not work in some cases.
As I always use - for assignments, I have not encoutered any problems.
Could somebody provide an example or explanation, why
Hi,
mydata[1:3, -(1:4)]
saves some typing.
cheers
Am 27.06.2011 16:51, schrieb Berry Boessenkool:
Hey again,
I didn't wat questions to get mangled up, so here's my second email.
In matlab, there is the simple possibility to access colums x till last of a
matrix using
mydata(1:3,
No, you are correct. It meets the letter of GPL. Took me a while to find
FunctionINIT.R though. As I wrote in the original, source is available, if only
by keying in function names and seeing their value. Was hoping for greater
clarity.
I of course have the paper and the help. I was trying to
On Mon, Jun 27, 2011 at 10:34 AM, Berry Boessenkool
berryboessenk...@hotmail.com wrote:
Hey all,
I learned that using the equals sign = to assign objects is generally OK,
but will not work in some cases.
As I always use - for assignments, I have not encoutered any problems.
Could
Yes.
Perhaps also worth mentioning in this context because they're generic
are head() and tail() -- ?head
Cheers,
Bert
On Mon, Jun 27, 2011 at 9:01 AM, Eik Vettorazzi
e.vettora...@uke.uni-hamburg.de wrote:
Hi,
mydata[1:3, -(1:4)]
saves some typing.
cheers
Am 27.06.2011 16:51, schrieb
G'day Gavin,
On Mon, 27 Jun 2011 16:36:57 +0100
Gavin Simpson gavin.simp...@ucl.ac.uk wrote:
[...]
I don't recall the GPL mentioning anything requiring that the source
code be helpful. The authors have most certainly fulfilled their
requirements under GPL, as has CRAN in distributing the
Look at package smatr. It fits RMA (called here standardized major axis)
regression, plots the line, and provides confidence limits.
--
David L Carlson
Associate Professor of Anthropology
Texas AM University
College Station, TX 77843-4352
Look at ?approx. For your example (of course your random numbers give
different results):
approx(f$x, f$y, c(-2, -1, 0, 1, 2))
$x
[1] -2 -1 0 1 2
$y
[1] 0.03757113 0.19007982 0.31941779 0.37066592 0.10227509
approx gives NA's if you try to interpolate outside the bounds of the data.
HI, Dear R-community,
I have one vector of the following string,
test-c(H.1.Y, N.0.E, G.1.N, E.0.P, W.2.G, W.4.G, W.3.W,
W.0.N, D.1.H)
test
[1] H.1.Y N.0.E G.1.N E.0.P W.2.G W.4.G W.3.W W.0.N D.1.H
*I want to change them into*
H{1}Y N{0}E G{1}N W{2}G E{0}P W{4}G W{3}W W{0}N D{1}H
I tried
Is this what you want:
test
[1] H.1.Y N.0.E G.1.N E.0.P W.2.G W.4.G W.3.W W.0.N D.1.H
sub(\\.([0-9]+)\\., {\\1}, test)
[1] H{1}Y N{0}E G{1}N E{0}P W{2}G W{4}G W{3}W W{0}N D{1}H
On Mon, Jun 27, 2011 at 1:18 PM, Changbin Du changb...@gmail.com wrote:
HI, Dear R-community,
I have one vector
Thanks, Jim!
It is exactly what I want, appreciated!
On Mon, Jun 27, 2011 at 10:22 AM, jim holtman jholt...@gmail.com wrote:
Is this what you want:
test
[1] H.1.Y N.0.E G.1.N E.0.P W.2.G W.4.G W.3.W W.0.N D.1.H
sub(\\.([0-9]+)\\., {\\1}, test)
[1] H{1}Y N{0}E G{1}N E{0}P W{2}G W{4}G
Hi,
I have the following data management issue. I am trying to combine multiple
years of ethnicity data into one variable called ethnic. The data looks
similar to the following
idethnic07ethnic08 ethnic09ethnic10
1 1 1 11
2
On Jun 27, 2011, at 12:56 PM, Christopher Desjardins wrote:
Hi,
I have the following data management issue. I am trying to combine
multiple
years of ethnicity data into one variable called ethnic. The data
looks
similar to the following
idethnic07ethnic08 ethnic09ethnic10
1
Regarding the subject, I want to thank the many respondents for clarifying the
nature of the relationship between R and the GPL, as well as giving help with
the structure of R-delivered source.
I want to emphasize I meant nothing at all harsh or accusatory in my email. I
did say I had access
Hello,
I am trying to simulate 10 relicates of 100-tables. Each table is a 2 x 3
and 80% pf the tables are true nulls and 20% are non-nulls. The nulls follow
the Hardy Weinberg distribution (ratio) 1:2:1.
I have the code below but the p-values are not what I am expecting. I want
to use the
SamiC s.cox.10 at aberdeen.ac.uk writes:
Hi, I want to fit a GEE with a negative binomial distribution. I have uesd
already a poisson glm and then neg binommial to deal with alot of
dispersion. In my neg binomial residuals i have some patterns so i have
implemented a GEE, but only with a
Hello,
thanks for your response.
Of course I already thought about this simple solution of the problem.
But I think this is not a nice workaround.
If I understand the manual correctly, density() already makes an
approximation. But this approximation is just evaluated at the wrong points.
See:
Thanks,
I have been using zuurs book but it only goes as far as poisson and binomial
GEE's. Initially I fitted a glm with poisson and this was over dispersed.
Then moved to binomial, but residual patterns are not great (ie variance).
Looks like some spatial correlation. However, in the GEE
On Jun 27, 2011, at 3:18 PM, Carsten Harlaß wrote:
Hello,
thanks for your response.
Of course I already thought about this simple solution of the
problem.
But I think this is not a nice workaround.
If I understand the manual correctly, density() already makes an
approximation. But this
Thank you for the answer, but how can i integrate the difference?
There is always an error with integrate().
Kindly regards, Scaramouch
--
View this message in context:
http://r.789695.n4.nabble.com/Overlapping-areas-under-normal-distributions-tp3625807p3628473.html
Sent from the R help mailing
Any ideas of what is happening when I attempt to load XLConnect? I think I
installed everything correctly. I posted what I get when I installed rJava
as well as when I attempted to load rJava directly. Thanks in advance for
any help you can give on this.
library(XLConnect)
Loading required
On Mon, 27-Jun-2011 at 01:22AM -0700, thomas.chesney wrote:
| Is there any easy way I can upgrade R without losing all of my
| packages, when I'm using a computer (Mac) on an organisation's
| network which blocks the package installer?
Keep a separate library for the packages you specifically
Hi, I am wondering if anyone can explain to me if cumulative incidence (CI) is
just 1 minus kaplan-Meier survival? Under what circumstance, you should use
cumulative incidence vs KM survival? If the relationship is just CI =
1-survival, then what difference it makes to use one vs. the other?
On Jun 27, 2011, at 4:31 PM, array chip wrote:
Hi, I am wondering if anyone can explain to me if cumulative
incidence (CI) is
just 1 minus kaplan-Meier survival?
First tell us what you think CI is defined as. I suspect it is not the
same. The KM estimator is cumulative product of
Hi,
I have a 1,785,421 rows x200 variables dataset with some missing values.
Approximately 55% of 1,785,421x200 are missing cells.
In this
ftp://ftp.stat.berkeley.edu/pub/users/breiman/Using_random_forests_v4.0.pdf
document , it is claimed that random forests can impute with great accuracy
even
On 27.06.2011 20:48, wwreith wrote:
Any ideas of what is happening when I attempt to load XLConnect? I think I
installed everything correctly. I posted what I get when I installed rJava
as well as when I attempted to load rJava directly. Thanks in advance for
any help you can give on this.
On 27.06.2011 20:28, scaramouch wrote:
Thank you for the answer,
If you thank me and ask me a follow up, please also send a message to
me, not only to the list.
Please also quote the original message, not all list members (and
particularly not myself) can remember the original stuff in
Hi David,
Thanks for responding, and plain text ...(didn't realized I was in rich text).
The endpoint is disease recurrence, I was producing a regular KM plot of
recurrence-free probability. Then someone recommend using cumulative incidence
is preferred because death was censored in the
I am taking a basic statistics course this summer, and while the majority of
the class is using a statistical package that came with the book, I am doing
everything in R for practical reasons. Forgive me if there is
documentation/instruction easily available on this topic, but Google and
forums
Dear R users,
Please consider the following example:
tmp - matrix(rnorm(100),10,10)
tmp[upper.tri(tmp, diag = TRUE)] - NA
heatmap(tmp, col = heat.colors(200), Rowv=NA, Colv=NA)
I have three questions:
(1) How would I change the axis labels from 1:10 to c(apples, pears,
bananas, kiwis,
#sampling elements
x-c(rep(0, 20), seq(1:37))
#number of simulations to perform
sims-100
#vector to store results
results-c()
#using for loop to perform simulations
for(i in 1:sims){
#take your sample
y-sample(x,3,replace=FALSE)
#check if all elements are zero
Hi robcinm,
You might also consider:
# data
x - c(rep(0, 20), 1:37)
# number of simulations
B - 1000
# result: TRUE/FALSE
out - replicate(B, {
y - sample(x, 3, replace = FALSE)
all(y == 0)
})
mean(out)
HTH,
Jorge
On Mon, Jun 27, 2011 at 5:08 PM, robcinm wrote:
I am
The first problem is that you are using a character string as the first
argument to agnes()
The help information for agnes says that its first argument, x, is
x: data matrix or data frame, or dissimilarity matrix, depending
on the value of the 'diss' argument.
Not a character
Here's a one-liner. Let's see their software do that!
sum(replicate(100, sum(sample(c(rep(0, 20), seq(1:37)), 3, replace =
FALSE)) == 0))
Am 27.06.2011 23:08, schrieb robcinm:
I am taking a basic statistics course this summer, and while the majority of
the class is using a statistical
Dear group,
When running the installation test:
testInstalledPackages(both,outDir='c:/Test')
I got the following message:
Running ‘testci.R’
comparing ‘testci.Rout’ to ‘testci.Rout.save’ ...
files differ in number of lines:
Please note the test does not result in 'OK' as do the other
As long as all the fields in each entry are the same, then the
following will work
df-data.frame(photo_id=NA,owner=NA,secret=NA,server=NA,farm=NA,title=NA,
ispublic=NA,isfriend=NA,isfamily=NA,tags=NA)
x-'photo id=5876248819 owner=13716719@N04 secret=faf9bb7f52
server=5264 farm=6
Hi everyone,
I am trying to find a way to filter a table; If I am given for example the
following table:
head(intra)
chr miRNA start end strand ACC hsa_ID
region region_start region_end gene_id transcrip_id
1 chr1 miRNA 1102484 1102578 + ACC=MI342;
Hi Dennis:
I tried your suggestions and I am getting the following errors:
Error in x$CP1 : $ operator is invalid for atomic vectors
In addition: Warning message:
In sqrt(diag(varcov)) : NaNs produced
Peter Maclean
Department of Economics
UDSM
- Original Message
From: Dennis Murphy
Hi,
Consistent with recent reports
(http://tolstoy.newcastle.edu.au/R/e14/devel/11/06/0434.html) I have found what
I believe is an error new to R 2.13.0 when reading a UTF-8 encoded file.
Minimal example:
=== filename: test.txt =
1 2 3 Я 4
Dear R Users,
I was using R to import several years traffic data, but every time after
data successfully imported (no error or warning) when I tried to save this
workplace or tackle these data, R GUI would automatically shut down. Would
you have any ideas that how I should deal with this
Hi all,
I have the next problem: I have a matrix with size 8,000,000x18. My personal
computer...blocks...so I have cut my original file into 100 different file.
I have written a function that should be run on each of this file.
So imagine
I need to read data from q1 to q100 file
On Jun 27, 2011, at 5:04 PM, array chip wrote:
Hi David,
Thanks for responding, and plain text ...(didn't realized I was in
rich text).
The endpoint is disease recurrence, I was producing a regular KM
plot of
recurrence-free probability. Then someone recommend using cumulative
the backage lmomco can plot the L-moments ratio diagram
would it be possible to draw the TL moments diagram also?
--
Osama Abdelaziz Hussien
Department of Statistics
Faculty of Commerce
Alexandria University
Egypt
__
R-help@r-project.org mailing list
This looks something like what you want.
http://r.789695.n4.nabble.com/Reading-in-a-series-of-files-using-a-for-loop-td906101.html
--- On Mon, 6/27/11, Trying To learn again tryingtolearnag...@gmail.com wrote:
From: Trying To learn again tryingtolearnag...@gmail.com
Subject: [R] Executing the
Use instead of
On Mon, Jun 27, 2011 at 6:01 PM, Nanami 13
deconstructed.morn...@gmail.com wrote:
Hi everyone,
I am trying to find a way to filter a table; If I am given for example the
following table:
head(intra)
chr miRNA start end strand ACC hsa_ID
Thank you David. Didn't realized someone posted related topic this morning.
John
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Sent: Mon, June 27, 2011 3:38:36 PM
Subject: Re: [R] cumulative incidence
try sqldf:
require(sqldf)
x
x
1 1
2 2
3 6
4 5
y
s e
1 1 4
2 3 5
3 8 10
sqldf(
+ select x.x
+ from x, y
+ where x.x between y.s and y.e
+ )
x
1 1
2 2
3 5
On Mon, Jun 27, 2011 at 5:33 AM, Salih Tuna saliht...@gmail.com wrote:
Thanks a lot. That solves the
On Jun 27, 2011, at 6:56 PM, jim holtman wrote:
Use instead of
Yes, and don't forget the comma!
--
David.
On Mon, Jun 27, 2011 at 6:01 PM, Nanami 13
deconstructed.morn...@gmail.com wrote:
Hi everyone,
I am trying to find a way to filter a table; If I am given for
example the
Dear R list.,
I am running a script to get a compact letter display.
library(lme4)
library(multcomp)
library(gplots)
### Mixed Effects Model #
data - read.table(AJmix.txt,header=TRUE, sep=\t)
attach(data)
y-cbind(positive,negative)
treatment-factor(treatment)
Hi all,
I'm trying to use a parallel script on a 20 node cluster, 8 processors
per node. Each node has a name, e.g., vm, vm0001, etc. Within
a foreach loop, I would like R to tell me what node it is actually
running on. How can this be accomplished? Thanks!
require(doMPI)
cl -
Here is another solution if you want to count a sequence like 0 as
three overlapping runs of 3 zeros.
x - sample(c(rep(0,20), rep(1,37)), 100, TRUE) # vector of 1M
z - rle(x) # get the runs
indx - which(z$values == 0 z$lengths 2) # find the runs
# now determine the number of 3
1 - 100 of 124 matches
Mail list logo