similar to: Simulating distribution of max of two die

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

Hi Is it possible to use the rep command (or maybe some other "shortcut") to create a vector with the following format: 4 3 4 2 3 4 1 2 3 4 Any help is much appreciated. Thanks Jacob Jacob L van Wyk Department of Statistics University of Johannesburg, APK Box 524 Auckland Park 2006 South Africa Office Tel: +27 11 559 3080 Fax: +27 11 559 2832 [[alternative HTML version

Hi I have installed the new version 2.7.1, and get the following messages. Is there some way I can "fix it" ? I use SciViews as well, so the problem may be with one of the libraries that is loaded at start-up: Start-up: Loading required package: datasets Loading required package: utils Loading required package: grDevices Loading required package: graphics Loading required package:

Hallo I would appreciate any help with the following. Problem: Suppose A, B and C are independent U(0,1) random variables. What is the probability that A(x^2) + Bx + C has real roots? I have done the theoretical work and obtained an answer of 1/9 = 0.1111. Now I want to show my students to get this in R with simulation. Below are two attemps, both giving the answer to be about 0.26. Could

I am trying to reproduce Figure 10.5 of Loader's book: Local Regression and Likelihood. The code provided in the book does not seem to work. I have managed (a while ago) to get the accompanied R-code for the figures in the book (file called lffigs.R) from somewhere - cannot find it on the web anymore. The code in the .R script file does not work either. Could anybody please direct me in

Hi I think there was a built-in function (perhaps C-code) called "band_den_bin" in the base library (used for kernel density estimators). Is there a newer version of this function. I have an old routine using this but it won't run in the newer version of R. Any help is much appreciated. Jacob Jacob L van Wyk Dept of Statistics University of Johannesburg, APK Campus Box

Hi I have 3 vectors, say g1, g2 and f, where g1 and g2 contain the row and column indices of a matrix, with corresponding elements in f. How can I create a matrix containing the values in f ? For example: g1 g2 f 3 4 10 2 1 30 4 4 50 How can then get a matrix that looks as follows: 0 0 0 0 30 0 0 0 0 0 0 10 0 0 0 50 I am still learning R. I think reshape will

Thanks to Greg and Michael for their advice. I have tried all these possibilities (and other); nothing seems to work; I persistently get the message: In open.connection(con, "r") : cannot open: HTTP status was '407 Proxy Authentication Required' It seems as if "proxy authentication required" is causing the problem. It seems I have to set up an .Renviron file, but I am

Could somebody please assist with the problem I have now. At our university certain changes have been made to users regarding accessing the internet. I have re-installed the (current) version of R, using custom settings (for Windows XP), choosing amongst other things "Internet2" somewhere. Now, when I wish to install new packages, I first have to select a CRAN mirror. R gives me a list

I have just installed the latest versions of R and Tinn-R (running Windows XP prof.) R 2.11.0 Tinn-R version 2.3.5.2 Everything seems fine, except for the following: I usually do this: Open Tinn-R and click on the R icon to open R - this splits the screen into two parts horiozontally, with Tinn-R on top and R in the bottom window. My problem is simply this: By highlighting one line, and click the

I am trying to find studentized bootstrap intervals for the skewness of a data set. I have tried the following (nerve.dat is a set of data containing observations on one variable) (using Windows XP): x <- scan("e:/Flashbackup2009/Nonparametrics/nerve.dat") n <- length(x) library(e1071) skewness(x) library(boot) sampleskew <- function(x,d) {return(skewness(x[d]))} bb <-

Hi Thanks again to Greg, and Prof Ripley. Greg, I was trying to send you an email, but has been unsuccessful. Speaking to someone from our IT deparment, I am trying to follow your advice to implement the settings to my computer (using Windows), by using Control Panel. I have selected System in Control Panel and chose System, Advanced tab, Environmental Variables, and created the user variables:

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Hi all, The kind of configuration that I use in my sip.conf to connect to various commercial SIP providers looks like this: [general] context=incoming-calls canreinvite=no qualify=yes register => jwinius:passwrd at sip.provider.com/0201234567 [provider] type=peer host=sip.provider.com fromuser=jwinius secret=passwrd This works. However, how would I

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Hi I would appreciate if anyone could help me with an elegant solution of the following: I have a matrix that contain a small number of -Inf values. How can form a new matrix from the old one that excludes all the rows with -Inf values ? Thank you. Jacob Dept of Statistics University of Johannesburg South Africa [[alternative HTML version deleted]]

Dear R users I have a big matrix like 6021 1188 790 290 1174 1015 1990 6613 6288 100714 6021 1 0.658 0.688 0.474 0.262 0.163 0.137 0.32 0.252 0.206 1188 0.658 1 0.917 0.245 0.331 0.122 0.148 0.194 0.168 0.171 790 0.688 0.917 1 0.243 0.31 0.122 0.15 0.19 0.171 0.174 290 0.474

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need