Carl Witthoft wrote:
Just jumping to conclusions since the site is down as I type (7PM EST).
What's the news from that part of Australia -- thus showing my
complete ignorance of DownUnder geography.
Hi Carl,
The big fires are hundreds of kilometers from Newcastle, although there
was a
William Revelle wrote:
At 6:41 PM -0500 2/9/09, David Winsemius wrote:
describe() in Hmisc provides much of the rest of what you asked for:
describe(pref900$TCHDL)
pref900$TCHDL
n missing uniqueMean .05 .10 .25 .50
.75 .90 .95
9061904469 16051
Hi,
#There are a lot of examples for barplot if you just type
example(barplot)
# i altered one slightly to put the values on top of each bar
mp - barplot(VADeaths) # default
#tot - colMeans(VADeaths) #changed this line
tot = colSums(VADeaths)#wether you need max, min, mean , colSums
Hi José,
I wonder if there are any useRs sharing day-to-day realizations/tricks
on twitter...
Seems like a good place for those things that are good findings, but
one is too lazy to blog about them...
I don't twitter, but I wrote an R package
to read and write tweets from the R command
line
Hello, everyone
How do I replace dot with empty space in ED4.Comdty?
I need to get ED4 Comdty
tried sub() in many different ways, like sub({.}, , ED4.Comdty)
etc but could not do it.
Thanks in advance,
Sergey
__
R-help@r-project.org mailing list
OK, got it now:
sub([.], , ED4.Comdty)
On Tue, Feb 10, 2009 at 10:01, Sergey Goriatchev serg...@gmail.com wrote:
Hello, everyone
How do I replace dot with empty space in ED4.Comdty?
I need to get ED4 Comdty
tried sub() in many different ways, like sub({.}, , ED4.Comdty)
etc but could not
Dear R-Help-Team,
I tried to use isoMDS-Function of the MASS-Package:
Message:
Fehler: konnte Funktion isoMDS nicht finden
(error: could not find function isoMDS)
so I tried to install the package MASS:
utils:::menuInstallPkgs()
versuche URL
On Mon, 2009-02-09 at 06:40 -0800, Vie wrote:
Hi,
Ive been trying to find a function that will allow me to pull out a number
between a minimum and maximum threshold.
I want a random decimal number between, for example, 0 and 0.5 or 0 and 0.7.
I've been searching everywhere for a function
Hi all,
A colleague of mine wrote a little function to send messages to twitter:
update - function(status){
method - http://twitter.com/statuses/update.xml?status=;
encoded_status - URLencode(status)
request - paste(method,encoded_status,sep = )
postForm(request,.opts = opts)
}
See also:
Hi
Did you try to issue
library(MASS)
before calling isoMDS?
Regards
Petr
r-help-boun...@r-project.org napsal dne 10.02.2009 10:22:57:
Dear R-Help-Team,
I tried to use isoMDS-Function of the MASS-Package:
Message:
Fehler: konnte Funktion isoMDS nicht finden
(error: could not find
Hi,
here is a practical problem we would like to solve. In a pneumatic post the
acceleration of the capsule is measured and plotted over time. From the
graph achieved we would like to derive some kind of statistic value that
describes the stress the capsule, or what is in it, is exhibited to.
Hello,
suppose I have a data frame:
mat
id age
1 NA NA
2 NA NA
3 NA NA
4 NA NA
5 NA NA
Then I attach the data frame:
attach(mat)
I assign some new values...
id - sample(100,5)
age - rnorm(5,mean=30)
Then I want to create a new data frame from the variables id and age which
still
Struggled a whole day, but still could not get DCOM working on my machine.
What I did is:
1. Installed R 2.8.1 under folder ../R/R-2.8.1
2. Installed R_Scilab_DCOM3.0-1B5 under folder .. \R\(D)COM Server
3. Installed packages: rscproxy_1.0-12, RExcelInstaller_3.0-10,
rcom_2.0-4...
4. run the
znmeb wrote:
On Mon, Feb 9, 2009 at 4:51 AM, Neil Shephard nsheph...@gmail.com wrote:
There are pro's and con's to each of the GNU/Linux flavours and its
really a
matter of deciding which you like/have invested time in learning.
Irrespective its still simple to install R from source
Hi
r-help-boun...@r-project.org napsal dne 10.02.2009 12:02:53:
Hello,
suppose I have a data frame:
mat
id age
1 NA NA
2 NA NA
3 NA NA
4 NA NA
5 NA NA
Then I attach the data frame:
attach(mat)
Look into docs what attach does. If you do not understand environments use
Also, try these alternatives:
sub(., , x, fixed = TRUE)
chartr(., , x)
On Tue, Feb 10, 2009 at 4:02 AM, Sergey Goriatchev serg...@gmail.com wrote:
OK, got it now:
sub([.], , ED4.Comdty)
On Tue, Feb 10, 2009 at 10:01, Sergey Goriatchev serg...@gmail.com wrote:
Hello, everyone
How do I
Jesús Guillermo Andrade wrote:
Dear fellows: This is the problem: I have 5 variables A, B, C, D and E
with a range from 1 to 100 with 0.1 steps. Depending on the different
values these have, the results of the formula change:
alitemp - ((Abase/llmcc$Clase)*PClase)+(((1/llmcc
Krehbiel, Dwight wrote:
I am getting the following error on iMac running Mac OS 10.5.6; this occurs
whenever I try to import a text file into R Commander:
Error in splitCmd(command) : unbalanced quotes
So can you tell us exactly what you did and what the resulting code in
the code
Hi ,
What I did was similar, install R, R(D)COM server an RExcelInstaller package.
Then just write installRExcel() in the console. And it works.
Something you should know is that there is a specific mailing list for RCOM
topic in here:
attach provides a copy of rather than aliases to the variables within a
data frame.
d = data.frame(x=0)
attach(d)
x
# 0, from the attached copy of d
x = 1
x
# 1, from the global anvironment
d$x
# 0, from d
x - 2
x
# 1, from the global environment
d$x
# 0, from d
get('x', pos=2)
# 2, from the
Dear R Users,
I have a CO2 time series. I want to fit this series seasonal cycle and trend
with fourth harmonic function,
and then compute residuals.
I am doing something like:
file-read.csv(co2data.csv)
names(file)
attach(file)
Thanks a lot everyone. You've been a great help!
--
View this message in context:
http://www.nabble.com/Choosing-a-random-number-between-x-and-y-tp21914106p21930844.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org
Well, I knew that attach() only creates a copy of the variables in the
search path.
What I wanted to ask was how to *retrieve* that copy...
mat
id age
1 NA NA
2 NA NA
3 NA NA
4 NA NA
5 NA NA
attach(mat)
id - sample(100,5)
age - rnorm(5,mean=30)
How can I make a new data frame
Thank you, Patricia,
So basically, both of us did the same thing, but installRExcel() just
worked in your case. Because I was installing (D)COM on my company's
computer ( I have temporary admin right), do you think it may be due to
that Excel won't allow the change caused by calling
dear Yogesh
It appears that your model based on parametric terms is too inflexible..
A better alternative to parametric harmonic terms is a spline-based
approach, may be cyclic splines.. Have a look to the mgcv package..
vito
Yogesh Tiwari ha scritto:
Dear R Users,
I have a CO2 time
Thanks, Gabor.
You are always very helpful.
On Tue, Feb 10, 2009 at 12:42, Gabor Grothendieck
ggrothendi...@gmail.com wrote:
Also, try these alternatives:
sub(., , x, fixed = TRUE)
chartr(., , x)
On Tue, Feb 10, 2009 at 4:02 AM, Sergey Goriatchev serg...@gmail.com wrote:
OK, got it now:
David Croll wrote:
Well, I knew that attach() only creates a copy of the variables in the
search path.
What I wanted to ask was how to *retrieve* that copy...
mat
id age
1 NA NA
2 NA NA
3 NA NA
4 NA NA
5 NA NA
attach(mat)
id - sample(100,5)
age - rnorm(5,mean=30)
How
Yacas was completely rewritten in java (Ryacas interfaces to the
C version) since the Ryacas project started so I would not exactly
characterize yacas as dead. The work that is going on in yacas
may not have high visibility but that does not mean there is none.
Also while Maxima is more
G'day David,
On Tue, 10 Feb 2009 13:23:50 +0100
David Croll david.cr...@gmx.ch wrote:
Well, I knew that attach() only creates a copy of the variables in
the search path.
What I wanted to ask was how to *retrieve* that copy...
mat
id age
1 NA NA
2 NA NA
3 NA NA
4 NA NA
5 NA
Hi all,
Recommendation technique is popular in recent web application. I was
wondering is there any R package used for recommender system?
Thanks in advance.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Berwin A Turlach wrote:
G'day David,
On Tue, 10 Feb 2009 13:23:50 +0100
David Croll david.cr...@gmx.ch wrote:
Well, I knew that attach() only creates a copy of the variables in
the search path.
What I wanted to ask was how to *retrieve* that copy...
mat
id age
1 NA NA
2 NA NA
Berwin A Turlach wrote:
G'day David,
On Tue, 10 Feb 2009 13:23:50 +0100
David Croll david.cr...@gmx.ch wrote:
Well, I knew that attach() only creates a copy of the variables in
the search path.
What I wanted to ask was how to *retrieve* that copy...
mat
id age
1 NA NA
2 NA
Hi
r-help-boun...@r-project.org napsal dne 10.02.2009 13:57:44:
G'day David,
On Tue, 10 Feb 2009 13:23:50 +0100
David Croll david.cr...@gmx.ch wrote:
Well, I knew that attach() only creates a copy of the variables in
the search path.
What I wanted to ask was how to *retrieve*
Dear Uwe,
Thanks for trying to field this question, but the error is unlikely to be
related to the data since splitCmd() is used internally to format commands
generated by the Rcmdr so that they fit in the Script window. Unfortunately,
the contents of the Script window are likely to be
botto b.otto at uke.uni-hamburg.de writes:
here is a practical problem we would like to solve. In a pneumatic post the
acceleration of the capsule is measured and plotted over time. From the
graph achieved we would like to derive some kind of statistic value that
describes the stress the
A powerful scheme for harmonic inversion of time signals known as
filter diagonalization method is available from MIT: http://ab-initio.mit.edu/wiki/index.php/Harminv
I don't know of any R interface, but it might be a good option for
your problem.
Cheers,
baptiste
On 10 Feb 2009, at
I am trying to plot two things on one graph, with a y-axis at the left for
one variable, and a y-axis at the right, for the other variable. I thought
I could do as follows
par(mar=rep(4.5, 4))
plot(0:10, 10:20)
par(new=TRUE)
plot(0:10, 100:110, axes=FALSE)
axis(4)
but this writes both y-axis
Try something like
rm(list=ls(pattern=y+))
John
2009/2/10 stephen sefick ssef...@gmail.com:
It depends on how the data is set up (I am not an expert), but I have
had good results with the subset function. subset(x, var!=3 var!=4)
this will take the subset of the dataframe x where var is
My main problem is that yacas can only factorize polynomials in one
variable.
Any CAS which is not able to factor a^2 - b^2 into (a+b)*(a-B)
for me is less than satisfactory.
On Feb 10, 2009, at 1:51 PM, Gabor Grothendieck wrote:
Yacas was completely rewritten in java (Ryacas interfaces to
Hi
I am actually running the White test for correcting Heteroscedasticity. I
used sandwich() car(), however the output shows the updated t test of
coefficients, with revised Standard Errors, however the estimates remained
same. My problem is that the residuals formed a pattern in the original
Dan Kelley wrote:
I am trying to plot two things on one graph, with a y-axis at the left for
one variable, and a y-axis at the right, for the other variable. I thought
I could do as follows
par(mar=rep(4.5, 4))
plot(0:10, 10:20)
par(new=TRUE)
plot(0:10, 100:110, axes=FALSE)
axis(4)
but this
Hi R-Users,
I have a time series of bivariate observations (x,y,t). I've
developed a few panel routines to explore certain aspects of the data,
and now would like to break the time series up into smaller chunks to
explore the temporal trends. I would like to know if anyone has any
Thank you, this does the job perfectly!
Regards, David
R new.df - data.frame(id,age)
Or, perhaps more automatic:
R new.df - do.call(data.frame, sapply(ls(pos=2), as.name))
Is this what you are after?
Cheers,
Berwin
__
Generally in the presence of heteroskedasticity of unknown form OLS
produces consistent estimates of your regression coefficients. The
estimates of standard errors are biased in the presence of
heteroskedasticity, White's procedure is a way of producing
consistent estimates of the standard
Hi all,
I am using ISOdatetime, and I just found out that when I do
ISOdatetime(1995,03,26,2,10,0) (or any other value under minutes or seconds)
I get NA
This does not happen with the same time in other dates, nor with different
hours in the same date.
Any hint why this happens?
Using R 2.8.1
Dear Kishore,
Yes, White's heteroscedasticity-consistent standard errors are just that --
standard errors for the OLS coefficients that are consistent in the presence
of heteroscedasticity. The coefficients themselves don't change. There is an
issue here: although the standard errors and OLS
Hello everybody,
I have a problem with the the function PVCM in the plm package:
I have an unbalanced panel with over 200 companies and 1 to 14 years for
each company. I already excluded all NAs:
dataset-na.omit(dataset)
Now I'm trying to fit a variable coefficients model:
model_vc -
pbarros wrote:
Hi all,
I am using ISOdatetime, and I just found out that when I do
ISOdatetime(1995,03,26,2,10,0) (or any other value under minutes or seconds)
I get NA
This does not happen with the same time in other dates, nor with different
hours in the same date.
Any hint why this
On Tue, 10 Feb 2009, pbarros wrote:
Hi all,
I am using ISOdatetime, and I just found out that when I do
ISOdatetime(1995,03,26,2,10,0) (or any other value under minutes or seconds)
I get NA
This does not happen with the same time in other dates, nor with different
hours in the same date.
Any
I tried using the command-line and indeed it does not buffer so that,
for example:
for (i in 1:5) {
print(x)
Sys.sleep(3) }
does as expected.
However, the following sleeps for 15 seconds and then produces no
output:
library(quantmod)
for (i in 1:5) {
getQuote(AAPL)
Sys.sleep(3) }
Why
Duncan Murdoch-2 wrote:
...Are you really sure you want a graph with two different scales on it?
Thanks, Duncan. I'll use mtext. It is common for people in my research area
to put several axes on plots. The top two panels from
library(oce)
data(ctd)
plot(ctd)
shows this. Actually,
Oops...forgot to print(getQuote(AAPL). Sorry.
On Feb 10, 2009, at 10:21 AM, Fuchs Ira wrote:
I tried using the command-line and indeed it does not buffer so
that, for example:
for (i in 1:5) {
print(x)
Sys.sleep(3) }
does as expected.
However, the following sleeps for 15 seconds and
Hello,
to plot the result of a singular non linear regression (using nls) I usually
use the function plotfit, for example:
r.PTG.V-nls(PTG.P~ fz1(Portata, a,b), data=dati, start=list(a=10, b=10),
nls.control(maxiter=200), algorithm='port',
trace=TRUE, na.action=na.omit, lower=list(a=0, b=10),
Dear All,
is there a way to superimpose points and/or lines on a surface plot?
Below I try to explain my problem.
Suppose I have the following surface plot (likelikood for the normal
variable when both parameters are unknown):
Hello,
I'm using the following for loop to find regression curves using a list of
functions (formList), a list of starting
values (startList), uppervalues (upperList) and lower values (lowerList).
A sample of the list of function I use in the loop is the following:
FormList - list(PTG.P ~
Thanks. Now I know how to handle it. It was killing me
Pedro
Prof Brian Ripley wrote:
On Tue, 10 Feb 2009, pbarros wrote:
Hi all,
I am using ISOdatetime, and I just found out that when I do
ISOdatetime(1995,03,26,2,10,0) (or any other value under minutes or
seconds)
I get NA
Hi All,
I am watching a strange behaviour of ISOdatetime. In my work
computer, I get NA when I try to do
ISOdatetime(1995,03,26,2,0,0)
[1] NA
But on other dates and/or times (hour) works OK
ISOdatetime(1995,03,25,2,0,0)
[1] 1995-03-25 02:00:00 GMT
In my home computer, I do not have this
Thanks a lot. I will need to change the time to GMT, then, I guess (this is
indeed GMT time).
Pedro
Daylight savings time, I guess. The clock went directly from 1:59 to
3:00 that day. Other systems behave strangely (er, handle the anomaly)
in different ways:
ISOdatetime(1995,03,26,1,10,0)
My guess is that your two machines have different setting for locale.
What does this produce on each of them:
sessionInfo()$locale
(Note: It would have been part of the information that you were asked
to provide per the posting guide.)
---
For questions about unexpected behavior or a
See ?persp, particularly the example that starts with
# (2) Add to existing persp plot - using trans3d() :
Uwe Ligges
Domenico Vistocco wrote:
Dear All,
is there a way to superimpose points and/or lines on a surface plot?
Below I try to explain my problem.
Suppose I have the following
Hi,
I have scenario in which I wish to check whether numeric values in one array
falls within the range of numbers defined over corresponding values in two
other vectors:
starts = c(12,45,67,110)
ends=c(24, 58,102,150)
trgroup=c(18,87)
The result should be 1,3 , indices of vector starts/ends.
lauramorgana at bluewin.ch lauramorgana at bluewin.ch writes:
Hello,
I'm using the following for loop to find regression curves using a list of
functions (formList), a list of
.. long non-reproducible code removed
And I got the following error message:
Error in chol2inv(object$m$Rmat()) :
?predict
On Tue, Feb 10, 2009 at 11:06 AM, lauramorg...@bluewin.ch
lauramorg...@bluewin.ch wrote:
Hello,
to plot the result of a singular non linear regression (using nls) I usually
use the function plotfit, for example:
r.PTG.V-nls(PTG.P~ fz1(Portata, a,b), data=dati, start=list(a=10,
or simultaneously estimate the coefficients and variance structure via
nlme::gls and its 'weights' argument...
On Tue, Feb 10, 2009 at 7:57 AM, John Fox j...@mcmaster.ca wrote:
Dear Kishore,
Yes, White's heteroscedasticity-consistent standard errors are just that --
standard errors for the
on 02/10/2009 11:02 AM T Joshi wrote:
Hi,
I have scenario in which I wish to check whether numeric values in one array
falls within the range of numbers defined over corresponding values in two
other vectors:
starts = c(12,45,67,110)
ends=c(24, 58,102,150)
trgroup=c(18,87)
The result
Nobody can help with this question (or tell me where I can find help) ?
Antje schrieb:
Hi there,
I'm not sure whether this is the right mailing list to put this
question... if not, please tell me where I can find help.
I have a spare dual-core machine which run Windows XP Professional and
sapply(trgroup, function(x) which(x=starts x =ends))
works for your example data, but there's probably a better way
HTH
KJ
T Joshi tejalonl...@gmail.com wrote in message
news:11417a880902100902w53664a3dq11aee64a963d8...@mail.gmail.com...
Hi,
I have scenario in which I wish to check whether
In this case, following your instructions and a bit of noodling with
the example solves the problem with:
persp(teta1, teta2, matrixMuSigma) - res
points(trans3d( teta1, sigmaMax, apply(matrixMuSigma, 1, max),
pmat=res), col=red)
--
David Winsemius
On Feb 10, 2009, at 11:50 AM, Uwe Ligges
On Tue, 10 Feb 2009, Keith Jewell wrote:
sapply(trgroup, function(x) which(x=starts x =ends))
works for your example data, but there's probably a better way
For large scale problems this works well:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/60815.html
HTH,
Chuck
HTH
Antje wrote:
Nobody can help with this question (or tell me where I can find help) ?
Antje schrieb:
Hi there,
I'm not sure whether this is the right mailing list to put this
question... if not, please tell me where I can find help.
I have a spare dual-core machine which run Windows XP
Did you try the forums on Knime? That might be a better bet.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Antje
Sent: Tuesday, February 10, 2009 12:24 PM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] R running as a server on
I, for one, don't fully understand the question, and while there are many on
this list that are much smarter than I am, they may also have trouble
understanding the question.
At first glance it looks like you are asking if R runs on Windows XP (the
answer is yes, but this seems to simple to be
All -
I am evaluating an arima(2,1,3) and arima(3,1,3) and notice the
log-likelihood of the restricted model is higher than the log-likelihood
of the unrestricted.
Since these are nested models, I thought the unrestricted model would have
a log-likelihood at least as large as that of the
Have you checked that that time exists in the time zone you are using? From
?ISOdatetime:
Note
...
Remember that in most timezones some times do not occur and some occur twice because of transitions to/from summer time. What happens in those cases is OS-specific.
You could try working out
I've been trying to use Rterm.exe, mainly because of Rgui.exe's lack of
support for incremental history searching. But on my system (Win XP, SP3),
it (rterm from R-2.8.1) frequently just stops responding to keyboard input.
This happens often (enough to make it unusable), and generally after a
Harsh wrote:
At the cost of sounding far-fetched and almost incredulous, I would
like to know if any R user is remotely considering the use of R on
Mobile devices, and Android in particular.
In addition to the other objections in the other replies, I add this
one: the assignment operator and
Yes you can run the Rserve application on a remote computer in e.g network
and then you can connect
it from your Desktop computer. I've done this already with a remote Linux
computer (Rserve is running
very stable on a Linux computer).
For example i clustered an R,G,B image on a linux computer
Hi,
this sounds awfully similar to what I reported in r-help thread '[R]
Windows Vista, Rterm LeftAlt + Tab issue' on Aug 21, 2008, cf.
https://stat.ethz.ch/pipermail/r-help/2008-August/171548.html
At least we are not alone ;) ...and I though it was an issue with
Windows Vista only, because
Hi ,
Does any one know how to split a character vector , I have a vector X that
looks like this and each row has 3 characters
X
ASK
DGH
ASG
AUJ
FRT
I would like to split the vector into 3 vectors that look like this
X1 X2 X3
A S K
D G H
A S G
A
Dear R users,
Is there any way to control the size of the box around the mean when
creating a Forest plot using the forestplot function?
regards
Mark
Mark T. Marino, MD
VP, Early Clinical Development
Mannkind Corp.
61 S. Paramus Road
Paramus, NJ 07652
201-983-5238 Office
Dear kayj,
Try this:
X-c(ASK, DGH, ASG, AUJ, FRT)
res-data.frame(do.call(rbind,strsplit(X,)))
colnames(res)-paste(X,1:ncol(res),sep=)
res
HTH,
Jorge
On Tue, Feb 10, 2009 at 1:50 PM, kayj kjaj...@yahoo.com wrote:
Hi ,
Does any one know how to split a character vector , I have a vector X
See ?strsplit
Hope this is helpful,
Dan
Daniel J. Nordlund
Washington State Department of Social and Health Services
Planning, Performance, and Accountability
Research and Data Analysis Division
Olympia, WA 98504-5204
-Original Message-
From: r-help-boun...@r-project.org
Nguyen Dinh Nguyen n.nguyen at garvan.org.au writes:
Dear Christian,
Every single time check update package, “cwhmisc” always requests updating.
I’m aware that the package was latest updated in CRAN on 20Nov2008.
Is there anything wrong with my R library or somethingelse?
I use R 2.8.1 on
kayj wrote:
Hi ,
Does any one know how to split a character vector , I have a vector X that
looks like this and each row has 3 characters
X
ASK
DGH
ASG
AUJ
FRT
I would like to split the vector into 3 vectors that look like this
X1X2 X3
A S K
D G
Another option is:
read.table(textConnection(gsub((\\w), \\1;, X)), sep = ;)
On Tue, Feb 10, 2009 at 4:50 PM, kayj kjaj...@yahoo.com wrote:
Hi ,
Does any one know how to split a character vector , I have a vector X that
looks like this and each row has 3 characters
X
ASK
DGH
ASG
hi: maybe it's minimizing the negative of the likeihood rather than
maximizing the likelihood ?
other than that, i don't see that being possible because of what you
said. see what
happens if you estimate an arima(1,1,3) ?
On Tue, Feb 10, 2009 at 12:51 PM, Stephen Collins wrote:
All -
I
Henrik,
Ah, I did briefly search the list archive, but hadn't come across your
earlier post, sorry.
I can confirm the conditions you outlined back then (ie. it's L-ALT+TAB, and
only that, that causes the problem, that on return to rterm I can use ctrl-C
to interrupt a running R call but then no
Hi,
I have this very simple two-line script :
test.r:
print (abc);
write(def,test.txt);
When run with rscript.exe --vanilla test.r (or just rscript.exe test.r), I
get the abc output in the shell but no file is created in the
directory
I noticed that the only time the test.txt
OK,
obviously that was a silly question...
SOLUTION: just specify the path in the write function...!
Sorry !
Applejus wrote:
Hi,
I have this very simple two-line script :
test.r:
print (abc);
write(def,test.txt);
When run with rscript.exe --vanilla test.r (or just
Marino, Mark wrote:
Dear R users,
Is there any way to control the size of the box around the mean when
creating a Forest plot using the forestplot function?
H,
forestplot
Error: object 'forestplot' not found
PLEASE do read the posting guide
Hi all,
I have data from an experiment with 3 independent variables, 2 are
within and 1 is between. In addition to the dependent variable, I have
a covariate that is a single measure per subject. Below I provide an
example generated data set and my approach to implementing the ANCOVA.
However the
Dear everybody!
Need help with graphics. I am runnig a simple lm and then using
allEffects from 'effects' package:
require(effects)
model-lm(Y~A+B, data=mydataframe)
I am trying to build (for each predictor - A and then B) a plot of means on Y.
I was successful doing it like this - in one swoop:
I am appealing to the general collective wisdom of this
list in respect of a statistics (rather than R) question. This question
comes to me from a friend who is a veterinary oncologist. In a study
that
she is writing up there were 73 cats who were treated with a drug called
piroxicam. None
Hi everyone,
We're very pleased to offer a two-day looking-at-data course.
July 30-31
Washington DC
Day one: static graphics with ggplot2
Day two: interactive graphics with rggobi and GGobi.
You can attend one day (for $295) or both days (for $550).
Student discounts are available.
All
In the biomedical arena, at least as I learned from Rosner's
introductory text, the usual approach to analyzing paired 2 x 2 tables
is McNemar's test.
?mcnemar.test
mcnemar.test(matrix(c(73,0,61,12),2,2))
McNemar's Chi-squared test with continuity correction
data: matrix(c(73,
I have a problem with ifelse(), I do not understand how it works.
X-c(2,2,1,1,0,0)
str(X)
num [1:6] 2 2 1 1 0 0
Y-ifelse(X0,1,0)
Y
[1] 1 1 1 1 0 0
Can some one explain what is going on, I do not understand what ifelse is
doing in this case. Can someone explain the output Y.
Thanks
--
I have a problem with ifelse(), I do not understand how it works.
X-c(2,2,1,1,0,0)
str(X)
num [1:6] 2 2 1 1 0 0
Y-ifelse(X0,1,0)
Y
[1] 1 1 1 1 0 0
Can some one explain what is going on, I do not understand what ifelse is
doing in this case. Can someone explain the output Y.
Thanks
--
Dear Colleagues,
When I try to install Rcmdr the following happens:
--
library(Rcmdr)
Error in structure(.External(dotTclObjv, objv, PACKAGE = tcltk), class =
tclObj) :
[tcl] invalid command name font.
Error : .onAttach failed in 'attachNamespace'
Error:
It is checking the series of values of X and returning a series of 1
and 0's;
1 if it meets your condition (X0) and
0 if it doesn't.
What did you expect?
A more complex invocation might look like this:
ifelse( X 0, 1:3, -1:-2)
[1] 1 2 3 1 -1 -2
Note the recycling of the elements
Hello,
ifelse checks condition whether X (or, ina fact every value in X) is greater
than 0. If it is 1 is assigned and if it isn't then 0.
2 and 1 are greater than 0, therefore the first four values of Y are 1. 0
isn't greater than 0, therefore the last two values of Y are 0.
2009/2/10 kayj
1 - 100 of 154 matches
Mail list logo