similar to: Regression Models

Hello to all, I am having trouble with intregrating a complicated uni-dimensional function of the following form Phi(x-a_1)*Phi(x-a_2)*...*Phi(x-a_{n-1})*phi(x-a_n). Here n is about 5000, Phi is the cumulative distribution function of standard normal, phi is the density function of standard normal, and x ranges over (-infty,infty). My idea is to to use quadrature to handle this integral. But

I would like to be able to grab x and y columns out of a dataframe and then grab all of the columns that are not equal to x or y. I am sure that I am missing something easy. ftbr_UTM_downstream <- (structure(list(site = c("Jennie_Creek_Main_Stem", "Wolf_Pit_Creek_Main_Stem", "Little_Rockfish_Main_Stem_North", "Big_Muddy_Creek_Main_Stem",

I have a model Q=K*A*(R^r)*(S^s) A, R, and S are data I have and K is a curve fitting parameter. I have linearized as log(Q)=log(K)+log(A)+r*log(R)+s*log(S) I have taken the log of the data that I have and this is the model formula without the K part lm(Q~offset(A)+R+S, data=x) What is the formula that I should use? Thanks for all of your help. I can provide a subset of data if necessary.

R 2.12.2 on Scientific Linux 6.4 #works chron(times.="15:00:00", format=c(times="h:m:s")) #doesn't work chron(times.="15:00", format=c(times="h:m")) From chron Manual: The times format can be any permutation of "h", "m", and "s" separated by any one non-special character. The default is "h:m:s". what am I

I have just upgraded to R 2.13 and have library(ggplot2) in my .Rprofile (among other things). when i start R I get an error message. Has something in the start up scripts changed? Is there a better way to specify the library calls in .Rprofile? Thanks for all of the help in advance. Error: Loading required package: grid Loading required package: proto Error in rename(x, .base_to_ggplot) :

s <- 1.00 max(s) returns 1 is there anyway that I can get it to return 1.00. I am using the results of this max statement in a grep statement and it returns the wrong numbers, I will provide more information and code if it would make more sense in context. -- Stephen Sefick ____________________________________ | Auburn University? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? | | Department of

Is there a way to use a continuous variable to pch in qplot? I believe this worked in previous version. I need to specify certain values of a shape for particular points so that multiple graphs all show the same shapes for the same streams. I have gone to the original data and added a pch column that I would like to use to specify the shapes to pch in qplot. Any help would be greatly

All, 1. I will try and make this clear and concise. Please let me know any information that would be helpful in figuring out this problem (I don't know the relevant information to post). I am on linux- see below for session information. 2. Problem: working directory: home an old version of a function is sourced into the R session and doesn't work working directory: Desktop the

Hi folks, I have a question about how to efficiently produce random numbers from Beta and Binomial distributions. For Beta distribution, suppose we have two shape vectors shape1 and shape2. I hope to generate a 10000 x 2 matrix X whose i th rwo is a sample from reta(2,shape1[i]mshape2[i]). Of course this can be done via loops: for(i in 1:10000) { X[i,]=rbeta(2,shape1[i],shape2[i]) } However,

Muy buenas tardes, les escribo con el fin de saber si alguno de ustedes, me puede dar una guia de como generar una simulacion para ilustrar el insesgamiento de los estimadores B0 y B1 de regresion lineal simple lo he logrado hacer para un tamaño de muestra fijo pero no he podido hacer las iterciones. En particular lo que requiero es hacer la simulacion que me grafique para varios diferentes

Saludos a todos. Tengo problemas para encontrar la funcion que te permite generar la formula de la regresion lineal entre dos variables y la linea de tendencia en el gráfico de dispersion de puntos, lo que si puedo realizar es el gráfico con los puntos distribuidos en los 2 ejes pero lo demás que mencioné no se como ingresarlo. Saludos, Gabriel Antonio Trujillo Paucar Estudiante de Biologia

These data represent stream channel cross-sectional surveys. I would like to be able to find the measurement on the tape (measure) where the Bank Full Depth (bkf_depths) is 0. This will happen twice because the channel has two sides. I thought fitting a loess line to these data and then predicting the measurment number would do it. I was wrong. Below is my failed attempt. My naive thought is

When I use this I am getting following warning at the time of opening R from that shortcut : "-max-mem-size=2048MB:too large and taken as 2047M" Why I am getting this? I have 3GB ram installed and using within Vista. Thanks -- View this message in context: http://r.789695.n4.nabble.com/How-to-allocate-more-memories-to-R-tp2271714p2272436.html Sent from the R help mailing list archive

I have a DF ID VALUE 100 120 101 100 102 100 103 104 105 .... when i calculate the sum of the values, it returned NA. should I populate the blank value as 0? Thanks, -- View this message in context: http://r.789695.n4.nabble.com/MISSING-VALUE-IN-R-tp2299586p2299586.html Sent from the R help mailing list archive at Nabble.com.

POS=sum(x[-1][x[-1]>0],na.rm=TRUE) is this the correct syntax? -- View this message in context: http://r.789695.n4.nabble.com/na-rm-TRUE-tp2299596p2299596.html Sent from the R help mailing list archive at Nabble.com.

Hi user's Does anybody work with wavelets on R? Please I need some help.....!!!! Atte Marize Sim?es -- View this message in context: http://r.789695.n4.nabble.com/Package-wavelets-tp2526023p2526023.html Sent from the R help mailing list archive at Nabble.com.

By any chance is there an R package that may contain functionality similar to the frestimate in Matlab/Simulink. Here is the URL for a description of the frestimate functionality: http://www.mathworks.com/help/toolbox/slcontrol/ug/frestimate.html Thank you again for any feedback.

Dear all, I'm just wondering how to submit a package to cran? I followed the instruction, using anonymous as username and my email address as password, but it didn't connect. Any hints? Yan [[alternative HTML version deleted]]

HI I've a unix machine which is of 16 GB. when i run any R process it takes only 2GB of Memory. How to increase Memory limit. It takes a lot of time to run the process for larger datasets ----- Thanks in Advance Arun -- View this message in context: http://r.789695.n4.nabble.com/Help-in-increasing-R-memory-in-unix-tp4616510.html Sent from the R help mailing list archive at Nabble.com.

The dataframe is id salary 100 500 101 600 102 700 103 800 how can i generate a subsets if salary>600? -- View this message in context: http://r.789695.n4.nabble.com/Filtering-in-R-tp2299603p2299603.html Sent from the R help mailing list archive at Nabble.com.