Hello John, thank you for the fast reply.
Thanks to your answer I was able to reproduce the "between" residuals by
simply writing:
sum(residuals(mod.ok)^2)
But I must admit that how to obtain the "within" was beyond me, so some
advice here would be of great help.
(Also, it might be worth adding t
Nicole Hackman wrote:
>
> Hello, I have a very simple data set i imported from excel including 96
> averages in a column along with 96 standard errors associated with those
> averages (calculated in excel). I plotted the 95 averages using r and I
> am
> wondering if it is possible to plot the s
Abelian-2 wrote:
>
> Dear all
> when i draw some figures by R, i try to use "jpeg" to display my
> result.
> However, it isn't still clear as i want.
> Especially when i compare the origial figures by R and the ".jpeg"
> file
>
JPEG is designed for photos (the "P" stands for "Photographic"). T
Hi
r-help-boun...@r-project.org napsal dne 19.02.2009 07:15:00:
> Hi all,
> I got one problem with using %in% condition while using in if()
condition
> where
> I used the condition as
>
>
> if(SubFinSpt$SPECIMENTYP %in% CAP$SPECIMENTYP)
SubFinSpt$SPECIMENTYP %in% CAP$SPECIMENTYP
gives you lo
Petr Pikal
petr.pi...@precheza.cz
724008364, 581252140, 581252257
r-help-boun...@r-project.org napsal dne 19.02.2009 07:47:40:
> Hi,
>
> I have a vector that look like this:
>
> > dat
> V1 V2 V3 V4 V5 V6
> 0.00377467 0.00377467 0.00377467 0
Hi Gundala: Below works but I was trying to figure out a way to not have
to add the last line and I gave up.
dat<-c(0.00377467,0.00377467,0.00377467,0.00380083,0.00380083,0.00380083,0.00380959,
0.00380959,0.00380959,0.00380083,0.00380083,0.00380083)
# MAKE A TEMPORARY MATRIX
temp <- matrix(dat
Hi,
I have a vector that look like this:
> dat
V1 V2 V3 V4 V5 V6
0.00377467 0.00377467 0.00377467 0.00380083 0.00380083 0.00380083
V7 V8 V9V10V11V12
0.00380959 0.00380959 0.00380959 0.00380
Hi all,
I got one problem with using %in% condition while using in if() condition
where
I used the condition as
if(SubFinSpt$SPECIMENTYP %in% CAP$SPECIMENTYP)
this "if()"condition is in "else" condition
and hear "*SubFinSpt$SPECIMENTYP*" having only one value but
"*CAP$SPECIMEN
I found the culprit. I had a number of variables in the SPSS file that
were a variable length string data type (255 characters). This seemed
to force R into creating 255-byte variables which eventually choked my
machine's memory...
On Feb 18, 5:34 pm, Uwe Ligges
wrote:
> dobomodewrote:
> > Hello
Hello
Can anyone tell me what I am doing wrong below? My Y and y_hat are the same.
A<-scale(stackloss)
n1<- dim(A)[1];n2<-dim(A)[2]
X<-svd(A)
Y<- matrix(A[,"stack.loss"],nrow=n1)
Y
y_hat <-matrix((X$u%*% t(X$u))%*%Y,nrow=n1,byrow=T)
y_hat
[[alternative HTML version deleted]]
__
Dear all
when i draw some figures by R, i try to use "jpeg" to display my
result.
However, it isn't still clear as i want.
Especially when i compare the origial figures by R and the ".jpeg"
file
the code is follwed
jpeg(filename=,width = 1200, height = 800, pointsize = 14)
..
dev.o
I went to this webpage (http://ubuntuforums.org/showthread.php?t=624644) and
followed the instructions to the letter on getting the latest Wine. I
installed WinBUGS again, but this time I cannot open it in Wine. It says
Black Box, Trap #101, and some text I can't copy/paste here. Is this the
la
Sorry. This is still unclear to me.
I generated a file called "Test.R" that contains the following lines:
commandArgs(TRUE)
cat("\n A = ",A,"\n")
cat("\n B = ",B,"\n")
cat("\n C = ",C,"\n")
First of all I have to clarify which command line we are talking about.
If I run the command "Rscript" from
Tobias, this looks great! I'm new to Linux, and I've only installed packages
from within R with install.packages(). I'm guessing that I should unzip
ROpenBUGS to the following path:
/usr/local/lib/R/site-library/
...just checking because I'm new and don't want to do it wrong. Thanks.
PS
This
Hello folks,
I am very excited to have discovered R and have been exploring its
capabilities. R's regression models are of great interest to me as my
company is in the business of running thousands of linear regressions
on large datasets.
I am using biglm to run linear regressions on datasets tha
Hello,
I know that two possible approaches to dealing with clustered data would be
GEE or a robust cluster covariance matrix from a standard regression. What
are the differences between these two methods, or are they doing the same
thing? Thanks.
--
View this message in context:
http://www.nabb
I realise that in the case of loglin the parameters are clacluated post
festum from the cell frequencies,
however other programmes that use Newton-Raphson as opposed to IPF work the
other way round, right?
In which case one would expect the output of parameters to be limited to the
particular cont
Try setting exclude=NULL if you think that the NA's need to be considered as a
value.
The default is to exclude both NaN and NA.
--
David Winsemius
-- Original message --
From: "Farley, Robert"
> I presume the behavior below (no Connector values in the second
-- Original message --
From: Dimitri Liakhovitski
> Thank you for providing advice on this graphics question.
>
> I am building an interaction.plot.
>
> d=data.frame(xx=c(3,3,2,2,1,1),yy=c(4,3,4,3,4,3),zz=c(5.1,4.2,4.4,3.5,3.3,-1.1,-
> 1.3)
You have different l
Dear list,
After thinking about it a little more I solved my question of why I was
calculating different residuals when using
residuals.lm(...,type="partial") and when dropping a single term and
recalculating the residuals. This is because the two variables are in a
sense competing with one
Check out ggplot2:
http://had.co.nz/ggplot2/
Particularly:
http://had.co.nz/ggplot2/geom_errorbar.html
If you need more help you'll have to provide a self-contained set of data/code.
Hopefully you'll soon be completely rid of this strange "excel" of
which you speak :Op
On Wed, Feb 18, 2009 at 9
It returns a chi-squared statistic with one degree of freedom.
--
David Winsemius
-- Original message --
From: Timthy Chang
>
>
> >See the cpower() and spower() functions in Frank Harrell's Hmisc package
> >on CRAN.
> >
> >HTH,
> >
> >Marc Schwartz
>
> How to
On Wed, 18 Feb 2009, maiya wrote:
I am fairly new to log-linear modelling, so as opposed to trying to fit
modells, I am still trying to figure out how it actually works - hence I am
looking at the interpretation of parameters. Now it seems most people skip
this part and go directly to measuring
>See the cpower() and spower() functions in Frank Harrell's Hmisc package
>on CRAN.
>
>HTH,
>
>Marc Schwartz
How to calculate the p-value of Log rank test ??
Thank you very much!
--
View this message in context:
http://www.nabble.com/Log-rank-test-power-calculations-tp15208266p22092630.html
Se
On Wed, 18 Feb 2009, Jong wrote:
Hi all,
I've compiled R (version 2.8.1) from the source code with "--with-
blas=" option in order to use external multi-threaded blas package.
However, it's very hard to verify that compilation is correct and R is
using multi-threads correctly or not. Which comm
I presume the behavior below (no Connector values in the second table)
is because I have na values in my dataframe. How do I get them to list?
I found references to na.action and na.pass, but I could not figure out
how I could use them with xtabs to list the na as a value. I'm sure
there's a
Dear list,
I would like to know how the function residuals.lm calculates the
partial residuals from an lm object with more than one predictor
variable. In other words what is residuals.lm(...,type="partial") doing
behind the scenes? According to the help file for residuals.lm
(?residuals.lm
Hello, I have a very simple data set i imported from excel including 96
averages in a column along with 96 standard errors associated with those
averages (calculated in excel). I plotted the 95 averages using r and I am
wondering if it is possible to plot the second column of standard errors
while
Hello,
I am attempting to "automate" a Bayesian normal linear regression
using Zelig. Basically, I have a list containing several zelig()
objects, each having a different formula, same data set, and same
model (normal.bayes).
My problem lies in the setx() method, where I am setting a numeric
par
On Wed, Feb 18, 2009 at 4:21 PM, Felipe Carrillo
wrote:
> Hi all:
> Using the example below, is there a way to add Y axis titles to each graphic
> instead of sharing the same title?
Not at the moment, no. It's on my to do list.
Hadley
>
> library(ggplot2)
>
> RT = matrix(c(814, 500, 424, 394,
Try this (or replace selected occurrences of * with ~ to get spaces
instead of juxtaposition):
plot(1:3, main = ~ "R=[" * 0.001 * mu * Mu * "]:A=[" * 750 * mu * Mu * "]")
On Wed, Feb 18, 2009 at 3:56 PM, gina patel wrote:
> I would like to add the greek letter mu to replace u in my title shown b
Stavros Macrakis wrote:
> On Thu, Feb 12, 2009 at 4:28 AM, Gavin Simpson
> wrote:
>
>> When I'm testing the speed of things like this (that are in and of themselves
>> very quick) for situations where it may matter, I wrap the function call in
>> a call
>> to replicate():
>>
>> system.time(re
James,
you're probably interested in image function rather than in
filled.contour. Type ?image to see the syntax.
Cheers,
Patrizio
2009/2/18 James Nicolson :
> Hi,
>
> Thanks for your help. I have looked at the beginners documentation and
> while there are options to configure various aspects of
I am fairly new to log-linear modelling, so as opposed to trying to fit
modells, I am still trying to figure out how it actually works - hence I am
looking at the interpretation of parameters. Now it seems most people skip
this part and go directly to measuring model fit, so I am finding very few
dobomode wrote:
Hello R-help,
I am trying to import a large dataset from SPSS into R. The SPSS file
is in .SAV format and is about 1GB in size. I use read.spss to import
the file and get an error saying that I have run out of memory. I am
on a MAC OS X 10.5 system with 4GB of RAM. Monitoring t
Hi all:
Using the example below, is there a way to add Y axis titles to each graphic
instead of sharing the same title?
library(ggplot2)
RT = matrix(c(814, 500, 424, 394, 967, 574, 472, 446),4,2)
colnames(RT) = c('repetition','alternation')
rownames(RT) = c('7-yrs','11-yrs','15-yrs','21-yrs')
r
Dear Tal,
I suppose that the "between" residuals would be obtained, for your example,
by residuals(mod.ok). I'm not sure what the "within" residuals are. You
could apply the transformation for each within-subject effect to the matrix
of residuals to get residuals for that effect -- is that what yo
Hi,
Thanks for your help. I have looked at the beginners documentation and
while there are options to configure various aspects of the plot none of
them seem to have the desired effect. I have managed to ensure that the
plot fills the space vertically with no margins, no axes etc (using
mai=c(0,0,
Thank you for providing advice on this graphics question.
I am building an interaction.plot.
d=data.frame(xx=c(3,3,2,2,1,1),yy=c(4,3,4,3,4,3),zz=c(5.1,4.2,4.4,3.5,3.3,-1.1,-1.3)
d[[1]]<-as.factor(d[[1]])
d[[2]]<-as.factor(d[[2]])
print(d)
interaction.plot(d$xx, d$yy, d$zz,
type="b", col=c("red
Hi all,
I've compiled R (version 2.8.1) from the source code with "--with-
blas=" option in order to use external multi-threaded blas package.
However, it's very hard to verify that compilation is correct and R is
using multi-threads correctly or not. Which command or operation in R
will run with
The call to replace is replacing the 1st 3 elements of a (your indexes in
Trans_CR) with the values and leaving the 4-6 elements alone. For what you
want, try:
A <- Trans_Prob_values[ match(a, Trans_CR) ]
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain
Alan Smith wrote:
Hello R Users and Developers,
I have a basic question about how R works. Over the past few years I have
struggled when I try to generate a new data frame that I believe should
contain numeric data in some columns and character data in others only to
find everything converted t
On 18/02/2009 4:24 PM, Christian Langkamp wrote:
Dear everyone
I would like to change values in vectors doing a translation. i.e. I have a
start vector giving me the levels in one vector (numbers 1 to x - rating)
and then I have a second vector giving me the values to be allocated (loss
probabil
I have yet another question concerning maps.
This time I want to create a colored map of number of students by state.
Can this be done using palettes in Rcolor brewer, is there some other way?
I've been looking through R archives, but am still really at a loss here, so
all help will be much apprec
The culprit is the cbind function. When given 2 vectors (not already something
else), cbind will create a matrix, not a data frame. A matrix can only have 1
type, so the numbers get converted to character. In your first example you
never do create a data frame, you just build a matrix (try st
Dear everyone
I would like to change values in vectors doing a translation. i.e. I have a
start vector giving me the levels in one vector (numbers 1 to x - rating)
and then I have a second vector giving me the values to be allocated (loss
probabilities), but the number of potential rating classes
Uwe Ligges wrote:
Hi all,
I've managed to get JAGS working on my Ubuntu Hardy Linux with a 32-bit
computer and AMD processors using R 2.8.1. JAGS is great. I've read
that
JAGS is the fastest, but that hasn't been my experience. At any rate, I
have more experience with WinBUGS under Windows
gina patel wrote:
I would like to add the greek letter mu to replace u in my title shown below.
main="R=[0.001uM]:A=[750uM]"
i tried using main=expression(R=[0.001~mu~M]:A=[750~mu~M])
plot(1:3, main=expression(paste("R=[0.001~",mu,"~M]:A=[750~",mu,"~M]")))
Ciao,
domenico
but this is not work
Hello dear R members.
I have been learning the Anova syntax in order to perform an SS type III
Anova with repeated measures designs (thank you Prof. John Fox!)
And another question came up: where/what are the (between/within) residuals
for my model?
Play code:
phase <- factor(rep
Hello R Users and Developers,
I have a basic question about how R works. Over the past few years I have
struggled when I try to generate a new data frame that I believe should
contain numeric data in some columns and character data in others only to
find everything converted to character data. Is
I would like to add the greek letter mu to replace u in my title shown below.
main="R=[0.001uM]:A=[750uM]"
i tried using main=expression(R=[0.001~mu~M]:A=[750~mu~M])
but this is not working at the moment.
any help is appreciated
thanks in advance
[[alternative HTML version d
I figured it out
thanks
On Wed, Feb 18, 2009 at 3:50 PM, Domenico Vistocco wrote:
> Alina Sheyman wrote:
>
>> I've downloaded the RcolorBrewer package, but when I try to run
>> mypalette<-brewer.pal(7,"Greens")
>>
>> (or any other command with brewer.pal)
>>
>> I get the following error message -
Alina Sheyman wrote:
I've downloaded the RcolorBrewer package, but when I try to run
mypalette<-brewer.pal(7,"Greens")
(or any other command with brewer.pal)
I get the following error message - Error: could not find function
"brewer.pal"
Does anyone know why that's happening? Is there smth els
diego Diego wrote:
Dear R experts:
I have a list (a very long one) and I need to create successively txt
outputs (on diferent files ideally) for the data of each component of the
list.
How can I do this?
Maybe this could help you:
list2Files <- list(1:3, letters[1:10], matrix(1:15, 5, 3))
I've downloaded the RcolorBrewer package, but when I try to run
mypalette<-brewer.pal(7,"Greens")
(or any other command with brewer.pal)
I get the following error message - Error: could not find function
"brewer.pal"
Does anyone know why that's happening? Is there smth else I need to
download?
Hi,
For me running winbugs through wine just works. Even when I do not specify
any directories.
The example they give in the bugs helpfile was my starting point.
Setup is suse 11.1, latest Wine, R, R2WinBUGS & winbugs.
I assume you first tried without specifying directories?
The directories yo
Harsh wrote:
Hi list,
I would like to use ggplot2 in creating a line plot with 4 lines (groups), 2
of which I want in colour and the remaining two as dotted lines.
### R code ###
library(ggplot2)
### create data
vals <- rnorm(400)
div<- c(rep("A",100),rep("B",100),rep("C",100),rep("D",100
Dear R experts:
I have a list (a very long one) and I need to create successively txt
outputs (on diferent files ideally) for the data of each component of the
list.
How can I do this?
Thanks in advance!!
[[alternative HTML version deleted]]
__
Duncan, Berwin, Martin,
Thanks for your thoughtful explanations, which make perfect sense.
May I simply suggest that the non-identity between last(cumsum) and
sum might be worth mentioning in the cumsum doc page?
-s
__
R-help@r-project.org
Dear R users,
My question is more methodology related rather than specific to R usage. Using
time on study as time in a cox model, eg:
library(Design)
stanf.cph1=cph(Surv(time, status) ~ t5+id+age, data=stanford2, surv=T)
#In this case the 1000-day survival probability would be:
stanf.surv1=su
G'day all,
On Wed, 18 Feb 2009 12:41:27 -0500
Stavros Macrakis wrote:
> Hmm. Why not use the same method to guarantee the same result?
Hmm, I did not look at the source code but, potentially, sum() could
use some tricks to reduce rounding errors further that would not be
available to cumsum(
Hi list,
I would like to use ggplot2 in creating a line plot with 4 lines (groups), 2
of which I want in colour and the remaining two as dotted lines.
### R code ###
library(ggplot2)
### create data
vals <- rnorm(400)
div<- c(rep("A",100),rep("B",100),rep("C",100),rep("D",100))
n<- rep(1:100
Hello R-help,
I am trying to import a large dataset from SPSS into R. The SPSS file
is in .SAV format and is about 1GB in size. I use read.spss to import
the file and get an error saying that I have run out of memory. I am
on a MAC OS X 10.5 system with 4GB of RAM. Monitoring the R process
tells m
Hi all,
I'm trying to fit a model using the shorthand coeff[factor] instead of
coding dummy variables. Is there a way to keep this notation when
specifying constraints? See example below:
x = runif(200)
b0 = c(rep(0,100),runif(100))
b1 = 1
fac <- as.factor(rep(c(0,1), each=100))
y = b0+b1*x+rnorm
francogrex wrote:
Hi all, do you know if an R program
If R program means pure R code, than now, R is an interpreted language.
For more information on the interpreter, you might want to start reading
the manual "R Internals".
Uwe Ligges
is compiled to machine language when
executed? And a
This is my situation:
I have a significant amount of data, and need to send it in pieces to R. I
need R to return certain parameters for further use.
I am sending files from C# (that are being queried from a database) into R.
Currently I am trying to use the R(D)-Com package to figure out how to
This is my approach:
If "cosa" is your data.frame . e.g.
>cosa
i1 i2 dis
[1,] 1 1 0.00
[2,] 1 2 0.93
[3,] 1 3 0.80
[4,] 1 4 1.00
[5,] 2 2 0.00
[6,] 2 3 0.02
[7,] 2 4 0.22
[8,] 3 3 0.00
[9,] 3 4 0.95
[10,] 4 4 0.00
# first crate a square matrix of 0
Hi all,
When running Rmdr using the demo data file using the following
commands:
data(mdrdata)
cvk<-10
nbr=2
res<-rmdr(mdrdata,10,2, randomize=TRUE)
I could not find the statistical results, but like this:
[1] 1
Cross Validation 1 Wed Feb 18 09:05:23 2009
The best set of loci is 13 17
Cross V
I am using auto.arima to find the best arima model but am a little confused
by the output. I want to choose the best model using the BIC criteria.
This is my code (straightforward where a is the data)
auto.arima(a,d=0,D=0,max.p=5,max.q=5,max.P=0,max.Q=0,max.order=5,start.p=0,start.q=0,start.P=0,st
hmm, further investigation shows that two different fits are used.
Why did nnet decide to use different fits when the data is basically
the same (2 factors in nn1 and binary in nn2)?
# uses an entropy fit (maximum conditional likelihood)
> nn1
a 57-3-1 network with 178 weights
inputs: make addres
Dear R users,
I use function kqr from package kernlab a large number of time and
every time it is used it prints an info message which slows the
process. Please note that it is not a warning message or an error
message and that there is no "info" option in the function that could
be set t
Hi all, do you know if an R program is compiled to machine language when
executed? And also is there any way to disassemble an R code/program: to see
how it is generating the machine instructions for the processor? Thanks
--
View this message in context:
http://www.nabble.com/R-code-compiled%2C-
I am using glm.nb, a ~b*c ( b is categorical and c is continuous). when I
run this model I get the warning message:
Warning messages:
1: In theta.ml(Y, mu, sum(w), w, limit = control$maxit, trace =
control$trace > :
iteration limit reached
2: In theta.ml(Y, mu, sum(w), w, limit = control$maxit,
On 18/02/2009 12:41 PM, Stavros Macrakis wrote:
Hmm. Why not use the same method to guarantee the same result? Or at
least document the possibility that cumsum(x)[length(x)] != sum(x)...
that seems like an easy trap to fall into.
Assuming equality of floating point numbers computed by two dif
For version 2.3-30 of chron which just appeared on CRAN this can
be simplified to:
library(chron)
tt <- times(0:47/48)
tt
chron("1/1/09", tt) # no rep needed
On Mon, Feb 16, 2009 at 6:04 AM, Gabor Grothendieck
wrote:
> Try this (and see R News 4/1 for more).
>
>> library(chron)
>> tt <- times(0:
To answer my own post, and for the archives (hopefully not that anyone
has to repeat what I had to do ;-), after much hair-pulling , frowning
at the screen and general dumb headedness the following slab of R code
achieves the results I wanted. It isn't elegant but does a job.
msr <- function(x) {
Can you show us the code used to get the data?
The usual methods that I can think of would have sorted the columns correctly
for you. The fact that this is not the case indicates that you are using a
different method, or doing something that looses the information along the way.
If you show u
Hmm. Why not use the same method to guarantee the same result? Or at
least document the possibility that cumsum(x)[length(x)] != sum(x)...
that seems like an easy trap to fall into.
-s
On Wed, Feb 18, 2009 at 11:39 AM, Martin Maechler
wrote:
>> "SM" == Stavros Macrakis
>>
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce
our*** R/Splus Fundamentals and Programming Techniques and R Advanced
Programming***courses at USA locations for March - April 2009.
* New York City ** March 19-20, 2009
R/Splus Fundamentals and Programming
I have a dissimilarity dataset with the form:
1 1 dissimilarity value
1 2 ...
1 3
1 4
2 2
2 3
2 4
...
I would like to do nonmetric multidimensional scaling with this data, but I
am having trouble using this format. I would like to either find a function
that accepts this format or find
> "SM" == Stavros Macrakis
> on Wed, 18 Feb 2009 10:00:40 -0500 writes:
SM> Nice! Glad to hear it. It sounds as though it is still possible for
SM> cumsum(x)[length(x)] to not be exactly equal to sum, though?
Well, possible, probably yes, platform-dependently;
However I vag
Just to be sure you're aware, there are packages for chromatograpy and mass
spec data in Bioconductor. Like xcms. Don't think any will directly address
your problem, but they might be useful.
Michael
On Tue, Feb 17, 2009 at 5:44 AM, bartjoosen wrote:
>
> Hi,
>
> I'm trying to match peaks betwee
Hi,
If the convex hull for *all* the data points is not ideal enough, is
it feasible to break the data into small subsets using clustering
methods such as kmeans() and compute the convex hull for each cluster?
Finally we are able to know the "borders" of all clusters using
chull(); I don't know ho
Thanks Prof Ripley, now I understood.
tiff(filename = "volc.tif", width=600, height=400, compression =
"none", bg = "white", res = 300, pointsize=3)
par(mar=c(3,3,2,2), cex=1)
image(t(volcano)[ncol(volcano):1,], xaxt="n", yaxt="n")
axis(1, at=seq(0, 1, 0.1), cex.axis=0.8, tick=T)
axis(2, at=seq(0,
Hi Jim,
> or sorting the columns of the above table if that is what you are using to
> plot.
How do you do that? Yes I am using that data exactly for the plotting.
- GV.
>
> Jim
>
>
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
On Wed, 18 Feb 2009, Lauri Nikkinen wrote:
Thanks. I upgraded to R 2.8.1 and tried
tiff(filename = "volc.tif", width=600, height=400, compression =
"none", bg = "white", res = 300)
image(t(volcano)[ncol(volcano):1,])
dev.off()
but this produces error
Error in plot.new() : figure margins too l
On 18/02/2009 10:00 AM, Iuri Gavronski wrote:
The funniest part is that we are close to carnival in Brazil...
vacation time... more than 30C... and I am the one to try to get back
to the track...
Any clue on how to rotate the object? Examples?
example(spin3d)
Also, still on the same object.
Lauri Nikkinen wrote:
Thanks. I upgraded to R 2.8.1 and tried
tiff(filename = "volc.tif", width=600, height=400, compression =
"none", bg = "white", res = 300)
image(t(volcano)[ncol(volcano):1,])
dev.off()
but this produces error
Error in plot.new() : figure margins too large
How this shoul
> Let Y be a normal multivariate function. For example, let Y have 4
> dimensions. I want to calculate
>
> P(Y1 < Z1, Y2 < Z2, Y3 > Z3, Y4 > Z4).
>
> There are R functions to do the calculation if all the inequalities
> are of the type "<" (the cdf). But is there an R function where the
The cdf
The funniest part is that we are close to carnival in Brazil...
vacation time... more than 30C... and I am the one to try to get back
to the track...
Any clue on how to rotate the object? Examples?
Also, still on the same object. If I change zoom, position, etc. with
the mouse, where to query thi
Thanks. I upgraded to R 2.8.1 and tried
tiff(filename = "volc.tif", width=600, height=400, compression =
"none", bg = "white", res = 300)
image(t(volcano)[ncol(volcano):1,])
dev.off()
but this produces error
Error in plot.new() : figure margins too large
How this should be modified to produce a
Thank you very much for your help
Alex
On Wed, Feb 18, 2009 at 1:26 PM, Pedro Silva wrote:
> --
>
> Message: 72
> Date: Tue, 17 Feb 2009 22:05:46 í (UTC)
> From: "Hans W. Borchers"
> Subject: Re: [R] Subset Regression Package
> To: r-h...@stat.math.ethz.ch
> Messa
Duncan Murdoch wrote:
>> There's no real difficulty there: axis takes an mgp arg as well.
Thanks for that. A good bit of practical advice, which I hadn't yet clicked
on. I won't comment on the thinking thing;)
Regards, Mark.
Duncan Murdoch-2 wrote:
>
> On 18/02/2009 7:50 AM, Mark Difford wro
Apologies, Jim Holtman has pointed out a couple of problems/queries with
my original email that I would like to make clear.
Firstly, I introduced a typo when trying to be helpful. In my email
below, I had incorrectly typed out one of the species codes I would
count:
1000
16220602
2011
240
> lm(y ~ x-1)
> solve(crossprod(x), t(x))%*%y# probably this can be done more
> efficiently
You could do
crossprod(x,y) instead of t(x))%*%y
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posti
If the data are cross-classified, then hwy would you want a hierarchical
linear model? You might try the lmer function for this instead.
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Luwis Tapiwa Diya
> Sent: Tuesday, Febru
lm does lots of computations, some of which you may never need. If speed
really matters, you might want to compute only those things you will really
use. If you only need coefficients, then using %*%, solve and crossprod will
be remarkably faster than lm
# repeating someone else's example
# lm(DAX
Gundala Viswanath wrote:
>
> Dear all,
>
> I have a binned data that looks like this:
>
>> dat
> (-1,9] (9,19] (19,29] (29,39] (39,49] (49,59] (59,69] (69,79]
> 10063374 79 1643443
> (79,89] (89,99]
>62
>
> I tri
I would think this could be approached by segmenting the probability "volume"
using identities such as these:
P(Y1 < Z1, Y2 < Z2, Y3 > Z3, Y4 > Z4) + P(Y1 < Z1, Y2 < Z2, Y3 > Z3, Y4 < Z4) =
P(Y1 < Z1, Y2 < Z2, Y3 > Z3, Y4 < Inf)
and
P(Y1 < Z1, Y2 < Z2, Y3 < Z3, Y4 Z3, Y4 < Inf) =
P(Y1 <
As it happens, I have also been looking into this. I began by
considering
Ken Clarkson's hull: http://www.netlib.org/voronoi/hull.html but
eventually
discovered that its alpha shapes don't seem to treat holes in regions,
only
simply connected regions. (I would be happy to hear to the con
1 - 100 of 134 matches
Mail list logo
