Hi,
I'm trying to perform biclustering to a matrix with genes as rows and
columns representing the frequency of each codon for each of its
aminoacids, like:
Tyr(TAT) Tyr(TAC) Iso(ATT) Iso(ATC) Iso(ATA)
Phe(TTT) Phe(TTC)
gene1 0.2 0.8 0.3 0.3
Hi all,
I just failed in loading a saved wordspace (13MB of size), and received
these errors:
Warning: stack imbalance in 'missing', 52 then 51
Warning: stack imbalance in 'if', 50 then 53
Warning: stack imbalance in 'as.environment', 57 then 59
Warning: stack imbalance in 'ls', 54 then 53
I used the information from this post to create a histogram with logarithmic
scale:
http://stackoverflow.com/questions/1245273/histogram-with-logarithmic-scale
However, the output from plot looks nothing like the output from hist. Does
anyone know how to configure the output from plot to resemble
|Dear list,
I have written some R-functions and build a R4CouchDB package for basic
R - couchdb interaction available at:
git://github.com/wactbprot/R4CouchDB.git
Regards
wactbprot
__
R-help@r-project.org mailing list
For some reason it now loads properly.
I tried it in several ways and wasn't able to reproduce the problem.
Sorry!
Tal
Contact
Details:---
Contact me: tal.gal...@gmail.com | 972-52-7275845
Read me: www.talgalili.com (Hebrew)
On 01/21/2010 03:53 AM, Lars Skjærven wrote:
Dear R users,
I'm making multiple plots within the same pdf page (par(mfcol = c(5,1)), and
want a legend for this at the bottom of all the plots. From previous mails
it has been suggested to use par(xpd=TRUE), increase the margin at the last
plot, and
I want n0 side by side then a gap and then n96 like in the first script
barplot(bio,beside=T)
I have done it in Excel but I had to put in some artifical zeros to produce
the gap and i will have to change the axis labels, but i was wondering if
there was a neater way of doing it in R
I've try it...
But don't work.
data1
1 2 3 4 5 6
5 0.2481509 0.2481509 1.240754 0.7444526 0.00 0.9926034
102.5269063 1.2634532 1.263453 2.5269063 1.263453 2.5269063
150.000 6.6856351 3.342818
On 21/01/2010 4:41 AM, Tal Galili wrote:
Hi all,
I just failed in loading a saved wordspace (13MB of size), and received
these errors:
Those messages indicate an internal error in R or a package that it's
using. If you can make such things reproducible, you should report them
as bugs,
Work's for me:
sapply(unique(union(names(data1), names(data2))),
function(n)Reduce('+', m[grep(n, names(m))]))
1 2 3 4 5 6
[1,] 0.2481509 0.2481509 1.240754 0.7444526 0.00 0.9926034
[2,] 2.5269063 1.2634532
Hello !
I have a quick question about loops.
For example, I have a 1 * 1 square and its inscribed circle tangent i, whose
radius, of course, is also 1.
The loop here is as the following:
generate n points, say 5, in the square randomly repeatedly until we have five
in total in the
On Jan 21, 2010, at 1:18 AM, Carrie Li wrote:
Dear R-helpers,
I have a question about giving index.
Suppose that I have a vector say, id=1:5, and each of them have some
given
index but some of them may share the same index
Now, I have another vector, and I want to give it a index based what
) data[sample(1:nrow(data), replace =
TRUE),])
r.indicator - sapply(resamples, indicator)
but it didn't work.
Can anybody help?
Thanks a lot.
Tomas
__ Informacia od ESET NOD32 Antivirus, verzia databazy 4792
(20100121) __
Tuto spravu preveril ESET NOD32
Hi,
One approach could be a while loop as follows. Note that your circle
should have radius 0.5 if I understood the problem correctly.
N - 5
npoints - 0
ntests - 0
points.in.circle - matrix(NA, ncol=2, nrow=N)
while (npoints N) {
test.point - runif(2, -0.5, 0.5) # generate new point in 2D
On Jan 21, 2010, at 7:28 AM, Wolfgang Amadeus wrote:
Hello !
I have a quick question about loops.
For example, I have a 1 * 1 square and its inscribed circle tangent
i, whose radius, of course, is also 1.
The loop here is as the following:
generate n points, say 5, in the square randomly
:
original biasstd. error
t1* 0.1774194 -0.001594390 0.01902264
Thanks a lot.
Tomas
__ Informacia od ESET NOD32 Antivirus, verzia databazy 4792
(20100121) __
Tuto spravu preveril ESET NOD32 Antivirus.
[1]http://www.eset.sk
Hi Gabor,
Perhaps I am missing something in my file names that I just can't see.
Whenever I have had problems reading a file into R, it always has been
because of some typo that I have made.
## This is what I am trying to run. The worksheet is the third one.
Originally it was named Paste
Hi, everyone:
I ask for help about translating a stata program into R.
The program perform cross validation as it stated.
#1. Randomly divide the data set into 10 sets of equal size, ensuring equal
numbers of events in each set
#2. Fit the model leaving out the 1st set
#3. Apply the fitted
On Jan 21, 2010, at 8:44 AM, Jerry Floren wrote:
Hi Gabor,
Perhaps I am missing something in my file names that I just can't see.
Whenever I have had problems reading a file into R, it always has been
because of some typo that I have made.
## This is what I am trying to run. The worksheet
Hi David,
Thanks for your response, and I should warn you that both my programming and
statistical experience is very limited.
I am pretty sure Perl is not installed. I'll take care of that first thing.
Thank you so much.
Jerry Floren
Minnesota Department of Agriculture
--
View this message
Dimitri Liakhovitski ld7631 at gmail.com writes:
One follow up question - the proposed solution was (notice - this time
I am introducing one NA in data frame x)
x-data.frame(factor=c(b,b,d,d,e,e),values=c(1,NA,10,20,100,200))
x$std.via.ave-ave(x$values, x$factor,
On Wed, Jan 20, 2010 at 4:37 PM, Dimitri Liakhovitski ld7...@gmail.com wrote:
Hello!
I have a data frame with a factor and a numeric variable:
x-data.frame(factor=c(b,b,d,d,e,e),values=c(1,2,10,20,100,200))
For each level of factor - I would like to divide each value of
values by the mean
Thanks a lot, everyone!
On Thu, Jan 21, 2010 at 9:44 AM, hadley wickham h.wick...@gmail.com wrote:
On Wed, Jan 20, 2010 at 4:37 PM, Dimitri Liakhovitski ld7...@gmail.com
wrote:
Hello!
I have a data frame with a factor and a numeric variable:
Hi,
On Thu, Jan 21, 2010 at 8:55 AM, zhu yao mailzhu...@gmail.com wrote:
Hi, everyone:
I ask for help about translating a stata program into R.
The program perform cross validation as it stated.
#1. Randomly divide the data set into 10 sets of equal size, ensuring equal
numbers of events
hello,
I am doing PCA in R using some habitat factors, and I used the function
result1=rda() and result2=princomp(),then pick up scores of the result1 and
result2 using scores(),but the scores are significantly different,i do not know
the meaning of it.
Best wishes!
Hi Thanks for all your help
Its a little difficult to follow those examples as all seem so different and
its hard to see how I do what I want to my data from the help files but i'll
try...
--
View this message in context:
Take a look at the validate.lrm function in the rms package.
Note that the use of threshold probabilities results in an improper
scoring rule which will mislead you. Also note that you need to repeat
10-fold CV 50-100 times for precision, and that at each repeat you have
to start from zero
Thanks Gabor, David, and David for all of your help. I did get it to work
with this code that I modified from the R Wiki from Gabor's original
suggestion:
START OF CODE
library(RODBC)
### setwd Set the working directory ###
setwd(C:/Documents and Settings/jfloren/My
On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris gpet...@uark.edu wrote:
I have found a good refernce to be S Programming by Venables and
Ripley.
I'll take a look at this book.
Since S and R are not completely the same, there are delicate
differences between S and R, which an S book may be
Hi
About your question, I not so clear what your specific context is. So I create
a general one for you, hope this can help.
a-function(x,y)sum(x*(10^seq(0,(y-1
You type
a(2,4) then you have
a(3,5) you have 555
If what you need here is only ,
you can use
a-
and
On Thu, 2010-01-21 at 21:56 +0800, echo_july wrote:
hello,
I am doing PCA in R using some habitat factors, and I used the
function result1=rda() and result2=princomp(),then pick up scores of
the result1 and result2 using scores(),but the scores are
significantly different,i do not know the
Based on the information below C:\TestR is not a file. Its a directory.
Try any of these three:
read.xls(file.choose(), sheet = 3, pattern = NAPT ID, as.is = TRUE)
read.xls(C:\\TestR\\MAP_2009Round2_2002Soiltest9-21-09.xls, sheet =
3, pattern = NAPT ID, as.is = TRUE)
setwd(C:\\TestR)
Hi,
I have a data set about age from 2006 to 2009 for certain individual. I
would like to draw the kernel density estimation to age for each year. But I
don't know how to draw the different plots in one graph.
The data looks like:
YearID Age
20061 20
20062 30
2007
I know my question has been asked many times. Please, forgive me for asking it
once more.
It's time to update R version.
I run R on Windows, Linux, and Mac OS/X.
I think I remember on Windows R comes with the un-installer. So I uninstall it
first and then I install the latest version.
On Mac OS
David Winsemius wrote:
On Jan 18, 2010, at 10:55 PM, David Winsemius wrote:
On Jan 18, 2010, at 8:58 PM, Jim Silverton wrote:
Hello everyone,
I tried searching for a manual for CAMAN package for mixture models,
but was
unable to find one. Does anyone have the pdf manual for this
I am so happy about learning how to read in multiple Excel files, that I have
to try and make another improvement. I know what I have been doing is
clumsy, but it works. Hopefully, someone can suggest a more elegant
solution. As a novice, I have been using MS-Word and mail merge to write my
code.
Hi everybody!
To use some functions, I have to transform my dataset into a list, where
each element contains one group, and I have to prepare a list for each
variable I have (altogether I have 15 variables, and many entries per
factor level)
Here is some part of my dataset:
SPECSHORBONE
Hi,
I have a file containing gene expressions written using the tillingArray
package. I used load() and attach() to get the data into R. Both of them works
fine. Now I want to see the contents of the file. How can I see the contents of
the file?
Thankyou for your time.
Regards,
Ambuj
## There are at least two options, one using mosaic and the other
using barchart.
## Here are both.
## this is your data
ld-matrix(c(25,25,50,10,30,60,15,35,50,30,30,40),nrow=4,byrow=T)
ntreat2-c(n0,n96,n0,n96)
rownames(ld)-ntreat2
typ-c(gr,fo,le)
colnames(ld)-typ
## now some rearrangements
mau...@alice.it wrote:
I know my question has been asked many times. Please, forgive me for asking it
once more.
It's time to update R version.
I run R on Windows, Linux, and Mac OS/X.
I think I remember on Windows R comes with the un-installer. So I uninstall it
first and then I install the
Hi Gabor,
Thanks for all of your help and the additional code. Sorry, about the
confusion on my part about the directory and file names.
Since my first posts, I have simplified file names, directory names, and the
spreadsheets by removing all the worksheets but the one needed. I tried this
I have a SpatialPolygonsDataFrame object (sp package).
It contains 40 polygons.
Now I want to join some of the polygones.
I cannot figure out how to do that? Any ideas?
The code below produces a map of Germany and I need to join some
counties.
library(sp)
library(RColorBrewer)
# get spatial
The error message says that it can't find perl on your system so
either you don't have perl or its not on your path. If you don't have
it at all install it. If its not on your path either add it to your
path or use the perl= argument to read.xls to tell it where to find
perl. One of the examples
On Wed, 20 Jan 2010, Robin Jeffries wrote:
I have been through the help file archives a number of times, and still
cannot figure out what is wrong.
I have a tab-delimited text file. 76Mb, so while it's large.. it's not
-that- large. I'm running Win7 x64 w/4G RAM and R 2.10.1
When I open this
On Jan 21, 2010, at 11:00 AM, mau...@alice.it mau...@alice.it wrote:
I know my question has been asked many times. Please, forgive me for
asking it once more.
It's time to update R version.
I run R on Windows, Linux, and Mac OS/X.
I think I remember on Windows R comes with the un-installer.
Just use
what=''
in your scan since all your data appears to be character.
Also use
comment.char='', quote=''
if there is the possibility of a misplaced quote or a comment
character. I have not problem reading in files of that size.
Also use 'count.fields' to see what your file looks like.
Dear R-helpers,
I have encountered the following unexpected behaviour of R-2.10.1, but not
R-2.9.0,
on both RHEL 4 and Ubuntu Karmic (precompiled via synaptic or built from
source).
I have a character vector from which I want to extract a certain pattern
that is surrounded
by junk as in:
nn
Dear group,
I want to create a plot similar to cdplot or cd_plot I vcd :
http://addictedtor.free.fr/graphiques/graphcode.php?graph=120
with lattice and support for multiple panels. Densities should not be
computed, these are stored in a data frame as in the example. My current
workaround with
Greetings to all,
I have presence/absence data for species in plots in different years. My
matrices looks like that:
YEAR1
sp1 sp2 sp3 sp4
11 1 0 1
21 1 1 1
31 1 0 1
41 1 0 1
51 1 0
Have you timed both solutions? Are you asking us to do you
analysis/debugging? It should be easy enough for you to setup a test
case and try it. That is what I would have to do to answer your
question and I don't use bloomberg.
On Wed, Jan 20, 2010 at 2:56 PM, anna lippelann...@hotmail.com
Hello,
I know there must be a simple soluton to this problem but it eludes me
currently.
My data is partitioned into two subsets, each subset has a common column
factor but with varying levels:
levels(fdf_ghc$AgeDemo)
[1] 26TO35 36TO45 46TO55 56TO65 66TO75 76TO85
levels(fdf_ghcnull$AgeDemo)
[1]
Hi,
On Thu, Jan 21, 2010 at 11:12 AM, AMBUJ ambuj_...@yahoo.co.in wrote:
Hi,
I have a file containing gene expressions written using the tillingArray
package. I used load() and attach() to get the data into R. Both of them
works fine. Now I want to see the contents of the file. How can I
Hi:
Perhaps this is what you were looking for:
densityplot(Age ~ Year, data = ages, group = ID, plot.points = FALSE)
HTH,
Dennis
On Thu, Jan 21, 2010 at 7:30 AM, Yi Du abraham...@gmail.com wrote:
Hi,
I have a data set about age from 2006 to 2009 for certain individual. I
would like to
Hi everyone,
I'm having some difficulty getting simple effects for the ols()
function in the rms package. The example below illustrates my
difficulty -- I'll be grateful for any help.
#make up some data
exD - structure(list(Gender = structure(c(1L, 2L, 1L, 2L, 1L, 1L, 1L,
2L, 1L, 2L, 2L, 2L, 1L,
One way is :
dataset = data.table(ssfamed)
dataset[, whatever some functions are on Asfc, Smc, epLsar, etc ,
by=SPECSHOR,BONE]
Your SPECSHOR and BONE names will be in your result alongside the results of
the whatever ...
Or try package plyr which does this sort of thing too. And sqldf may
Ista Zahn wrote:
Hi everyone,
I'm having some difficulty getting simple effects for the ols()
function in the rms package. The example below illustrates my
difficulty -- I'll be grateful for any help.
#make up some data
exD - structure(list(Gender = structure(c(1L, 2L, 1L, 2L, 1L, 1L, 1L,
2L,
On Thu, Jan 21, 2010 at 9:05 AM, Dieter Menne
dieter.me...@menne-biomed.de wrote:
Dear group,
I want to create a plot similar to cdplot or cd_plot I vcd :
http://addictedtor.free.fr/graphiques/graphcode.php?graph=120
with lattice and support for multiple panels. Densities should not be
Hello R-help group,
I have a question about merging lists. I have two lists:
Genes list (hSgenes)
namechrstrandstartendtransStarttransEnd
symboldescriptionfeature
ENSG02239721111874144121187414412
DEAD/H box polypeptide 11
but I have thousands of results so it would be really hand to find away of
doing this quickly
its a little difficult to follow those examples
Given your data in data.frame DF, maybe add the following to your list to
investigate :
dat = data.table(DF)
dat[, cor(Score1,Score2),
Trying to install package 'filehash' I get the following error on
FreeBSD 9.0-CURRENT (amd64) with R version 2.11.0 (2010-01-15 r50990):
---
R CMD INSTALL filehash_2.0-1.tar.gz
* installing to library '/usr/local/lib/R/library'
* installing *source* package
Hi,
Is there anyway to extract the AIC and BIC from the summary statistics
in fGarch package?
g - garchfit(data)
summary(g) will display the AIC and BIC, but I can't find a way to
store those in a variable (in order to loop through different
parameters, investigating their effect).
Thank you,
Hi,
1. Please keep replies on list -- note that you have to click reply all,
not just reply
2. For your new problem:
On Thu, Jan 21, 2010 at 1:22 PM, AMBUJ ambuj_...@yahoo.co.in wrote:
Hi Steve,
Thankyou for the reply. I still did not get the data from the
file.I executed the
Hi David,
str(g) gives you a ton of output, and the @fit slot has a $ics
component, part of which has the promising name of AIC...
(g...@fit)$ics[1]
HTH,
Stephan
David Rubins schrieb:
Hi,
Is there anyway to extract the AIC and BIC from the summary statistics
in fGarch package?
g -
Dear all,
I am having difficulty to built a model of quarter sales of spirits data, and
deciding which is the best model. The yfit2, yfit3, and yfit4 lines was not
appeared right at the end. The data and script is enclosed with this email.
I am using the harmonic regression model to
Tena koe Jerry
I'm not sure exactly what you wish to do, but it would seem you could
create a character vector (myVars) containing all your merge field names
and then use a for() loop. Something like:
for (myV in myVars)
{
x116 - subset(napt, Analysis_Soil %in% myV)
if (nrow(x116) 8)
{
On Jan 21, 2010, at 12:38 PM, b k wrote:
Hello,
I know there must be a simple soluton to this problem but it eludes me
currently.
My data is partitioned into two subsets, each subset has a common
column
factor but with varying levels:
levels(fdf_ghc$AgeDemo)
[1] 26TO35 36TO45 46TO55
I found this paper on ANOVA on unequal error variance. Has this be
incorporated to any R package? Is there any textbook that discuss the
problem of ANOVA on unequal error variance in general?
http://www.jstor.org/stable/2532947?cookieSet=1
__
Classification: UNCLASSIFIED
Caveats: NONE
Dear R users,
How can I retrieve an evaluated gradient value from deriv function
provided below? I want to retrieve b0 value. Appreciate your help.
Kyong
junk1-deriv(~(1/b1)*k-b0/b1,b0,c(k,b0,b1),formal=T)
junk1(k=0,b0=-14.0236,b1=2.44031)
[1]
I got the usual regression output from a model i did in R.
Im looking for alternative test criteria that can be done in R which are not
in the standard lm output.
I need to:
a) asses the fit of a logistic regression model and of a linear regression
model
b) compare the fit of models against
Try this:
attr(junk1(k=0,b0=-14.0236,b1=2.44031), 'gradient')
On Thu, Jan 21, 2010 at 6:05 PM, Park, Kyong H Mr CIV USA AMC
kyong.ho.p...@us.army.mil wrote:
Classification: UNCLASSIFIED
Caveats: NONE
Dear R users,
How can I retrieve an evaluated gradient value from deriv function
I have a 3d scatter plot and I need to add a trend line. I can add a
regression plane, but I need a line, not a plane.
Thanks.
--
View this message in context:
http://n4.nabble.com/3d-trend-line-tp1053686p1053686.html
Sent from the R help mailing list archive at Nabble.com.
On Thu, Jan 21, 2010 at 2:16 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
This paper was a prelude to his first book 'Exact Statistical Methods for
Data Analysis'.
He uses what is called a generalized p-value approach to inference, and for
the
book he wrote commercial software. AFAIK, no
Hi mlcarte3 (it is you?)
?contour
Try algo google R graph gallery
bests
milton
On Thu, Jan 21, 2010 at 2:43 PM, mlcarte3 mlcar...@ncsu.edu wrote:
I have a 3d scatter plot and I need to add a trend line. I can add a
regression plane, but I need a line, not a plane.
Thanks.
--
View this
Hi All,
please aplogize me if my qustion is a bit OT here,
but maybe is there someone that uses R from inside python
using rpy or rpy2 interface.
i'm at my first experience using rpy2.
i need to port torpy some R code .. but i'm having little problems
i hope yopu can help me.
The R code is :
On 21/01/2010 2:43 PM, mlcarte3 wrote:
I have a 3d scatter plot and I need to add a trend line. I can add a
regression plane, but I need a line, not a plane.
Which line do you want? There's more than one choice. For example, you
could regress y and z on x, or you could find the principle
Hey all,
I am just born today in the r language and I discover that it the language that
I need for my research.
I don't too much yet about this language. I want to know if I can find in this
language and install the
package data set related to The Framingham Heart Study, see references
Hello,
I'm fairly new to R and I can't work out how to produce a double
inequality like (LaTeX) $0 \leq x \leq 1$ in the legend of a graph. If
I try
legend(50, 0.1, legend = c(expression(0 = x = 1), c(2 = x = 3)), pch =
c(1,1), col = c(2, 3))
then I get an error message unexpected '=' in
Thank you Peter. I am really new to this. The spreadsheet I am working with
has 12,379 rows with the first row consisting of the variable names and
12,378 rows of data. There are seven columns, and the 7th column is the only
one with numerical data (Results).
I need to match up the variable
On 21/01/2010 3:50 PM, Mark Seeto wrote:
Hello,
I'm fairly new to R and I can't work out how to produce a double
inequality like (LaTeX) $0 \leq x \leq 1$ in the legend of a graph. If
I try
legend(50, 0.1, legend = c(expression(0 = x = 1), c(2 = x = 3)), pch =
c(1,1), col = c(2, 3))
then I
On Thu, Jan 21, 2010 at 2:41 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
On Thu, Jan 21, 2010 at 12:29 PM, Peng Yu pengyu...@gmail.com wrote:
On Thu, Jan 21, 2010 at 2:16 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
This paper was a prelude to his first book 'Exact Statistical
Read FAQ 3.3 as a starting place.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of Peng Yu
Sent: Thursday,
can you do Logistic regression in R, if so how do you do it and how do you
test the fit of a model?
--
View this message in context:
http://n4.nabble.com/Logistic-regression-tp1059870p1059870.html
Sent from the R help mailing list archive at Nabble.com.
Hi There,
If you are starting with S or R, don't be worried.
If you just started (I suppose with R), don't be worried too.
Just enjoy and be happy :-)
bests
milton
On Thu, Jan 21, 2010 at 10:44 AM, Peng Yu pengyu...@gmail.com wrote:
On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris
Thank you very much Duncan. I see I also had the mistake of a c
instead of expression in my question.
Mark
On Fri, Jan 22, 2010 at 7:59 AM, Duncan Murdoch murd...@stats.uwo.ca wrote:
On 21/01/2010 3:50 PM, Mark Seeto wrote:
Hello,
I'm fairly new to R and I can't work out how to produce a
?glm. Run the examples.
This might also help:
http://nlp.stanford.edu/manning/courses/ling289/logistic.pdf
DM
On Thu, Jan 21, 2010 at 1:17 PM, DispersionMap frenc...@btinternet.comwrote:
can you do Logistic regression in R, if so how do you do it and how do you
test the fit of a model?
--
Is there a way to display something like the quadratic equation, i.e. with a
numerator and demonator on a plot?
I think there is a way to create the equation in Latex, but wasn't sure if this
format would translate to plotting.
\begin{equation}
x = \frac{b +/- sqrt(b^2 - 4*a*c)}{2*a}
Hello R specialists,
I am using the runmean function to get the moving averages of a time serie.
Here is the code that I enter:
ma-runmean(na.remove(myTimeSerie),20,align=right,endrule=NA)
I get the vector that I want except for one detail: the first 20 values are
supposed to be NAs but I get a
Tena koe Jerry
myVars would be the unique values of Anlysis_Soil. I guess Anlysis_Soil
is a factor, in which case
myVars - levels(Anlysis_Soil)
The for() loop then steps through each of these in turn.
HTH
Peter Alspach
PS If you'd like more details it might be better to contact me
Thank you very much Peter. I'm off until Monday, but now I have an
interesting project to think about over the weekend. It looks pretty simple,
and it sure will beat over 2,000 pages of code to cut and paste.
Thanks again,
Jerry Floren
Minnesota Department of Agriculture
--
View this message
Hi:
Try this:
plot(c(0, 1), c(0, 1))
text(0.5, 0.5, expression(x == over(-b %+-% sqrt(b^2 - 4 * a * c), 2 * a)))
I've found it useful to print out the different plotmath expressions:
pdf('plotmath.pdf')
demo(plotmath)
dev.off()
HTH,
Dennis
On Thu, Jan 21, 2010 at 2:10 PM, Jason Rupert
Hi,
I have some code I run interactively through the R interpreter and it works
fine. I then run it as a script with Rscript and I get an error. The error
is coming when Rscript builds a model matrix.
Here is the Rscript code :
#!/usr/lib/R/bin/Rscript --verbose
require(MASS)
options(contrasts
Hello everyone, I am using the na.omit() function on a 785 x 79 matrix. When
I restrict the matrix to 10 columns for example it returns me the right
result with the na rows deleted: 756 x 10 matrix. But, if I enter the whole
matrix it returns me a 64 x 79 matrix, ununderstandable, I have no idea
[Apologies in advance if this is too statistics and not enough R.]
I've got an experiment with two sets of treatments. Each subject either received
all treatments from set A or all treatments from set B.
I can compute the N pairwise correlations for all treatments in either set using
cor(). If
Hello all,
I have data from a csv that i need to coerce into a multi dimensional
array. As you can see from my sample code i have yet to master building
efficient code. I can sort of do what i want but its very cumbersome code
and i know there is a better way to do it, i am just not
Tena koe Anna
Presumably because only 69 rows of your matrix have no missing values in
any column?
HTH
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of anna
Sent: Friday, 22 January 2010 12:01 p.m.
To:
No, the NA values are in the same rows and there are only 28 of missing value
rows...
--
View this message in context:
http://n4.nabble.com/Na-omit-on-matrix-does-the-matrix-have-to-have-a-limited-size-tp1068000p1069099.html
Sent from the R help mailing list archive at Nabble.com.
Hi:
Before you do all of that, might one inquire about which goal you are trying
to
achieve? Perhaps a 27-dimensional array may not be the wisest choice of data
structure for your goal.
Dennis
On Thu, Jan 21, 2010 at 3:06 PM, LCOG1 jr...@lcog.org wrote:
Hello all,
I have data from a csv
Please see the following example. I can not write '' to a csv file
successfully. Could somebody let me if it is possible to write '' to
a csv file with the default setting of write.csv?
my_home$ Rscript main_quote.R
x=rbind(
+ \A\
+ , \B\
+ )
x
[,1]
[1,] \A\
[2,] \B\
write.csv(x,
Hi Peng,
May be write.table(...,quote=F)
bests
milton
On Thu, Jan 21, 2010 at 6:29 PM, Peng Yu pengyu...@gmail.com wrote:
Please see the following example. I can not write '' to a csv file
successfully. Could somebody let me if it is possible to write '' to
a csv file with the default
On Thu, Jan 21, 2010 at 5:32 PM, milton ruser milton.ru...@gmail.com wrote:
Hi Peng,
May be write.table(...,quote=F)
bests
milton
Would this cause some trouble when filed is missing, etc.? Although I
don't have a concrete example, I think that it might happen. So I
prefer to have quote=T.
1 - 100 of 119 matches
Mail list logo