similar to: how to do this sum?

Dear all, I would like to fit a von mises mixture to my data using the EM Algorithm. Unfornately, existing packages (including circstats) doesn't provide such method. So, does anybody have written such code ? Big thanks, Christophe, (R newbie)

Hi all, I am trying to use GAM to work on some data... Are there any resources providing hands-on tutorial/guide on how to do GAM on data in R? Specifically, I am not sure about which model to choose, and smooth models with which effective degree-of-freedom shall I use... I knew there is a book titled: GAM: an introduction using R. Unfornately our local library does not have it... so that's

Hi all, I don't understand the behavior of is.integer(): > x <- integer() > is.integer(x) [1] TRUE > x <- 10 > is.integer(x) [1] FALSE > x <- 1:10 > is.integer(x) [1] TRUE Why is.interger() returns FALSE if x has only one element? And how can someone check if x is an integer but contains only one value? (R 1.5.0 on Linux i386) Thanks, Sven

Hi,R_Help! I have done something like this: >x<-seq(412,612,1) >plot(x,dpois(x,512),col="blue",type="l") >plot(x,dnorm(x,512,sd=sqrt(512)),col="red",type="l") And now,I want to plot the two lines in the same picture or with the same "plot" action? What should I do? And any advises? Thank you in advance! 仭仭仭仭仭仭仭仭仭仭仭仭仭仭仭仭伖併佇伹

Here is an R problem I am struggling with: My dataset is organized like this... subject session variable_x variable_y 01 1 1 <interger values> 01 1 2 01 1 3 01 2 1 01 2 2 01

Dear R-community, my data set looks like 'mat' below. Plant<-c(NA,1,1,1,NA,NA,NA,NA,NA,1); Value1<-rnorm(1:10); Value2<-rnorm(1:10); mat<-cbind(Plant,Value1,Value2); I receive data from two different sites. One site is identified by an interger number, the other site has no data in column Plant=NA. My pb: I'm trying to assign labels "A" or "B" to

Hi All, I am trying to run several linear regressions and print out the summay and the anova reslts on the top of each other for each model. Below is a sample progarm that did not work. is it possible to print the anova below the summary of lm in one file? thanks for your help ###################################################### data<-read.table("data.txt", header=T,

In the R Document, the usage of the acf() is as follow: acf(x, lag.max = NULL, type = c("correlation", "covariance", "partial"), plot = TRUE, na.action = na.fail, demean = TRUE, ...) But now I want to get the result picture like: plot(x,type="l") or plot(x,type="p") How can I do this with acf function? 仭仭仭仭仭仭仭仭仭仭仭仭仭仭仭仭佒伮 伬侎仯仭

Hi All, I'm a little stumped by the following problem. I've got a dataset with the following structure: idxy ix iy country (other variables) 1 1 1 c1 x1 2 1 2 c1 x2 3 1 3 c1 x3 . . . . . 3739 55 67 c7 x3739 3740 55 68 c7 x3740 where ix and

>> In other words, abandoning overflow detection makes the >> duplication of types redundant, while requiring it would be a >> great burden on CPUs that don't have overflow exception hardware. > >Yes, you're right. This has been a desired change for quite some time >now. Unfortunately, its a huge impact to nearly every part of LLVM. We >will

Sorry bit of a Newbie question, and I promise I have searched the forum already, but I'm getting a bit desperate! I have over-dispersed, zero inflated data, with variance greater than the mean, suggesting Zero-Inflated Negative Binomial - which I attempted in R with the pscl package suggested on http://www.ats.ucla.edu/stat/R/dae/zinbreg.htm However my data is non-integer with some pesky

Chris Lattner wrote: > On May 5, 2009, at 8:09 PM, Talin wrote: > > >> I wanted to mention, by the way, that my need/desire for this hasn't >> gone away :) >> >> And my wish list still includes support for something like uintptr_t >> - a >> primitive integer type that is defined to always be the same size as a >> pointer, however large or

Hi list I loaded the content of a file dureetasks.txt to variable a. This file contains an interger per line. when I print a vector, it displays correctly. however, when I try to print the histogram, I get this error message > a=read.table("dureetasks.txt") > hist(a) Error in hist.default(a) : 'x' must be numeric Can you help please? regards -- PhD candidate in

Hello All, Thanks in advance for all help I have the following vector of strings that I want to split.. > y = c("1/1","2/2", "3/3", "4/4") > y [1] "1/1" "2/2" "3/3" "4/4" > z = strsplit(y,"/") > z [[1]] [1] "1" "1" [[2]] [1] "2" "2" [[3]] [1]

Hello, I am a beginner of R. I want to plot the vectors Y vs. X where X has continuous interger values. How can I use the name of each element as the tick label? Thanks a lot! Haiyong

Hello, LLVMers. I have a quick question for constant folding. When I treated a code twice by changing the data type from 'floating point' to 'integer,' What I found was: optimization is done better for the code having int's. E.g., for this simple code: ------------------------------------------------------------- double foo() { double i,j; double sum = 0; for (i=0;

Hola colegas: He estado usando R Studo para obtener unos reportes que espero me ayuden a comprender un ejemplo que estoy tomando de la red. Usando sink("Tabla.txt") , al pedir summay(n), obtengo: sr tl 1 : 1 1: 4 10 : 1 2:13 18 : 1 3:10 30 : 1 4: 4 47 : 1 7: 9 67 : 1 (Other):34 d Min. : 1.000 1st Qu.: 2.000 Median : 5.000

Dear all, I tried to run WinBUGS model in R by using package R2WinBUGS, but I failed because of the massage : cannot calculate DIC for model in log file of WinBUGS14. In fact, the model works in WinBUGS and I can get summaries of the model, such as density, stats..., except DIC , because DIC may not be appropriate in this kind of model due to some reasons(according to the manual of WinBUGS14),

Hi, I would like to implement macro debug info support in LLVM. Below you will find 4 parts: 1. Background on what does it mean to debug macros. 2. A brief explanation on how to represent macro debug info in DWARF 4.0. 3. The suggested design. 4. A full example: Source -> AST -> LLVM IR -> DWARF. Feel free to skip first two parts if you think you know the background.

Hello, I'm posting this to receive some comments/hints about a rather statistical than R-technical question ... . In an anova of a lme factor SSPos11 shows up non-significant, but in the t-test of the summay 2 of the 4 levels (one for constrast) are significant. See below for some truncated output. I realize that the two test are different (F-test/t-test), but I'm looking for for a