Here are some ideas you might like to consider
par(mar = c(5,4,2,4)+0.1, yaxs = "r")
Sample <- rgamma(1000,2.5,.8)
hist(Sample, main = "", freq = FALSE, ylim = c(0,1))
pu <- par("usr")[1:2]
x <- seq(pu[1], pu[2], len = 5000)
y <- pgamma(x, 2.5, 0.8)
par(new = TRUE)
plot(x, y, type = "l", axes = F
Dear R Helpers,
I have trouble applying reShape and reshape although I read the documentation
and several posts, so I would very much appreciate your help on the two points
below.
I have a dataframe
df = data.frame(Name=c("a", "a", "a", "b", "b", "c"), X1=c("12", "13", "14",
"20", "25", "30"
Hi,
Is it possible to overlay a cummulative distribution function on a
histogram of a gamma distribuition.
I have a gamma function
Sample = rgamma(1000,2.5,.8)+1.5
hist(Sample)
regards
--
View this message in context:
http://www.nabble.com/Overlay-cdf-tp23515551p23515551.html
Sent from th
On Tue, 12 May 2009, Mark Na wrote:
H R-helpers,
I am using the following code to make a dataframe from an Excel spreadsheet:
library(RODBC)
channel <- odbcConnectExcel("Spreadsheet.xls")
Data <- sqlFetch(channel, "Tab1")
odbcClose(channel)
One column (several, actually) in the spreadsheet co
I don't understand. The help page for "maxLik" says that if the
"grad" argument is 'NULL', "numeric gradient will be used."
If this does not answer your question, please read the posting
guide "www.R-project.org/posting-guide.html" and provide commented,
minimal, self-contained, r
Hi all,
I have posted similar questions regarding this topic, but I just can't seem
to get over the hump and find a straightforward way to do this. I have
attached my file as a reference.
Basically, the attached file is a 5 degree by 5 degree grid of the the world
(2592 cells), most of them are N
Hi,
I would like to create one postscript file with two different xyplots (which
use different datasets).
Is there a (relitively easy) way to do this?
Many thanks!!!
Liat.
--
View this message in context:
http://www.nabble.com/creating-a-postscript-file-with-two-xyplots-tp23515202p23515202.html
joaks1 wrote:
I need to perform some calculations with some extremely small numbers
The R package Brobdingnag does this. It uses logarithmic
representation, as recommended by others in this thread, but wraps it
all up for you in a custom numeric class so you can use them as any
other numeri
Anyone has some clue to this question?
Thanks in advance!
lehe wrote:
>
> Hi,
> I was wondering how to specify the line type for line instead of for bar.
> Here is my code:
> "plotCI(x=mcra1avg, uiw=stdev1, type="l",col=2,lty=2)"
> This way, I will have the bar line as dashed "lty=2" and red
Farrel Buchinsky wrote:
Is R an appropriate tool for data manipulation and data reshaping and data
organizing? I think so but someone who recently joined our group thinks not.
The new recruit believes that python or another language is a far better
tool for developing data manipulation scripts th
or use nls.lm as in
install.packages("minpack.lm")
library(minpack.lm)
x <- c(2, 8, 14, 20, 26, 32, 38, 44, 50, 56, 62, 68, 74)
y <- c(100, 99, 99, 98, 97, 94, 82, 66, 48, 38, 22, 10, 1)
res <- function(p, x, y) y - ff(p,x)
ff <- function(p, x) 100*exp(p[1]*(1-exp(p[2]*x))/p[2])
aa <- nls.lm(par
On 13/05/2009, at 2:22 PM, Debbie Zhang wrote:
Dear R users,
Does anyone know how to graph the function below?
sqrt(2)Γ(n/2)/[sqrt(n − 1)Γ((n − 1)/2]
Yes, pretty well everyone knows how to graph it.
Please help.
No. Do your own homework.
cheers,
?curve
just create an R expression for the equation and then plot it. I am not
sure exactly what your expression is supposed to be.
On Tue, May 12, 2009 at 10:22 PM, Debbie Zhang wrote:
>
>
> Dear R users,
>
> Does anyone know how to graph the function below?
>
> sqrt(2)Ã(n/2)/[sqrt(n - 1)Ã((n
Dear R users,
Does anyone know how to graph the function below?
sqrt(2)Î(n/2)/[sqrt(n â 1)Î((n â 1)/2]
Please help.
Debbie
_
Want to stay on top of your life online? Find out how with Windows Live!
http://windowslive.ni
I posed this question imprecisely earlier and did not specify that the
pairing needs to be one-to-one. A precise version of the question
follows:
Given two numeric vectors of possibly unequal length, I'd like to pair
each element of the shorter vector with an element of the longer
vector
Dear Dmitry,
Take a look at ?nls and its examples.
HTH,
Jorge
On Tue, May 12, 2009 at 5:44 PM, Dmitry Gospodaryov
wrote:
> I have the data:
> for x: 2, 8, 14, 20, 26, 32, 38, 44, 50, 56, 62, 68, 74,
> for y: 100, 99, 99, 98, 97, 94, 82, 66, 48, 38, 22, 10, 1.
> y depends on x by equation: y =
hello,
I have an code for similar plot in 2D.
HTH
Cleber Borges
--
#
#
#
t
Hi
MikSmith wrote:
> Hi
>
> I've been playing with a 3x2 graphics device using the default size as it
> appears on screen. This has given me tall thin plots which I can resize by
> dragging the window and increasing the window width. However I was wondering
> if I can force R to produce square p
Greetings -
Am hoping to use the strucchange package to look for structural breaks
in some messy regression data. A series of preliminary analyses indicate
that BLUE for these data will involve some weighting the data (estimates
of a particular population parameter) by a function of the varian
I have the data:
for x: 2, 8, 14, 20, 26, 32, 38, 44, 50, 56, 62, 68, 74,
for y: 100, 99, 99, 98, 97, 94, 82, 66, 48, 38, 22, 10, 1.
y depends on x by equation: y = 100*exp(b*(1-exp(c*x))/c),
where b and c are coefficients. I need to find coefficients in this
equation for given data. How can I do
So, given these vectors, what's the right answer?
a <-c(1,2,3)
b<-c(2, 3, 4 ,5, 6)
2,2,3 (distance =1) or 2,3,4 (distance=3)?
I.e., is recycling of elements allowed or must the pairing be one-one?
if the first, then:
df <- expand.grid(a,b)
df$dist <- with(df, (Var1-Var2)^2 )
head(df[o
Dear Jacy,
If AT and TA always one after the other, you might consider the following as
an alternative:
res <- apply(a, 2, function(x) c(x[1], sum(x[2:3]), x[4] ))
rownames(res) <- rownames(a)[-3]
res
#[,1] [,2] [,3] [,4]
#AA159 13
#AT5 13 21 29
#TT48 12 16
HTH
In the first reply, what was calculated was the overall means by group (amino
acids). It does not work for a larger database.
I am quite really new to R, and I worked on your question just to learn how
to manipulate data with R.
The following seems to work. The code could be made a lot more elegan
Hi, Stephen:
Have you discussed this with any of your professors? With a
little luck, you might find the right prof to work with who could help
you select an ecology journal and write an article for that journal
giving an overview of your package. I suggest you think in terms of a
2-
Anyone has some clue to this question?
Thanks in advance!
lehe wrote:
>
> Hi,
> I was wondering how to specify the line type for line instead of for bar.
> Here is my code:
> "plotCI(x=mcra1avg, uiw=stdev1, type="l",col=2,lty=2)"
> This way, I will have the bar line as dashed "lty=2" and red "c
Hi,
I got it compiled.
> It says "recompile with -fPIC" But I am already
> doing so, am I not?
It appeared so, but in fact that wasn't what happened.
When I looked at the date of CConverters.o, I realized
the date was several days old. I deleted the entire
tree and did everything once more (fr
Can you provide commented, minimal, self-contained, reproducible code.
You can check out 'duplicated' to remove duplicates.
On Tue, May 12, 2009 at 7:06 AM, Xin Shi wrote:
> Dear:
>
>
>
> I am trying to merge two tables by a common variable. However, there are a
> few same variables which are in
Afternoon everyone,
I've spent the last week or so looking at all the experimental design
packages I can find in R. AlgDesign, design.conf and BHH2 being the best one
I could find.
Unfortunately none of these do a particularly good job for complex designs,
in particular for conjoint or disc
H R-helpers,
I am using the following code to make a dataframe from an Excel spreadsheet:
library(RODBC)
channel <- odbcConnectExcel("Spreadsheet.xls")
Data <- sqlFetch(channel, "Tab1")
odbcClose(channel)
One column (several, actually) in the spreadsheet contains integers in
its first few rows b
Have you worked through "vignette('dlm')"? Vignettes are nice
because they provide an Adobe Acrobat Portable Document Format (pdf)
file with a companion R script file, which you can get as follows:
(dlm. <- vignette('dlm'))
Stangle(dlm.$file)
The first of these two lines opens
Hi,
i am trying to find an appopriate R function which will estimate the bias
associated with the apparent error rate of my logistic discriminant model
(groups = 2, covariates = 3). I have read that bootstrapping can be used for
this. Does anyone have any ideas on how I can go about doing this?
I
rgl will rotate and zoom, but I don't know if you can select points.
On Tue, May 12, 2009 at 4:35 PM, Abby Drake
wrote:
> It looks like Rcmdr may be able to select points on 3D scatterplots
> however when I try to use it's 3dscatter plot function I get the error
> message: could not find functio
It looks like Rcmdr may be able to select points on 3D scatterplots
however when I try to use it's 3dscatter plot function I get the error
message: could not find function "nice"
If I copy the code:
scatter3d(data$X, data$Z, data$Y, surface=FALSE, residuals=TRUE, bg="white",
+ axis.scales=TRUE,
Jonas Stein wrote:
Thank you, Baptiste and Charlie.
I found some examples wich look great on:
http://www.texample.net/tikz/examples/
Perhaps try the pgfSweave package on r-forge?
http://r-forge.r-project.org/projects/pgfsweave/>
[..]
accomplished using the java utility eps2pgf and require
Or with scatterplot3d:
The regular identify() usage with the iris data:
s3d <- scatterplot3d(iris[,1:3], color=as.numeric(iris[,5]))
identify(s3d$xyz.convert(iris[,1:3]))
or for locator: the values you get are the 2D representation of the
current coordinate system and you can use locator dir
I am trying to apply this technique (M.J Anderson 2001) to a dataset
of aquatic insect abundances. There is a sample in the unrestored and
restored segement of a stream for every time period. I would like to
compare the centroids of the distance matrices for the treatments up
(unrestored) and dn
Hi,
I'm making a database for my company, and I would like to automate some tasks.
One of the tasks that is demanding a lot of time is making up the task schedule.
I'm wondering if it's possible to automate this, and how this can be done.
As the database is made up in postgresql, I can access my d
Hi,
I cannot compile R with shared library.
I am using Redhad Linux on Dell hardware.
Here''s what I am doing:
I set ?PICFLAGS in config.site file:
CPICFLAGS=-fPIC
FPICFLAGS=-fPIC
and issue:
./configure --enable-R-shlib --prefix=$HOME/newR
make
configure finishes with no complaints, the las
Sorry for these multiple postings.
I solved the problem using na.omit() to drop records with missing values
for the time being. I will worry about imputation, etc. later.
I calculated the sum of squared errors for 3 models, linear regression,
neural networks, and support vector machines. This
On 5/12/2009 3:08 PM, Abby Drake wrote:
Hello Everyone,
I am new to R and need some help.
I have a matrix of x,y,z coordinates that I would like to
interactively plot in 3D and then using the cursor select points on
the plot and have the coordinates sent to a matrix. I am using the rgl
package
Thank you, Baptiste and Charlie.
I found some examples wich look great on:
http://www.texample.net/tikz/examples/
>> Perhaps try the pgfSweave package on r-forge?
>> http://r-forge.r-project.org/projects/pgfsweave/>
> [..]
> accomplished using the java utility eps2pgf and requires the PGF/TiKZ
Hello Everyone,
I am new to R and need some help.
I have a matrix of x,y,z coordinates that I would like to
interactively plot in 3D and then using the cursor select points on
the plot and have the coordinates sent to a matrix. I am using the rgl
package to plot the data at the moment because it
Jun,
sapply() does the trick! Thank you Jun! I really appreciate your help.
-Matt
On Tue, May 12, 2009 at 1:20 PM, Jun Shen wrote:
> Forgot one thing, make sure your data is a list or data frame.
>
>
> On Tue, May 12, 2009 at 12:19 PM, Jun Shen wrote:
>
>> I don't think get(factor[i]) will
On Tue, May 12, 2009 at 2:35 PM, Zeljko Vrba wrote:
> Aha, so you DON'T have only red and green rectangles in the picture, there is
> also white in between, and some pale delimiter lines. Nevertheless, both
> things ease the job slightly, and the approach I described should work. You
> can even
If the matching need not be one-to-one, then you can just compute the
Euclidean distances between the two vectors, then in each row (or
column, which ever corresponds to the shorter vector) and find the
smallest. This should be fairly easy to do.
Andy
From: Thomas S. Dye
>
> Given two numeric
Using traditional ANOVA, you'd have to drop either cases or time
points with missing data. Using linear mixed effects analysis, you'd
be able to use all the data. LME also has the benefit of *not*
assuming sphericity, which is good for data like yours (many measures
across few cases) where the trad
jim holtman schrieb:
Try this:
key <- rownames(a)
key[key == "AT"] <- "TA"
do.call(rbind, by(a, key, colSums))
something like
paste(sort(strsplit(key, split="")[[1]]), "")
might be more general.
__
R-help@r-project.org mailing list
https:/
As a follow-up to my email below:
The input data frame to nnet() has dimensions:
> dim(coreaff.trn.nn)
[1] 50888
And the predictions from the neural network (35 records are dropped -
see email below for more details) has dimensions:
> pred <- predict(coreaff.nn1)
> dim(pred)
Forgot one thing, make sure your data is a list or data frame.
On Tue, May 12, 2009 at 12:19 PM, Jun Shen wrote:
> I don't think get(factor[i]) will work. The problem is get only sees a
> character string "data$focus" instead of doing "extracting focus from data".
> In your case isn't lapply (or
I don't think get(factor[i]) will work. The problem is get only sees a
character string "data$focus" instead of doing "extracting focus from data".
In your case isn't lapply (or sapply) good enough?
sapply (data, summary)
try ?lapply for details
Jun
On Tue, May 12, 2009 at 10:55 AM, rapton wrot
Hello
I'm plotting 6 graphs using "mfrow = c(2, 3)". In these plots, only
graphs in the first column have titles for the y axis, and only the ones
in the last row have titles for the x axis.
I'd like all plots to be of the same size, and I'm trying to keep them
as near each other as possible, but
Given two numeric vectors of possibly unequal length, I'd like to pair
each element of the shorter vector with an element of the longer
vector such that the sum of squared differences between the pairs is
minimized. Can someone point me to an R function or an algorithm for
accomplishing th
Hi Mike,
MikSmith wrote:
Hi
I have a dataset of orientation and elevation in a dataframe which I am
plotting in circular. Orientation is on the range 0-360, however I need to
reduce this to 0-180 as the values over 180 are not necessarily correct. I
can do
striae$orientation[striae$orientation
Hello,
I have a data set with many variables, and often I want to run a given
function, like summary() or cor() or lmer() etc. on many combinations of one
or more than one of these variables. For every combination of variables I
want to analyze I have been writing out the code by hand, but given
Hello,
I am using triangle.class in ade4 package and a factor. I would like to use
this graphic parameter (pch) to set a different kind of point for each group
in my factor (lakes). I have tried pch=1:4, I guess it does what it is
supposed to do but not what I want. I can't use colors.
Thank you
Hello Dieter
I got them from a Matrix on with I use the applyfunction tu compute the mean
columns by columns
print(apply(matSD,2,mean))
So it's a kind of vector.
I'll check with the cat function
Thank you, have a nice day
Thomas
--
View this message in context:
http://www.nabble.com/How-to
Dear:
I am trying to merge two tables by a common variable. However, there are a
few same variables which are in both of two tables. How can I take them away
when I merge the two tables?
Thanks!
Xin
__
R-help@r-project.org mailing li
Greetings,
I am using rpart for classification with "class" method. The test data is
the Indian diabetes data from package mlbench.
I fitted a classification tree firstly using the original data, and then
exchanged the order of Body mass and Plasma glucose which are the
strongest/important varia
Hi
I've been playing with a 3x2 graphics device using the default size as it
appears on screen. This has given me tall thin plots which I can resize by
dragging the window and increasing the window width. However I was wondering
if I can force R to produce square plots or set the actual aspect ra
Thanks for showing me how to use genD correctly and that it matches deriv
here. Like you say, there must be a problem with the manual way of doing
it, and I will look at it closely. Thanks again.
spencerg wrote:
>
> Hi, Paul:
>
>
> Your example is so complicated that I don't want to
Totally awsom!
Thank you very much
Thom
--
View this message in context:
http://www.nabble.com/How-to-display-data-content-on-a-only-row--tp23497437p23502952.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mai
Hi
I have a dataset of orientation and elevation in a dataframe which I am
plotting in circular. Orientation is on the range 0-360, however I need to
reduce this to 0-180 as the values over 180 are not necessarily correct. I
can do
striae$orientation[striae$orientation>180]-180)
to extract thes
Hello,
I'm trying to do a comparsion on a large scale say 10L bottle of liquid and a
small scale bottle of liquid 0.5L, I have 5 different samples from each and
they are measured over the space of 8 days as % viability and the % viability
decreases over time. However not all 10 samples got me
>
>
> I am a PhD Student at the University of
> >Agriculture,Abeokuta Nigeria.Also a Statistician at
> National
> > Horticultural Research institites ibadan Nigeria.
> >
> > I currently working on Discriniant analysis and i need
> a
> > help on hoe to solve Lanche Bruch's "hold out
> Procedure"
I am analyzing the voice telephone traffic of some operators. In this type of
phenomenon is useful to consider the anomalies in the duration in minutes of
conversation but it is very important to take account of another variable:
average time of conversation, the relationship between the duration
Uwe,
John's comment about the difficulties with finding polynomial roots is even
more forceful for a system of polynomials. There are likely numerous roots,
some possibly real, and some possibly multiple. Homotopy methods are
currrently the state-of-art for finding "all" the roots, but beware th
Dear R-community,
Dear Prof. Therneau,
I would like to fit an AFT-model with time-dependent covariates and
right-censored data.
Searching the mailing list for information on the subject, I found some old
posts which said it didn't work back then.
My questions:
(1) Has this kind of fitting al
Hi,
I'm having trouble understanding how to construct a random number generator
for a parametric bootstrap. My aim is to bootstrap a Likelihood Ratio
statistic (under the null) for a linear model. The function at this point
is given by
boot.test.n01 <- function(data, indeces, maxit=20) {
y1 <- f
Greetings:
There is a LinkedIn R Group. There are job postings, news, as well
as other discussions.
http://www.linkedin.com/groups?about=&gid=77616&trk=anet_ug_grppro
Regards,
Ajit
http://www.linkedin.com/groups?about=&gid=77616&trk=anet_ug_grppro
__
On 5/12/2009 8:43 AM, Tim Carnus wrote:
Dear List,
I am trying to plot a similar graph to attached from minitab manual in R.
I have a response Y and three components which systematically vary in their
proportions. I have found in R methods/packages to plot ternary plots (eg.
plotri
Greetings:
There is a LinkedIn R Group. There are job postings, news, as well as
other discussions.
http://www.linkedin.com/groups?about=&gid=77616&trk=anet_ug_grppro
Regards,
Ajit
--
_
Ajit Gemunu de Silva
Oakland CA 94619
skype: ajit_de_silva
[[alt
I am exploring neural networks (adding non-linearities) to see if I can
get more predictive power than a linear regression model I built. I am
using the function nnet and following the example of Venables and
Ripley, in Modern Applied Statistics with S, on pages 246 to 249. I have
standardized vari
Hello,
I was wondering if there is a function in R that imports tables directly
from a HTML document. I know there are functions (say, getURL() from {RCurl}
) that download the entire page source, but here I refer to something like
google document's function importHTML() (if you don't know this fun
Hi,
Â
I'm actually Iâm performing a TSLS linear multiple regression on annually
data which go from 1971 to 1997. After performing the TSLS regression, I tried
to extract the R squared value using âoutput$r.squaredâ function and to
perform autocorrelation (Durbin Watson and Breus
John C Nash wrote:
Finding polynomial roots is not a problem where one wants a quick and
dirty code. There are a lot of pitfalls, especially if there are roots
that are multiples, and there has been a lot of work on this problem.
See http://en.wikipedia.org/wiki/Category:Root-finding_algorith
On 12-May-09 12:00:50, Steve Candy wrote:
> Hi All
> This is a general stats problem that I am dealing with using R,
> so any help is greater appreciated.
>
> I have two lognormal distributions with means M1 and M2.
> If we have:
>
> H0: log(M1/M2)=0
>
> H1: log(M1/M2) !=0 equivalent to log(M1/M
Amor,
You have the possibility to use the duplicated() function.
A[!duplicated(A[,var]),]
Atenciosamente,
Leandro Lins Marino
Centro de Avaliação
Fundação CESGRANRIO
Rua Santa Alexandrina, 1011 - 2º andar
Rio de Janeiro, RJ - CEP: 20261-903
R (21) 2103-9600 R.:236
( (21) 8777-7907
( lean...@ces
John,
Nothing about the data set you posted should give issues in any version of
the
package. The "points <=0" usually arises when the survival curve drops to
zero,
and so log(-log(0)) is being trimmed off the graph by the plot command. This
is
what should happen. As Thomas pointed out t
Finding polynomial roots is not a problem where one wants a quick and dirty
code. There are a lot of pitfalls, especially if there are roots that are
multiples, and there has been a lot of work on this problem. See
http://en.wikipedia.org/wiki/Category:Root-finding_algorithms .
And Uwe may not
On 5/12/2009 8:29 AM, amor Gandhi wrote:
> Hi,
>
> I have the following data and I would like to delete douple names, it is
> almost similar to SAS PROC SORT nodupkey! Is there any function in R does
> this?
>
> x1 <- rnorm(11,5,1)
> x2 <- runif(11,0,1)
> nam <-paste("A", c(1:4,4,5:9,9), sep=
Hi,
I was wondering how to specify the line type for line instead of for bar.
Here is my code:
"plotCI(x=mcra1avg, uiw=stdev1, type="l",col=2,lty=2)"
This way, I will have the bar line as dashed "lty=2" and red "col=2", and
the line connecting the centers of the bars is also red "col=2" but solid
*I´m writing to ask you how can I do Survivals Curves using Time-dependent
*covariates? Which packages I need to Install?*
This is a very difficult problem statistically. That is, there are not many
good ideas for what SHOULD be done. Hence, there are no packages. Almost
everything you find
Dear List,
I am trying to plot a similar graph to attached from minitab manual in R.
I have a response Y and three components which systematically vary in their
proportions. I have found in R methods/packages to plot ternary plots (eg.
plotrix) but nothing which can extend it to res
Thanks, Baptiste and Zeljko. I am attaching here an example of the
picture of the rectangles.
Paul
On Tue, May 12, 2009 at 1:23 PM, baptiste auguie wrote:
>
> Depending on the nature of your pdf file, it may be possible to use
> the grImport package. I've never used it before, but a quick test
Dear List,
I am trying to plot a similar graph to attached from minitab manual in R.
I have a response Y and three components which systematically vary in their
proportions. I have found in R methods/packages to plot ternary plots (eg.
plotrix) but nothing which can extend it to response surface
Thom_249 wrote:
>
> I got them from a Matrix on with I use the applyfunction tu compute the
> mean columns by columns
> print(apply(matSD,2,mean))
>
matSD = matrix(round(rnorm(20),2),nrow=4)
cat(matSD)
print(matSD)
dput(matSD) # How to send this matrix to r-help
newMat = apply(matSD,2,mean)
Hi,
I have the following data and I would like to delete douple names, it is almost
similar to SAS PROC SORT nodupkey! Is there any function in R does this?
x1 <- rnorm(11,5,1)
x2 <- runif(11,0,1)
nam <-paste("A", c(1:4,4,5:9,9), sep=".")
mydata <- data.frame(x1,x2)
crownames(mydata) <- nam
Depending on the nature of your pdf file, it may be possible to use
the grImport package. I've never used it before, but a quick test
seems promising,
> # create a test picture
> colorStrip <-
> function (colors, draw = T)
> {
> x <- seq(0, 1 - 1/ncol(colors), length = ncol(colors))
>
To color the error bars in ROCR the same way as the performance curve,
you need to add one more argument (plotCI.col='red') to your plot
call:
plot( perf2,avg="threshold",lty=2,col=2, spread.estimate="stddev", plotCI.col=2)
The use of 'plotCI.col' is an example for the general mechanism of
ROCR t
On May 12, 2009, at 3:50 AM, utkarshsinghal wrote:
Hi All,
Can anybody explain why the following three ways of extracting
residuals from a glm object are giving me different outputs:
> idv = runif(1000,0,1)
> dv = rbinom(1000,1,0.5)
> d = data.frame(idv,dv)
> fit = glm(dv~idv, data=d, famil
Hi All
This is a general stats problem that I am dealing with using R, so any help is
greater appreciated.
I have two lognormal distributions with means M1 and M2.
If we have:
H0: log(M1/M2)=0
H1: log(M1/M2) !=0 equivalent to log(M1/M2)=log(1+P) where P<0 or P>0.
If we calculate the power fo
Thanks Tobias!
A new question: if I want to draw an average ROC from cross-validation, how to
make the bar color same as the line color? Here is my code:
"plot( perf2,avg="threshold",lty=2,col=2, spread.estimate="stddev",barcol=2)"
Even I specify "barcol=2", the color of bars are still black, th
On Tue, May 12, 2009 at 12:20:56PM +0100, Paul Smith wrote:
>
> I have got several pdf files with rows of colored rectangles: red
> rectangles should be read as 0; green rectangles as 1. No other color
> exists. Is there some way to have R reading the colored rectangles to
> a matrix or data frame
> I would like to have no ticks on a scale that represents a factor. The
> tick.number argument from scales does not work in such a situation, as
> the help page as well as this simple (fairly stupid) code show:
>
> require(lattice)
> fact<-gl(4,1,labels=LETTERS[1:4])
> y<-c(1,4,3,2)
> xyplot(y~fa
You could also use rSymPy to symbolically differentiate it. For example,
based on the semi-automatic differentiation example on the rSymPy home
page:
http://code.google.com/p/rsympy/#Semi-Automatic_Differentiation
just replace the indented lines there with the indented lines here
(I've also added
Dear All,
Perhaps, what I am asking is impossible, but I am asking it anyway.
I have got several pdf files with rows of colored rectangles: red
rectangles should be read as 0; green rectangles as 1. No other color
exists. Is there some way to have R reading the colored rectangles to
a matrix or d
dkumar = function(x,a,b) a*b*(x^(a-1))*((1-x^a)^(b-1))
pkumar = function(x,a,b) 1-(1-x^a)^b
(I've based this entirely on the wikipedia entry on the Kumaraswamy
distribution [http://en.wikipedia.org/wiki/Kumaraswamy_distribution],
so best to check both my replication of the formula there and the
ac
Dear Debbie,
On Tue, May 12, 2009 at 5:37 AM, Debbie Zhang wrote:
>
> Dear R users,
>
> Does anyone know how to write function for Kumaraswamy distribution in R?
> Since I cannot write dkumar, pkumar, etc. in R.
>
>
>
> Please help.
>
>
>
> Thanks a lot,
>
> Debbie
>
Check the CRAN Task View on
Hi,
I would like to have no ticks on a scale that represents a factor. The
tick.number argument from scales does not work in such a situation, as
the help page as well as this simple (fairly stupid) code show:
require(lattice)
fact<-gl(4,1,labels=LETTERS[1:4])
y<-c(1,4,3,2)
xyplot(y~fact,scales=l
On 12/05/2009 6:18 AM, Titus von der Malsburg wrote:
Hi list,
I have a function that detects saccadic eye movements in a time series
of eye positions sampled at a rate of 250Hz. This function needs
three vectors: x-coordinate, y-coordinate, trial-id. This information
is usually contained in a
1 - 100 of 136 matches
Mail list logo
