Dear R users,
in my application i can transfer images to R with the help of Rserve. The
images come from
a java application. When i plot a greyscale image (values 0-255) with images
(imageMatrix...as grey)
the image is created with inverse colours. My first question is how can i
plot the image wi
Mark Wardle a écrit :
> 1. Which version of Filemaker? NB: Framemaker is a different program
> (desktop publishing), so do be a little precise!
Dunno. The file is named "export.fm7" ; one might be tempted to infer
Filemaker 7.
"Framemaker" is a typo
> 2. If it is an ancient version, then I s
I will be out of the office starting 12/10/2007 and will not return until
15/10/2007.
Otherwise, I will reply to your email when I am back.
Have a great day!
Hwee Pin
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
2007/10/11, Ray Brownrigg <[EMAIL PROTECTED]>:
> On Thu, 11 Oct 2007, Lukas Gudmundsson wrote:
> > Hello,
> > I am trying to draw geographical maps with the maps package. However
> >
> > if I try to access the data following error occurs:
> > > require(maps)
> > > map()
> >
> > Fehler in zip.file.e
I think the problem is that tempdir() has been set somewhere invalid.
However, does help() work at all, since that uses the same mechanism?
Here are some diagnostic hints:
Start an R session, printout tempdir() and see if it looks right. Search
for the directory on your file system and see if y
I have been using random forest on a data set with 226 sites and 36
explanatory variables (continuous and categorical). When I use
"tune.randomforest" to determine the best value to use in "mtry" there
is a fairly consistent and steady decrease in MSE, with the optimum of
"mtry" usually equal t
At 22:48 12/10/07, you wrote:
>Hello,
>
>I would like to know if there is a clever way to avoid the problem
>illustrated below within the xyplot function.
>
>x <- seq(1:10)
>y <- seq(1:10)
>pr1 <- xyplot(x ~ y)
>
>u <- seq(1:12)
>v <- seq(1:12)
>pr2 <- xyplot(u ~ v, col = "red", more = FALSE)
>
>p
have you tried the function image.plot in the package fields?
On Fri, 12 Oct 2007, zhijie zhang wrote:
> Dear friends,
> Anybody has ever met the problem to add a legend to a figure generated by
> image()? I have three variables,x,y and z.
> x and y are the coordinates, and z is the third value
Still one problem remains. The list I have is a list where every element
refers to a file, where the actual matrices are. This means problems with
dealing with them. How do I turn the list of references into an actual
vector with real matrices as components?
See below for earlier details.
Svemp
Here is one way. Not sure what you wanted done with some of the other
variables, so I just chose the first one; you could do max/min:
> z <- by(h, h$BROOD, function(x){
+ # take first value of elements you don't want to change
+ data.frame(BROOD=x$BROOD[1], TICKS.mean=mean(x$TICKS),
TICKS
On 10/12/07, David Afshartous <[EMAIL PROTECTED]> wrote:
> All,
> Sorry for overly simplistic question, but I can't seem to remember how to
> create the basic plot shown in Figure 1.1 of Pinheiro & Bates (2004; p.4).
> The y-axis delineates a factor (Rail) while the x-axis displays the
> distributi
On Fri, 2007-10-12 at 11:16 -0600, D. R. Evans wrote:
> The standard chisq.test() and fisher.test() functions, when applied to
> two distributions (to determine whether the same underlying
> distribution applies to both) requires one to pre-bin the
> distributions.
>
> Is there a library function
On 10/12/2007 1:16 PM, D. R. Evans wrote:
> The standard chisq.test() and fisher.test() functions, when applied to
> two distributions (to determine whether the same underlying
> distribution applies to both) requires one to pre-bin the
> distributions.
>
> Is there a library function (either buil
Dear all,
I am developing a small suite of cgi-bin programs on my server that help
provide some things that statistical/psychometric beginners may find
helpful and which are aren't widely available. Some early e.g.s are:
http://www.psyctc.org/stats/R/multirater.html
http://www.psyctc.org/stats/R
The standard chisq.test() and fisher.test() functions, when applied to
two distributions (to determine whether the same underlying
distribution applies to both) requires one to pre-bin the
distributions.
Is there a library function (either built-in or in a package) that
acts more like the ks.test(
Desmond Campbell wrote:
>
> Dear Ben Bolker,
>
> Thanks for replying and offering advice, unfortunately it doesn't solve my
> problem.
>
> 1) The mshapiro.test() in the mvnormtest package appears only applicable
> for datasets containing 3-5000 samples, whereas my dataset contains
> 100,000
Trying to find a quick/slick/easily interpretable way to
collapse a data set.
Suppose I have a data set that looks like this:
h <- structure(list(INDEX = structure(1:6, .Label = c("1", "2", "3",
"4", "5", "6"), class = "factor"), TICKS = c(0, 0, 0, 0, 0, 3
), BROOD = structure(c(1L, 1L, 2L
Hi everyone,
anybodyâs got an idea why the following script doesnât produce
batch-histograms?
Iâm using Windows XP and R 2.5.1.
Hereâs the script:
matrix<-read.csv("C:\\Stadtwerke_mit_Umlage.csv", header=TRUE,sep=";",dec=".")
Stadtwerke<-colnames(matrix)
Bereich_blau<-66.67
Perhaps the VAR and irf functions in the vars library will allow you to do
what you want.
good luck,
spencer
On 10/12/07, Martin Ivanov <[EMAIL PROTECTED]> wrote:
>
> Dear R users,
>
> I need perform structural analysis on a no intercept VAR model.
> Unfortunately the functions irf.VAR and d
Dear R users,
I need perform structural analysis on a no intercept VAR model. Unfortunately
the functions irf.VAR and dfev that come with the MSBVAR package only work with
objects output by the reduced.form.var function, which seems to only evaluate
VAR models with intercept. Is there a way to
Hello!
I am trying to embed a plot of a curve(say x^2) on a matrix that
I am viewing using the image(matrix) command.
I was wondering if someone could give me some idea
of how to do this.
Thanks,
J.
[[alternative HTML version deleted]]
__
R-
On Friday 12 October 2007 17:31:22 David Afshartous wrote:
DA > All,
DA > Sorry for overly simplistic question, but I can't seem to remember how
to
DA > create the basic plot shown in Figure 1.1 of Pinheiro & Bates (2004;
p.4).
DA > The y-axis delineates a factor (Rail) while the x-axis display
I did ?par and it looks switching to mfrow will fix my problem so
thanks but disregard my previous message and
I apologize for not checking ?par first. I never thought it would be
that simple.
This is not an offer (or solicitation of an off
On 10/12/07, Bert Gunter <[EMAIL PROTECTED]> wrote:
> ?levelplot() is the standard lattice package version of image and
> automatically generates a legend.
>
> Note: if the image is expressed as a matrix, X, of nrow x ncol values, then
> (as in image()) then levelplot should be invoked with:
>
> r
Hello all together,
I (R Beginner) think the answer for my question is easy but anyway I
have not found the solution.
I am using hist for a frequency histogramm. But the divisions of the
xaxis is to large
How can I change the division (I am not sure if this is the right
name for what I mea
?levelplot() is the standard lattice package version of image and
automatically generates a legend.
Note: if the image is expressed as a matrix, X, of nrow x ncol values, then
(as in image()) then levelplot should be invoked with:
rowindx <- seq.int(nrow(X))
levelplot(t(X)[,rev(rowindx)],scales=
Hi,
I'm running
> sessionInfo()
R version 2.6.0 (2007-10-03)
i386-apple-darwin8.10.1
locale:
en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
I have two libraries:
(1) /Users/mk/R/i386-ap
why don't you just output to postscript instead of using dev.copy:
postscript("yourfile.ps")
...your script..
dev.off()
On 10/12/07, Thomas Schwander <[EMAIL PROTECTED]> wrote:
> Hi everyone,
>
>
>
> anybody's got an idea why the following script doesn't produce
> batch-histograms?
>
> I
On Fri, 2007-10-12 at 18:47 +0800, Samuel wrote:
> Hi,
>
> I'm quite fresh to R, and a layman of English as well. I hope I can make you
> understood.
>
> Now I have two vectors A and B. Is there any quick way to know whether B is
> a subset of A? and If B is a subset of A, can I know easily which
Samuel wrote:
> Hi,
>
> I'm quite fresh to R, and a layman of English as well. I hope I can make you
> understood.
>
> Now I have two vectors A and B. Is there any quick way to know whether B is
> a subset of A? and If B is a subset of A, can I know easily which elements
> in A (the index of A) equ
Dear friends,
Anybody has ever met the problem to add a legend to a figure generated by
image()? I have three variables,x,y and z.
x and y are the coordinates, and z is the third values. we can use image(x,
y, z,...) to generate a figure according to the z-values, but the problem is
the figure le
Hello,
I found a working but not satisfactory solution. I did compare the
file structures of "maps" and "mapdata" and found that the \data
subdirectories of both packages have one main difference:
maps\data contains a Rdata.zip file with some xxx.r and yyy.rda files.
mapdata\data does not contain
Dear all,
I am trying to understand the output from a binomial lmer object and why
the scaling of a fixed effect changes the variance components.
In the model p2rec is cbind(number recruits2,number recruits 1), Pop is
populations (five level factor) and ja is year (covariate running from
1955-
If you do something like length(coef(lm(y~.+v3:v4 + v5:v6, data=dat)))
to get a quick empirical estimate of required number of coefficients,
you will find that you have 35 coefficients, so 32 observations cannot
provide a solution at all. And indeed, nTrials=35 is the first size at
which optFederov
Hallo!
Is there a package in R that does Q-type factor analysis?
I know how to do principal component analysis, but haven't found any
application of Q-type factor analysis.
Thx,
Julia
--
Pt! Schon vom neuen GMX MultiMessenger gehört?
Der kanns mit allen: http://www.gmx.net/de/go/multimesse
Dear R users,
I need perform structural analysis on a no intercept VAR model. Unfortunately
the functions irf.VAR and dfev that come with the MSBVAR package only work with
objects output by the reduced.form.var function, which seems to only evaluate
VAR models with intercept. Is there a way to
Try
par(mar=c(3,4,2,2), mfrow=c(5,2))
On 10/12/07, Leeds, Mark (IED) <[EMAIL PROTECTED]> wrote:
> I am constructing plots ( regular not lattice ) and my initial command
> is
>
> par(mar=c(3,4,2,2), mfcol=c(5,2))
>
> and then I create 10 plots on the page. It looks great but the plots on
> the pa
All,
Sorry for overly simplistic question, but I can't seem to remember how to
create the basic plot shown in Figure 1.1 of Pinheiro & Bates (2004; p.4).
The y-axis delineates a factor (Rail) while the x-axis displays the
distribution of a continuous variable (time) according to each level of the
f
Powell, Jeff wrote:
> Hello,
>
> I would like to know if there is a clever way to avoid the problem
> illustrated below within the xyplot function.
>
> x <- seq(1:10)
> y <- seq(1:10)
> pr1 <- xyplot(x ~ y)
>
> u <- seq(1:12)
> v <- seq(1:12)
> pr2 <- xyplot(u ~ v, col = "red", more = FALSE)
>
Hi,
I'm quite fresh to R, and a layman of English as well. I hope I can make you
understood.
Now I have two vectors A and B. Is there any quick way to know whether B is
a subset of A? and If B is a subset of A, can I know easily which elements
in A (the index of A) equals to B's elements accordin
Dave,
> I have been using random forest on a data set with 226 sites and 36
> explanatory variables (continuous and categorical). When I use
> "tune.randomforest" to determine the best value to use in "mtry" there
> is a fairly consistent and steady decrease in MSE, with the optimum of
> "mtry
Hello,
I would like to know if there is a clever way to avoid the problem
illustrated below within the xyplot function.
x <- seq(1:10)
y <- seq(1:10)
pr1 <- xyplot(x ~ y)
u <- seq(1:12)
v <- seq(1:12)
pr2 <- xyplot(u ~ v, col = "red", more = FALSE)
prts <- list(pr1, pr2)
for(i in prts) pr
Colleagues,
I am analyzing data collected during oceanographic cruises. We have
conducted many cruises over the last decade. On each cruise we visit ~50
stations. At each station (termed EventNum)we lower an instrument that
measures depth, temperature, salinity and oxygen every few seconds as it
I am constructing plots ( regular not lattice ) and my initial command
is
par(mar=c(3,4,2,2), mfcol=c(5,2))
and then I create 10 plots on the page. It looks great but the plots on
the page go in the order
16
27
38
49
510
Where the numbers denote decile breakdowns.
Is there
Have you looked at layout() ?
Hadley
On 10/12/07, Leeds, Mark (IED) <[EMAIL PROTECTED]> wrote:
> I am constructing plots ( regular not lattice ) and my initial command
> is
>
> par(mar=c(3,4,2,2), mfcol=c(5,2))
>
> and then I create 10 plots on the page. It looks great but the plots on
> the page
Dear Prof. Wood,
Just another quick question. I am doing model selection following Wood
and Augustin (2002). One of the criteria for retaining a term is to see
if removing it causes an increase in the GCV score. When doing this, do
I also need to fix the smooth parameters?
Thanks,
Julian B
On 10/11/2007 6:32 PM, Deepayan Sarkar wrote:
> On 10/11/07, Karim Rahim <[EMAIL PROTECTED]> wrote:
>> Thank you for your reply.
>>
>> In this graphics context, I'm not sure what the incident or reflected
>> light rays are.
>>
>> May I ask for an example of using a colour key with the volcano data
Hey there!
I would like to justify the stability of the cluster of a subset of my data by
comparing it to another cluster of another subset. Does there exist a
quantitative similarity measure that can be applied?
I am open for any suggestions,
thx for your help,
Julia
--
_
On 10/12/07, Duncan Murdoch <[EMAIL PROTECTED]> wrote:
> On 10/11/2007 6:32 PM, Deepayan Sarkar wrote:
> > On 10/11/07, Karim Rahim <[EMAIL PROTECTED]> wrote:
> >> Thank you for your reply.
> >>
> >> In this graphics context, I'm not sure what the incident or reflected
> >> light rays are.
> >>
> >
Hi all,
I want to constrain observed correlations to be equal. E.g., I want to
find the ML cor(x,y) and cor (w,z) given that cor(x,y) must be equal
to cor(w,z).
I've received some suggestions that were sent directly to me, but
based on these responses, I'm afraid I wasn't clear enough in my
orig
> On 10/12/07, Ben Bolker ufl.edu> wrote:
> >
> > Trying to find a quick/slick/easily interpretable way to
> > collapse a data set.
> >
Another alternative for SQL fans is the sqldf package. I used the MySQL driver
here since SQLite does not support standard deviation.
sqldf("select BROOD,
On 10/12/07, "Julia Kröpfl" <[EMAIL PROTECTED]> wrote:
>
> Is there a package in R that does Q-type factor analysis?
> I know how to do principal component analysis, but haven't found any
> application of Q-type factor analysis.
Q-mode factor analysis is not a separate "type" of factor analysis b
#Hello,
# I have a question about the addition of values in specific columns and
rows of a Data frame.
# Below I have created two data frames, X.df and "Y.df".
## creation of X.df data frame
X<- matrix(0,16,3)
X.df<-data.frame(X)
X.df[,1] <- c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4)
X.df[,2] <- c(1,2,3
My guess is that there's an easier way but this gives what you want.
newY.df<-aggregate(Y.df$Counts, list(Y.df[,1],Y.df[,2]), FUN=sum)
names(newY.df)<-names(X.df)
temp.df<-merge(newY.df, X.df,
by=intersect(names(X.df),names(newY.df)),all=TRUE)
almost.df<-aggregate(temp.df$Counts, list(temp.df[
On Fri, 2007-10-12 at 17:41 +0200, Birgit Lemcke wrote:
> Hello all together,
>
> I (R Beginner) think the answer for my question is easy but anyway I
> have not found the solution.
>
> I am using hist for a frequency histogramm. But the divisions of the
> xaxis is to large
>
> How can I cha
There is a dataset 'm', which has 3 columns: 'index', 'old1' and 'old2';
I want to create 2 new columns: 'new1' and 'new2' on this condition:
if 'index'==i, then 'new1'='old1'+add[i].
'add' is a vector of numbers to be added to old columns, e.g. add=c(10,20,30
...)
Like this:
index old1
runner said the following on 10/12/2007 4:46 PM:
> There is a dataset 'm', which has 3 columns: 'index', 'old1' and 'old2';
>
> I want to create 2 new columns: 'new1' and 'new2' on this condition:
> if 'index'==i, then 'new1'='old1'+add[i].
> 'add' is a vector of numbers to be added to old colu
My data is the following:
Time
Resistance
Temperature
5
2000
4
10
2200
8
15
2500
14
20
2900
20
25
3000
29
30
3100
38
35
3500
46
40
3800
47
45
3900
50
50
4000
51
I would like to create a scatter plot with Time on the x axis, Resistance on
the y axis
basically something along these lines:
plot(Time, Resistance, bty='c')
par(new=TRUE)
plot(Time, Temperature, axes=FALSE, ylab='', xlab='')
axis(4)
On 10/12/07, Keith Cox <[EMAIL PROTECTED]> wrote:
> My data is the following:
>
>
>
>
> Time
>
> Resistance
>
> Temperature
>
>
> 5
>
> 2000
>
> 4
>
>
What problem are you actually having with 'diff'? Now if you are using
'diff', you will get a vector that is shorter by one than the
original. Now do you want to do do something like:
Xbar = Sum{c(Depth[1], diff(Depth))*temp}/Sum(c(Depth[1], diff(Depth))
On 10/12/07, Thomas Miller <[EMAIL PROTE
Due to an intermittent Broadband connection, I cannot check email as
often as I'd like. Please be assured that all your emails are being
received & I will respond to those that require attention as my ISP
Connection allows. CIncinnati Bell cannot advise as to when this will
be fixed, so I ask that
On 10/12/07, Ben Bolker <[EMAIL PROTECTED]> wrote:
>
>Trying to find a quick/slick/easily interpretable way to
> collapse a data set.
>
> Suppose I have a data set that looks like this:
>
> h <- structure(list(INDEX = structure(1:6, .Label = c("1", "2", "3",
> "4", "5", "6"), class = "factor"
> > Here's a solution that takes the first element of each factor
> > and the mean of each numeric variable. I can imagine there
> > are more general/flexible solutions. (One might want to
> > specify more than one summary function, or specify that
> > factors that vary within group should be d
See the EBImage package on Bioconductor.org. It builds on top of
ImageMagick. /Henrik
On 10/12/07, Bio7 <[EMAIL PROTECTED]> wrote:
>
> Dear R users,
>
> in my application i can transfer images to R with the help of Rserve. The
> images come from
> a java application. When i plot a greyscale imag
Hello all
I'm trying to do a simulation that involves identifying the minimum
point between two peaks of a (usually) bimodal distribution. I can do
this easily if there are only two peaks:
CnBdens<-density(Ys/Xs) #probability density function for ratio of Ys
to Xs
for(p in
65 matches
Mail list logo
