Displaying 20 results from an estimated 1200 matches similar to: ".Internal(optim)"
2003 May 18
1
Fisher LDA and prior=c(...) argument
hello,
I am using LDA and QDA function of MASS library. I understand Fisher LDA
is a method non-probabilistic in nature, so I wonder what happens when I
try to predict my test set examples as in:
> fit <- lda(labels~., data=train.table, prior=c(.5,.5))
> pred <- predict(fit, data=test.table, prior=c(.5,.5))
Specifically I ask this because in my problem there are 700 examples
2003 Jul 07
1
'postscript' command within a function
hello all,
I am trying to print a ps file as part of a function as in:
func <- function (..., filename="temp.ps") {
# some stuff
[...]
# plot
eval( cat("postscript(\"",filename,"\")\n", sep="") )
plot(...)
abline(...)
dev.off()
# more stuff
[...]
}
but it does not work. Nor it does with 'paste' instead
2003 Jun 25
2
dendrograms
Hello all,
I am using libraries (mva,cluster) to produce dendrograms. With 1000
examples the dendrogram gets too crowded, and i am wondering whether there
is an option (which i cannot find) to set the number of leaf nodes, like
in matlab, and return the plot and the assignment map examples -> leaf
nodes. Any suggestion is appreciated. Thanks
Edo
2008 Nov 19
1
How to look within .Internal ?
In the optim() function there is a syntax :
res <- .Internal(optim(par, fn1, gr1, method, con, lower,
Here how can I see the inside-codes of ".Internal" function ?
Regards,
--
View this message in context: http://www.nabble.com/How-to-look-within-.Internal---tp20581897p20581897.html
Sent from the R help mailing list archive at Nabble.com.
2005 Apr 20
2
how to get code of a .Internal() function?
Hello,
I'm working with the function optim() from stats
package,
and inside this function is called the function
.Internal(optim(....)) and I want to get the code of
this function which would help me to understand why
the Nelder-Mead algorithm doesn't converge with my
data.
I'm working under Windows XP.
Could you reply to this adress because I don't belong
to the mailing list
2006 Mar 08
1
power and sample size for a GLM with Poisson response variable
Craig, Thanks for your follow-up note on using the asypow package. My
problem was not only constructing the "constraints" vector but, for my
particular situation (Poisson regression, two groups, sample sizes of
(1081,3180), I get very different results using asypow package compared
to my other (home grown) approaches.
library(asypow)
pois.mean<-c(0.0065,0.0003)
info.pois <-
2003 May 25
1
LDA once again
hi there,
i have one more question about LDA. just to make surei understand,
suppose we have two classes, then if i specify a prior=c(.3,.7) in
lda(...) this will affect my between classes covariance matrix as in:
SB = (.3*m1 - .7*m2) %*% inv(Sigma) %*% t(.3*m1 - .7*m2)
[is Sigma affected ?] and the threshold to decide which class to assign
'test' data = log(.3/.7)
if i specify a
2006 Jul 04
1
problem getting R 2.3.1 svn r38481 to pass make check-all
Hi,
I noticed this problem on my home desktop running FC4 and again on my
laptop running FC5. Both have previously compiled and passed make
check-all on 2.3.1 svn revisions from 10 days ago or so. On both these
machines, make check-all is consistently failing (4 out of 4 attempts on
the FC 4 desktop and 3 out of 3 on the FC 5 laptop) in the
p-r-random-tests tests. This is with both default
2011 Apr 20
1
Error in dimnames(x) for Poisson EWMA model
I am attempting to run a Poisson EWMA model using Patrick Brandt's source code. I get the following error when I run the code:
Error in dimnames(x) <- dn :
length of 'dimnames' [1] not equal to array extent
Dimnames(x) looks like this:
[[1]]
NULL
[[2]]
[1] "mip" "div" "nom" "unity" "mood"
2007 Feb 14
1
Problem with the 'hist' function
Hi, I am using the following R version:
> R version 2.4.1 (2006-12-18)
> Copyright (C) 2006 The R Foundation for Statistical Computing
> ISBN 3-900051-07-0
I believe I found a bug in the 'hist' function, when
'probability=TRUE'. I looked in the archives and I came across
problems with the 'hist' functions (e.g., bug PR# 944, posted in
2001), however, a
2007 Aug 16
1
R 2.5.1.
installed from the binary image (http://cran.r-project.org/bin/
macosx/) on an intel mac 10.4.10 throws three warnings:
-----------------
R version 2.5.1 (2007-06-27)
Copyright (C) 2007 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or
2004 Jun 30
1
funny plotting
Hi, I just wanted to plot a boxplot with a nice curve going through it,
I thought this would be a simple task but for some reason I can't get
the two graphs on the same page accurately. Enclosed is the code showing
the two plots seperately and together. I would have thought it should
work if I could use boxplot() then plot() overlayed but it won't allow
the command add=TRUE (which
2011 Nov 12
1
Please Help
HiI want to construct a logliikelood function in RHere is the situationy=number of particles emitted in 1 hr period~pois(30)p=probability of detection of radiation particlesx=number of particles detected by a radiation detector~pois(30p)where p~beta(a,1)I have to calculate the loglikehood for a for the range a(2,50)I wish to simulate 100 random samples for each aHere is my code:-m=481n=100x =
2005 Feb 11
3
DVD Shrink and wine 20050111 vs 20041019
hallo
i'm running mepis (debian unstable) and if i install wine from the
regular (non-wine specific) deb repositories, I get wine 20041019.
if I install from wine.sourceforge, I get 200501011.
with the former, the application seems to lock up much more often.
with the latter, the ability to play the movie (e.g. to see which audio
stream to include) is disabled.
i see from this screenshot
2002 Sep 17
3
[ smbfs ] Does it support ACLS?
Hello there friends & foes,
I've had if with M$ windows! Next week I'll sabotage my shop (how 'bout
password protecting dualboots on WS) and force LINUX down my users' throats
;-) (oh, don't I like this tone... makes me feel Billy G.) Oops, showstopper!
ACL enabled network filesystem, where is it? I can't go with alpha code nfs4,
so is there a patch to the smb
2011 Sep 06
2
Generalizing call to function
Hello guys,
I would like to ask for help to understand what is going on in
"func2". My plan is to generalize "func1", so that are expected same
results in "func2" as in "func1". Executing "func1" returns...
0.25 with absolute error < 8.4e-05
But for "func2" I get...
Error in dpois(1, 0.1, 23.3065168689948, 0.000429064542600244,
2010 Jun 21
1
glm, poisson and negative binomial distribution and confidence interval
Dear list,
I am using glm's to predict count data for a fish species inside and outside
a marine reserve for three different methods of monitoring.
I run glms and figured out the best model using step function for each
methods used.
I predicted two values for my fish counts inside and outside the reserve
using means of each of the covariates (using predict() )
therefore I have only one value
2003 Jun 09
2
A stupid question...
su-2.05a# ping 192.168.1.1
PING 192.168.1.1 (192.168.1.1): 56 data bytes
ping: sendto: No route to host
ping: sendto: No route to host
^C
--- 192.168.1.1 ping statistics ---
2 packets transmitted, 0 packets received, 100% packet
rl1: flags=8843<UP,BROADCAST,RUNNING,SIMPLEX,MULTICAST> mtu 1500
inet 192.168.1.1 netmask 0xffffff00 broadcast 192.168.1.255
inet6
2011 Mar 18
1
[Patch suggestion] Adding 3rd arg to tempfile() to set extension
The other day I was working on an example which used tempfile() to create
file for use by the graphics device. And while I love tempfile()---as it is
portable and clever and the files get cleaned by R and all that---I noticed
one missing feature I would like to see: beside a starting name pattern, and
an optional directory, an 'file extension' argument would be nice to have.
As e.g. in
2005 Mar 17
1
adding winehq debian repository to apt-get | synaptic
I added the winehq debian repository to apt-get | synaptic
(I am running mepis 3.3) and would like to install those rather than the
debian official unstable builds
but
since the winehq repository is also marked as part of 'unstable' distro
I can't find a way to distinguish the wine / winelib provided by winehq
from the official ones (which would need the many winelib-* packages)