similar to: sm.density

Hello to everyone Does someone know where is the bug in sm.logit function in sm library? The estimated probabilities are always 1/2 I'm new to this list, so I don't know if this list is about this kind of things. Sorry for my english, and thanks in advance Oscar -------------- next part -------------- An HTML attachment was scrubbed... URL:

I'm trying to figure out a good way to smooth binary data. The ideal approach appears to be the "sm.logit" function in library "sm", but I haven't had success with it. Below, is some code illustrating what I've come up with so far, but I'm hoping there is a better approach. I'm using R 1.2 (development) under Windows 98. -Bill library(MASS) data(birthwt)

Hi R-users,   I wrote a code to evaluate double sum as follows:   ff2 <- function(bb,eta,z,k) { r <- length(z) for (i in 1:r) { sm1 <- sum((z[i]*bb/2)*(psigamma((0:k)+eta+1,deriv=0)/(factorial(0:k)*gamma((0:k)+eta+1))))  sm2 <- sum((besselI(z[i]*bb,eta)*log(z[i]*bb/2) - sm1)/besselI(z[i]*bb,eta))  sm2 } ff2(bb,eta,z,10)     but it gave me the following message:   >

Hi, I'm attempting to add a "Julian Day" column to a data frame. Here is my code and the resulting data frame: vic.data <- read.table("C:/VIC/data/vic.data.csv", header=F) names(vic.data) <- c("year", "month", "day", "precip", "evap", "runoff", "baseflow", "Tsup",

Sure there is, Install debian, follow my howto and you will have success. Just, your using an .local domain, and thats a reserved name for apples mDNS (zeroconf) And should not be used. ( same for .lan ) https://wiki.samba.org/index.php/FAQ#Can_I_Use_the_.local_Top-level_Domain_for_My_AD_DNS_Zone.3F So the info is good, thats not the problem, finding it, is. Can you post your /etc/hosts

RID solved my problem. But while reading docs I saw new things and I changed my smb.conf completely. I have read almost every parameter but i'm still not %100 sure. Can you do me a last favor? Please can you tell me do I have any problem with new smb.conf? Kernel: Linux 4.14.13-1-ARCH Filesystem: zfs-linux 0.7.5.4.14.13.1-1 Thank you so much for your help. --------------------- [global]

Dear List, I have been trying to use library sm to do density estimation on a 3D dataset. I am using the current MacOS X binary of sm from CRAN. If I do this on a 2D dataset, sm.density returns a list including the component "estimate" which contains the density estimate over a uniform grid. When doing this with 3D data, although I get a nice plot (even when I don't ask for one),

I am running the sm.density.compare function amd I am getting the following error:   my code is  > sm.density.compare(LBSTRESN,COHORT,xlab="Units = umol/L"subset = LBTEST=="Creatinine")   Error in if (from == to) rep.int(from, length.out) else as.vector(c(from,  :   missing value where TRUE/FALSE needed   I do not understand the error and I have had no help when searching

Dear all, Let say I have a sets of numbers: rno1 = rnorm(1000) rno2 = rnorm(1000) If I write request as follow: sm.density.compare (rno, rno3, xfit=min(rno), max(rno2)) why I receive a lot of curves in my plot, while I have only two data sets? regards Przemek -- View this message in context: http://r.789695.n4.nabble.com/sm-density-compare-a-lot-of-curves-tp4253039p4253039.html Sent from the

Hi! I have two dimensional dataset which has and I need to decide if a point lies in some "confidence level". If a point has low confidence/density it can be anomaly which I need to find. For example: #load library library(sm) #get some data x.locs = c(74, 74.5, 75, 77,74.5) y.locs = c(64, 63.5, 63, 61,61.5) points = cbind(x.locs, y.locs) #plot it plot(points) #get points density

I've been using the sm.denstiy.compare function and would like to know some specifics: a) What kind of bootstrap sampling is used, with or without replacement? b) A p-value is returned after running the following script (for example): library(sm) y <- rnorm(100) g <- rep(1:2, rep(50,2)) sm.density.compare(y, g, model="equal") to what statistic does this p-value correspond?

Dear R users, I used sm.density function in the sm package and kde2d() in the MASS package to estimate the bivariate density. Then I calculated the Kullback leibler divergence meassure between a distribution and the each of the estimated densities, but the asnwers are different. Is there any difference between the kde2d and sm.density estimates? if there is a difference, then which is the best

> > (Um. So what device were you using?) In terms of hardware, I'm using a Dell Latitude CPx (not quite sure on the RAM, but the hard disk is 2GB and the server its linked to is much larger...I'm running on the unit's hard drive though). The printer is an EPSON LP-9300. I'm printing directly from the window one gets in R by typing plot(...). Another quick note, its

Hi, I'm using R 1001 for Windows NT and the sm library. I'm trying to create plots for my data set like Bowman and Azzalini have in Figure 1.8 (p. 9) of their book for my data (i.e. a contour plot for each group in my data set and its all plotted on 1 plot). The problem I'm having is that R is not drawing closed contour lines for each group. Sometimes it does; other times it

Is it possible to generate a density plot comparing several (6 or 7) groups of data and have the x-axis in date format (e.g.: "%d%b", 10Mar)? When I try to input the data in date format I get an error saying "this function only allows 1-d data". Thank you, Claudia [[alternative HTML version deleted]]

Dear fellows, The two questions are on sm.density.compare(). I compare kernel density estimates of two arrays of data. I'd like to scale up y-axis so that I can show better the differences in y values. English is not my first language so I'll try to explain it. I would like to stretch y-axis a bit longer but not to change the range of y values. How can I do this? Second, I'm

The trailing backslash seems to be the culprit. After looking at rw1020\bin\INSTALL, I noticed $pkg =~ s/\/$//; which I assume is to strip the trailing slash from a package name before obtaining the package name and cleaning out its directory I added $pkg =~ s/\\$//; immediately after this line. Rerunning Rcmd INSTALL C:\rw1011\src\library\bqtl\ with that change seemed to work OK

Ok, one last try. R1.2.0, WinNT 4.0: [R is installed in: R_HOME=D:\Programme\R\rw1020] Can anyone give a hint what the following error message (when running Rcmd INSTALL) means and what to do about it? cd D:\Programme\R\rw1020\src\library D:\Programme\R\rw1020\bin\rcmd install testfunctions make: Entering directory `/cygdrive/d/Programme/R/rw1020/src/gnuwin32' dlltool -k --as as --dllname

Colleagues ---------------------------------- System info: R version rw1020 on NT ESS using emacs ver. 20.4 ~/.emacs ; Default function template (setq ess-insert-function-templates t) ---------------------------------- When loading a new function into the R workspace using ctrl-CD I would like the default function template to come up, rather than a blank buffer. I thought I had the Default

Dear all, For a population of about 200 I have a continuous variable and an ordinal variable. The question I would like to ask is whether the continuously increases (or decreases) as the rank of the ordinal variable increases. I was thinking that a Spearmen's rank correlation or or a chi squared trend might be appropriate. I don't have any experience dealing with ordinal variables so I