Dear all
I want to run a simple multiple linear regression model, say:
mod - lm(y~x1+x2+x3+x4)
I have programmed a custom algorithm to calculate the p-values, because
there is clustering present.
I ordered the variables according to their p-value and want to stepwise
eliminate the variables
On Mon, Mar 21, 2011 at 02:28:32PM -0700, armstrwa wrote:
Hi all,
Forgive me for this basic question. I've been doing some research and
haven't been able to figure out how to best do this yet.
I have 75 variables defined as vector time series. I am trying to create a
script to automate
Hello,
I'm wondering if someone may be able to help me, and do apologize if there is a
simple and obvious solution for this. I am somewhat new to R, and have been
searching for a simple solution for a couple of days.
I am interested in finding a tool that allows me to plot a stacked bar
On 2011-03-21 14:39, joe82 wrote:
Hello All,
I need help with my dataframe, it is big but here I am using a small table
as an example.
My dataframe df looks like:
X1 X2X3
1 2011-02 0.00 96.00
2 2011-02 0.00 2.11
3 2011-02 2.00 3.08
4 2011-02 0.06 2.79
5
Hello,
I am running large simulations, which unfortunately I can't really
replicate here because the code is so extensive. I rely heavily on
mclapply, but I realize that I'm losing data somewhere.
There are two worrisome symptoms:
1) I am getting 'NULL' as a return value for some (but not
On 2011-03-21 14:16, armstrwa wrote:
Hi,
I am running a correlation analysis on a temporal dataset. I was wondering
if you would receive the same tau and p values running the function:
MannKendall(x), where x is the dependant variable that changes with time
as you would running:
require(ppls)
data(BOD)
X-BOD[,1]
y-BOD[,2]
Xtest=seq(min(X),max(X),length=200)
dummy-X2s(X,Xtest,deg=3,nknot=20)
Z-dummy$Z
Ztest-dummy$Ztest
size-dummy$sizeZ
P-Penalty.matrix(size,order=2)
lambda-200
number.comp-3
penalized.pls(Z,y,P=lambda*P,ncomp=number.comp)$coefficients # By default
kernel=F
Hello
I want to filter a data frame called 'filtEXP' for rows where column
'max' 6 OR column 'FC' 4 and am trying various permutations of the
below to no avail:
filtEXP2 - filtEXP[filtEXP$max 6,] || filtEXP[filtEXP$FC 4,]
Any ideas?
W dniu 22 marca 2011 10:27 u¿ytkownik £ukasz Rêc³awowicz
lukasz.reclawow...@gmail.com napisa³:
penalized.pls.kernel(Z,y,M=lambda*P,ncomp=number.comp) # But using
directly, coefficients are different.
I see me error P=!M, but still results are the same...
p - ncol(Z)
Minv - diag(p) + P
M -
Try this
x - seq(0,15,0.01)
plot(x+1, dexp(x, 0.2), c(0,15), c(0,0.8),type=l, ylab=f(x))
lines(x+1, dexp(x, 0.8))
grid()
abline(v=1, col=lightgray, lty=dotted,)
Best
ep
--
View this message in context:
http://r.789695.n4.nabble.com/Exponential-distribution-tp3394476p3395829.html
Sent from
Hi,
Have installed Ubuntu packages for R by adding:
deb http:///bin/linux/ubuntu hardy/
to my /etc/apt/sources.list file and using:
sudo apt-get update
sudo apt-get install r-base
I understand that to enable shlib I need to download source files to my home
directory. When I run
apt-get source
On 2011-03-21 10:37, Savitri N Appana wrote:
Thank you for your suggestion Allan. I should have paid attention to
the posting instructions.
Pls find below the sample code from the ?splsda in the caret package.
Note: It used to work fine in R v2.8.1, but this error shows up now,
given that
subset(filtEXP,max 6 | FC 4)
Sent from my iPad
On Mar 22, 2011, at 5:27, Dry, Jonathan R jonathan@astrazeneca.com
wrote:
Hello
I want to filter a data frame called 'filtEXP' for rows where column
'max' 6 OR column 'FC' 4 and am trying various permutations of the
below to no
Hi Denis,
how about this:
to.drop-x1
update(mod,as.formula(paste(.~.-,to.drop,sep=)))
hth.
Am 22.03.2011 08:12, schrieb denis.ay...@unibas.ch:
Dear all
I want to run a simple multiple linear regression model, say:
mod - lm(y~x1+x2+x3+x4)
I have programmed a custom algorithm to
On 03/22/2011 06:30 PM, Chandra Salgado Kent wrote:
Hello,
I'm wondering if someone may be able to help me, and do apologize if there is a
simple and obvious solution for this. I am somewhat new to R, and have been
searching for a simple solution for a couple of days.
I am interested in
On 03/22/2011 07:54 AM, Jim Silverton wrote:
I am doing a histogram with 2 superimposed densities. However, the density
of one of the graphs is not coming out..its being erased.. Any ideas on how
to fix this problem?
Hi Jim,
It's probably a scale problem, but an example to show us what is
Greetings,
I've just put out the latest version of The Popularity of Data Analysis
Software at http://r4stats.com/popularity. This update includes complete data
for 2010, the addition of number of blogs for each software, more coverage of
Statistica, and, where possible, measures regarding the
Is there an easy way to control the font size for the contour lines
in plots of kde objects in package ks?
The label size seems not obey the cex parameter (which probably
would be the R way of doing it).
__
R-help@r-project.org mailing list
Hi Laszlo,
the see also part of ?assign comes in use here, it leads to get
get(paste(df,i,sep=.)
Cheers
Am 21.03.2011 17:14, schrieb Bodnar Laszlo EB_HU:
Hello everyone,
I'd like to ask you a question again, basically focusing on referring to
different objects.
Let's suppose we
On 22 March 2011 at 02:17, nick@R wrote:
| Hi,
|
| Have installed Ubuntu packages for R by adding:
| deb http:///bin/linux/ubuntu hardy/
| to my /etc/apt/sources.list file and using:
| sudo apt-get update
| sudo apt-get install r-base
|
| I understand that to enable shlib I need to download
Dear R users,
I like to see plots with the tickmarks covering all range of data.
For example, instead of
plot(rnorm(100)),
which gives y axis tick marks between -2 and 2, while the data has a wider
range
I usually do this:
x=rnorm(100); plot(x, ylim=range(pretty(range(x)))
where y axis tick
Hi,
I am trying to calculate Principal Component Scores per id per year using
the psych package. The following lines provide the scores per obeservation
pca = data.frame(read.table(textConnection(id year A B C D
1001 1972 64 56 14 23
1003 1972 60 55
Hi Umesh,
I use the plyr package for this sort of thing:
library(plyr)
daply(dataframe, .(ped), myfun)
Best,
Ista
On Tue, Mar 22, 2011 at 3:48 AM, Umesh Rosyara rosyar...@gmail.com wrote:
Dear R-experts
Execuse me for an easy question, but I need help, sorry for that.
From days I have been
Hi Chandra,
You could use ggplot2:
library(ggplot2)
ggplot(dat, aes(Group_size, Number, fill=Sex)) +
geom_bar(stat=summary, fun.y=mean)
Best,
Ista
On Tue, Mar 22, 2011 at 7:30 AM, Chandra Salgado Kent
c.salg...@cmst.curtin.edu.au wrote:
Hello,
I'm wondering if someone may be able to help
Thx Petr. It worked likea charm. Regards.
__
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 commented, minimal, self-contained,
Hi Deepayan,
Meanwhile I found another solution by simply creating more steps within the
colorscale, so one could not notice when the colours mismatch one or two
units. Your solution is more elegant though.
Joep
--
View this message in context:
Dear All,
I am an Italian researcher in Economics. I work with large sample data. I
need to increase the memory in R-project in order to upload a file .dta.
How can I do this?
Thank you.
graziella
--
View this message in context:
Thank you all !!
On Tue, Mar 22, 2011 at 3:53 AM, Peter Ehlers [via R]
ml-node+3395720-1679132054-222...@n4.nabble.com wrote:
On 2011-03-21 14:39, joe82 wrote:
Hello All,
I need help with my dataframe, it is big but here I am using a small
table
as an example.
My dataframe df
Dear all. I have to plot a the marginal population density for a heat map that
represents the population density of a city. I have been able to plot the heat
map in the lower left corner, the marginal density in x in the upper left
corner and the marginal density in y in the lower left corner.
Hi,
I am looking for a way to study some phase data with a circular distribution
measured in rad.I would like to do a two way ANOVA (if possible mixed, with
inter and intrasubject).I haven´t found a package that does that in R?Does
sombeody know if there is one or how to do the analysis.Thanks
R experts,
Hello, I am trying to sample a vector 1:40 without replacement such that no
element in the new vector
is within 7 units of either of its immediate neighbors. This is part of a
larger program I am working
on. The following code works well about 65 % of the time (14/40). The problem
Dear List,
I have a data frame with approximately 50 rows that looks like this:
Date time value
…
19.07.1956 12:00:00 4.84
19.07.1956 13:00:00 4.85
19.07.1956 14:00:00 4.89
19.07.1956 15:00:00
1) What kind of file is a .dta?
2) Do you have more RAM than R is using? Check this with ?memory.limit
3) Are you on a 64 bit OS with 64 bit R?
I can tell you that starting R with --max-mem-size=
may help, but you probably want to look into R interfaces with databases
first.
On Tue, Mar 22, 2011 at 3:05 PM, Tonja Krueger tonja.krue...@web.de wrote:
Dear List,
I have a data frame with approximately 50 rows that looks like this:
Datetimevalue
19.07.1956 12:00:00 4.84
19.07.1956 13:00:00 4.85
On Tue, Mar 22, 2011 at 11:05 AM, Tonja Krueger tonja.krue...@web.de wrote:
Dear List,
I have a data frame with approximately 50 rows that looks like this:
Date time value
…
19.07.1956 12:00:00 4.84
19.07.1956 13:00:00 4.85
Dear All,
I need to introduce new distribution for the error components of a frontier
model. For this reason, I use the package Frontier 4.1, but I need to
introduce other distribution function than Normal and Half-Normal (that are
those given by default).
I hope that someone help me.
Thanks.
Got it working with a package called 'xlsReadWritePro' (which requires a
license after 30 day's...) This is my code:
setwd(J:/Stage/Datasets2/Datasets/outData)
masterTable-read.table(AR1000900A_N_241110_(Mapping250K_Nsp)_2,Mapping250K_Nsp,CNprobes.tab
_SNP_IDs.xls,sep=\t, dec=,, fill=T, header=T)
Thank you Jonathan.
A file .dta is a dataset uploaded on Stata 10, that is another Statistical
Program of data processing.
I think that you help me.
Thank you very much.
graziella
--
View this message in context:
http://r.789695.n4.nabble.com/memory-increasing-tp3396511p3396679.html
Sent from
Dear R community members:
The print edition of the R Cookbook is now available:
http://www.amazon.com/dp/0596809158/
The Cookbook is also available from the publisher in a variety of e-book
formats, including ePub, mobi, and PDF
(http://oreilly.com/catalog/9780596809157).
The book is a
Thank you very much Peter. It works fine now
Best,
Savi
Peter Ehlers ehl...@ucalgary.ca 3/22/2011 5:49 AM
On 2011-03-21 10:37, Savitri N Appana wrote:
Thank you for your suggestion Allan. I should have paid attention
to
the posting instructions.
Pls find below the sample code from
I have a 45 x 16 data frame consisting of dissimilarities among 10
colors, giving in each
column the 45 = 10*9/2 pairwise judgments for one of 16 subjects. The
rownames
identify each pair of colors, e.g, AC = (A,C), and the pairs are
ordered by columns
in the lower triangle of each distance
I don't know if this is what you're looking for but it describes how to do a 2
way repeated measures with R. You problem may be different and I lack the
stats knowledge to know that. If that's the case I apoligize:
1. I suggest you try posting on the r-sig-geo list, as I would guess
that circular distributions are a relevant topic of concern there
(spatial statistics) for which expertise would be available.
2. Also check out the lme4 package and glmer. You might be able to
write your own family function
Dear Graziella
On 22 March 2011 16:08, graziella graziella.bona...@unical.it wrote:
I need to introduce new distribution for the error components of a frontier
model. For this reason, I use the package Frontier 4.1,
The software Frontier 4.1 is not related to R. However, the R
package frontier
On Mar 22, 2011, at 10:17 AM, graziella wrote:
Dear All,
I am an Italian researcher in Economics. I work with large sample
data. I
need to increase the memory in R-project in order to upload a file
.dta.
How can I do this?
Hello,
First, you can try to split your data frame in this way:
list.year-split(pca, unique(pca$year))
And then apply the principal component analysis over the list list.year.
Regards,
Carlos Ortega
www.qualityexcellence.es
On Tue, Mar 22, 2011 at 12:35 PM, mathijsdevaan
filter(), in the stats package, can do moving
averages (with any weights).
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Tonja Krueger
Sent: Tuesday, March 22, 2011
On Mar 22, 2011, at 12:06 PM, David Winsemius wrote:
On Mar 22, 2011, at 10:17 AM, graziella wrote:
Dear All,
I am an Italian researcher in Economics. I work with large sample
data. I
need to increase the memory in R-project in order to upload a file
.dta.
How can I do this?
Hello,
I forgot to mention that I am looping over ~70K objects. If I do
mclapply on the first 200, its fine (i.e. doesn't give NULL values); if
I go up to 2K (or over all of them), then I start to see NULL values.
Also the function I call uses commands 'restrict', 'gaps' and 'width'
from
Michael -
I think this does what you want:
helm.raw -
read.table(http://euclid.psych.yorku.ca/datavis/Private/mdshelm.dat,header=TRUE,
row.names=1)
trans =
c('A'='RPur','C'='Red','E'='Yel','G'='Gy1','I'='Gy2','K'='Green','M'='Blue','O'='BlP','Q'='Pur1','S'='Pur2')
cnames =
The simple thing to do is to put a sanity counter in the 'repeat'
statement and if you have been through it a certain number of times,
then exit. Anytime you have a loop that might run forever, you should
have some sanity/watchdog timer on it.
On Tue, Mar 22, 2011 at 11:02 AM, Hosack, Michael
Hi all,
I am trying to write a script that will automate the task of running a
Kendall's Tau correlation test on 75 time series that I am interested in.
The code I have written is:
for(i in 1:length(gagehandles)){
dates-get(paste(gagehandles[i],_amsd,sep=))
lt;!DOCTYPE html PUBLIC quot;-//W3C//DTD HTML 4.01 Transitional//ENquot;gt;
Thank you very much, Peter.nbsp; That does make it clearer.
Billy
Peter Ehlers [via R] wrote:
On
2011-03-21 14:16, armstrwa wrote:
gt; Hi,
gt;
gt; I am running a correlation analysis on a temporal
lt;!DOCTYPE html PUBLIC quot;-//W3C//DTD HTML 4.01 Transitional//ENquot;gt;
Thanks, Petr.nbsp; Your insight has helped me out a lot.
Billy
Petr Savicky-2 [via R] wrote:
On Mon, Mar 21, 2011 at 02:28:32PM -0700, armstrwa wrote:
gt; Hi all,
gt;
gt; Forgive me for this basic
On 3/22/2011 12:47 PM, Phil Spector wrote:
Michael -
I think this does what you want:
*Perfectly*
helm.raw -
read.table(http://euclid.psych.yorku.ca/datavis/Private/mdshelm.dat,header=TRUE,
row.names=1)
trans =
Hi,
This might do what you want:
iter - 0
repeat {
iter - iter + 1
ss - numeric(40)
ss[1] - sample(1:40,1)
for (i in 1:39) {
## calculate all possible step sizes that will give a new value in
the 1:40 range
pmove - sample((1 - ss[i]):(40-ss[i]))
## drop all step sizes that
If you care about confidence interval coverage, type I error, or predictive
accuracy, trying different models in this way is not the way to go.
Frank
agent dunham wrote:
Dear all,
I want to improve my adj - R sq. I 've chequed some established models and
they introduce two times the
Thank you Jim, I was not aware of watchdog timers. R is my first computer
language. I searched through the archives and found the function setTimeLimit
and applied it with low values for cpu and elapsed time. It provides a way out
of the loop which is what I needed. It's not a counter, but I
Thank you so much Martyn. Now I can most likely avoid having to rerun my
program
multiple times.
Mike
-Original Message-
From: Martyn Byng [mailto:martyn.b...@nag.co.uk]
Sent: Tuesday, March 22, 2011 1:21 PM
To: Hosack, Michael
Cc: r-help@R-project.org
Subject: RE: [R] Infinite loop
armstrwa wrote:
for(i in 1:length(gagehandles)){
dates-get(paste(gagehandles[i],_amsd,sep=))
Error in get(paste(gagehandles[i], _amsd, sep = )) :
variable names are limited to 256 bytes
It didn't have a problem with the variable names before, so I'm not sure
what's going on
Dear List,
I've used R fr a while now, but I'm totally unfamiliar with the compiling
process. In particular, I'd like to compile R for AIX 5.3. servers. I'll
appreciate your guidance on the best place to get started on this.
Thanks in advance,
Lars.
[[alternative HTML version deleted]]
Hi there
I recently got some affymetrix human gene 1.0 st array. I only get out the
probe level data and I wonder if there is a function that is able to average
the probe level data for each gene?
Best wishes
Kristian Unger
Helmholtz Zentrum M?nchen
Deutsches
On Mar 22, 2011, at 2:35 PM, Unger, Kristian, Dr. wrote:
Hi there
I recently got some affymetrix human gene 1.0 st array. I only get
out the probe level data and I wonder if there is a function that is
able to average the probe level data for each gene?
The BioConductor list is probably
I have 2 questions concerning the EM algorithm. Is it true that the
EM algorithm gives unique answers for the means and variances of a mixture
of 2 normals? I am using mixtools and I am surprised that it works better
than a Bayesian program I wrote.
If so can someone say why the mixing
I have a dataframe that looks like this:
str(chr)
'data.frame': 84 obs. of 7 variables:
$ county: Factor w/ 3 levels Broome,Nassau,..: 3 3 3 3 3 3 3 3 3 3 ...
$ item : Factor w/ 28 levels Access to healthy foods,..: 21 19 20
18 16 3 2 6 17 8 ...
$ value : num 8644 15 3.5 3.9 7.7 ...
$
In general when you want to split up your data and do the same thing on each
piece then combine the results back together it is good to look at the plyr
package. But for this specific case you should look at the lmList function in
the nlme package which may do exactly what you want with the
Bert Gunter gunter.berton at gene.com writes:
1. I suggest you try posting on the r-sig-geo list, as I would guess
that circular distributions are a relevant topic of concern there
(spatial statistics) for which expertise would be available.
2. Also check out the lme4 package and glmer.
On 2011-03-22 03:57, Erich Neuwirth wrote:
Is there an easy way to control the font size for the contour lines
in plots of kde objects in package ks?
The label size seems not obey the cex parameter (which probably
would be the R way of doing it).
Erich,
Try the argument 'labcex'.
Peter
Hi R-users,
I'm trying to built a plot of two series of data, but thees series
result superimposed. The R-code is like this:
goal - c(125, 143, 81, 26, 2, 3)
numgoal - 0:5
lambda - sum(goal*numgoal)/sum(goal)
plot(numgoal, goal, type=h)
x - 0:5
y-dpois(x, lambda)
att-y*380
Thank you very very much... It's perfect to my goal...
--
View this message in context:
http://r.789695.n4.nabble.com/Exponential-distribution-tp3394476p3396803.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org
Well im thoroughly frustrated after 25 minutes of checking and rechecking my
path. What do i not know about loading a mdb that is keeping me from
loading my data.
i have loaded the Hmisc library and pointed it too my data using mdb.get and
continue to get the following error.
Error in
This is a reminder that the deadline for abstract submission and early bird
registration for the R User Conference 2011 is April 1st.
useR! 2011
http://www.R-project.org/useR-2011
A conference centred on the use of R for data analysis and statistical
computing.
August 16-18, University of
Dear all,
I want to improve my adj - R sq. I 've chequed some established models and
they introduce two times the same variable, one transformed, and the other
not. It also improves my adj - R sq.
But, isn't this bad for the collinearity? Do I interpret coefficients as
usual?
After four years of using R, I finally have run into a problem to which I can't
find a solution in the guides or forums and thus I am making my first post.
Our lab has Fortran code for population modeling. I have been using R as a
wrapper to process the raw data, generate a batch script that
Hi,
I have a data.frame(zscores) that looks like this:
gA gB
g1 0.20.6
g2 0.3Na
My problem is that I need to use a function and the output is a vector
of only the non NA values, so shorter than the list I would obtain
dropping the data.frame.
What is the cleanest way to
I believe that you are misinterpreting the error message.
The mdb-get function relies on external programs named
mdb-tables and mdb-schema to do its work. Those programs
must be installed on your computer (and in directories in
your search path), for the function to work.
I don't know what
Hi:
Try this instead of your original lines() statement:
lines(x + 0.1, att, type=h, col=red)
or
plot(numgoal, goal, type=h, lwd = 3)
x - 0:5
y-dpois(x, lambda)
att-y*380
lines(x + 0.1, att, type=h, col=red, lwd = 3)
legend(x=topright, legend=c(Osservate, Attese), col=c(black,red),
ncol=1,
Colleagues
R: 2.12.2
Windows 7 (64-bit)
I am experiencing unexpected behavior in the following situation:
FILE1:
windows()
plot(1, 1,)
FILE2:
windows()
PLOT- code.to.create.a.lattice.object()
On 3/22/2011 6:37 AM, Muenchen, Robert A (Bob) wrote:
Greetings,
I've just put out the latest version of The Popularity of Data Analysis
Software at http://r4stats.com/popularity. This update includes complete data for
2010, the addition of number of blogs for each software, more coverage of
Some systems (both hardware and software) may have watchdog timers;
I was just using it to illustrate a point as to one way of breaking
out. In the case of this loop, it probably run too fast to try and
put a timer on it; counting the iterations is good enough. There are
some ways of creating
I've used RODBC to read in ms access files... or if you're as lazy as
me you could use the following below (it can handle some other ms
office file types too and thinks it can recognize file types but as
has been pointed out in this list, using it with excel probably means
trouble)
read.mso -
I don't doubt that R may be the most popular in terms of discussion group
traffic, but you should be aware that the traffic for SAS comprises two
separate lists that used to be mirrored, but are no longer linked
Usenet -- news://comp.soft-sys.sas (what you counted)
listserve -- SAS-L
Josh -
There are two ways that I know of to use Access in R on
Windows:
1) odbcConnectAccess or odbcConnectAccess2007 in the ODBC package
and
2) Configuring the Access ODBC connector in
Control Panel - Administrative Tools - Data Sources(ODBC)
and using connectODBC in the
On Tue, Mar 22, 2011 at 11:35 AM, Unger, Kristian, Dr.
un...@helmholtz-muenchen.de wrote:
Hi there
I recently got some affymetrix human gene 1.0 st array. I only get out the
probe level data and I wonder if there is a function that is able to average
the probe level data for each gene?
In
Without seeing your code it could be just about anything.
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained reproducible code.
--- On Mon, 3/21/11, Jim Silverton jim.silver...@gmail.com wrote:
From: Jim Silverton
Hello, I've been searching on the web for a few hours and seem to be stuck on
this. The code pasted below generates a histogram of subject responses in four
different conditions in an experiment. This version of the graph is one I'm
using for internal consistency checking, so I've set it up
I am getting into some more posts about having to have mdb-tools installed on
my system (Windows XP) and have located the download but its not an executable
install file and it looks like from what you said these tools are for
non-windows based systems anyhow?
Is the mdb.get function the
I am trying the RODBC package and the odbcConnectAccess function but am getting
an error:
Error in
odbcConnectAccess(T:/Models/LandUse/GenericLandSupplyModel/Projects/2008BaseYear/EugeneUGB/
Reporting/Summary.mdb) :
odbcConnectAccess is only usable with 32-bit Windows
I am on a 64 bit
Thank you, Ista. It helps.
Best Regards
Umesh R
_
From: istaz...@gmail.com [mailto:istaz...@gmail.com] On Behalf Of Ista Zahn
Sent: Tuesday, March 22, 2011 8:58 AM
To: Umesh Rosyara
Cc: R mailing list
Subject: Re: [R] help need on working in subset within a dataframe
Hi
On Mar 22, 2011, at 4:14 PM, Dennis Fisher wrote:
Colleagues
R: 2.12.2
Windows 7 (64-bit)
I am experiencing unexpected behavior in the following situation:
FILE1:
windows()
plot(1, 1,)
FILE2:
windows()
Hi,
thanks for your reply.
* Gabor Grothendieck ttebguraqv...@tznvy.pbz [2011-03-18 17:51:03 -0400]:
On Fri, Mar 18, 2011 at 2:19 PM, Sam Steingold s...@gnu.org wrote:
2. An alternative that won't generate a warning message but involves a
double read is:
library(zoo)
toPOSIXct -
Well, a custom panel function is what you need (or one that may
already exist somewhere: try googling on high low intervals in R
graphs or some such).
So if you haven;t already done so, try Paul Morrell's Chapter on
lattice plots from his book for how panel functions work:
On Tue, Mar 22, 2011 at 6:49 PM, Sam Steingold s...@gnu.org wrote:
Hi,
thanks for your reply.
* Gabor Grothendieck ttebguraqv...@tznvy.pbz [2011-03-18 17:51:03 -0400]:
On Fri, Mar 18, 2011 at 2:19 PM, Sam Steingold s...@gnu.org wrote:
2. An alternative that won't generate a warning message
Dear All,
I would like to generate random protein sequences using a HMM model.
Has anybody done that before, or would you have any idea which package
is likely to be best for that?
The important facts are that the HMM will be fitted on ~3 million
sequential observations, with 20 different states
Thank you all for the help. After I used proper subsetting, it all worked
fine.
And yes, I am open to receiving a prize for the most ancient version of R in
use...
Best,
alina
--
View this message in context:
http://r.789695.n4.nabble.com/problem-with-xyplot-tp3391822p3398070.html
Sent from
Subject: Re: [R] Rapache ( was Developing a web crawler )
From: m...@biostatmatt.com
To: marchy...@hotmail.com
CC: r-help@r-project.org
Date: Sun, 6 Mar 2011 13:51:53 -0500
On Sun, 2011-03-06 at 08:06 -0500, Mike Marchywka wrote:
You might also consider the Deducer package. You can build up a plot by
point and click and then have a look at (and amend) the code and learn the
syntax of ggplot2, which is a nice alternative to the lattice package.
The website of the Deducer package (www.deducer.org) is a good start.
--
Dear John, Peter and others,
So, I now have a query at an even more elementary level and that is
regarding my results from anova.mlm() not matching the car package's
Manova(). Specifically, I have been trying the following out with regard
to a simple one-way MANOVA setup. So, I try out the
I have a dataframe with many rows like this:
df
X1 X2 X3 X4 X5 X6 X7 week d
sim1 FALSE TRUE TRUE TRUE TRUE TRUE TRUE1 0.3064985
sim1 is the rowname, X1..X7,week,d are the column names. X1..X7 are factors,
booleans in this case.
I need to add another
You need a list object indeed of a vector, try this:
rbind(df, dreps = c(rep(list(TRUE), 7), 5, 0))
On Tue, Mar 22, 2011 at 10:12 PM, Alexy Khrabrov delivera...@gmail.com wrote:
I have a dataframe with many rows like this:
df
X1 X2 X3 X4 X5 X6 X7 week d
1 - 100 of 120 matches
Mail list logo
