similar to: Logistic regression

I realize this may be more of a math question. I have the following optim: optim(c(0.0,1.0),logis.op,x=d1_all$SOA,y=as.numeric(md1[,i])) which uses the following function: logis.op <- function(p,x,y) { ypred <- 1.0 / (1.0 + exp((p[1] - x) / p[2])); res <- sum((y-ypred)^2) return(res) } I would like to add weights to the optim. Do I have to alter the logis.op function by

Hello all, I'm a beginner in R working on a script that will produce a set of models (linear, polynomial and logistic) for each location in a dataset. However, the self-starting logistic model often fails - if this happens I would like to just skip to the next iteration of the loop using tryCatch. I've looked at a few examples and read the help file, but didn't understand tryCatch

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

Hi I?m trying to fit a nonlinear model to a derivative of the logistic function y = a/(1+exp((b-x)/c)) (this is the parametrization for the SSlogis function with nls) The derivative calculated with D function is: > logis<- expression(a/(1+exp((b-x)/c))) > D(logis, "x") a * (exp((b - x)/c) * (1/c))/(1 + exp((b - x)/c))^2 So I enter this expression in the nls function:

Hello I want to know if the summary of the logistic model with survey Pr (> | t |) to test if the coefficient of the model is significant, ie is the p_valor wald test for the model coefficients, for I am interested to know if the three levels of the variable educational level are significant to the model (significance of handling 0.2), I present below the results of my model

Hi, I have a problem when checking the package 'probdistr' (on probability distributions). I got this warning * checking index information ... WARNING Vignettes with missing or empty \VignetteIndexEntry: [1] "probdistr-chi" "probdistr-contextra" "probdistr-discrete" [4] "probdistr-discrextra" "probdistr-exp"

Hi, First at alll, I beg your pardon because maybe I explained bad my questions (because my low level english) I have asterisk 0.5.0, asterisk-oh323-0.5.6, openh323-1.12.2 and pwlib 1.5.2 compiled and installed. I have modules alsa 0.9.8 compiled and installed My PC has an audio card ac97 chipset intel i810 in its motherboard. I want to use asterisk in this way: |PC linux | --ethernet--

I have a problem ordering by descending magnitude a POSIXt object. Can someone help please and let me know how to work around this. My goal is to be able to order my data by DATE and then by descending TIME. I have tried to include as much info as possible below. The problem stems from trying to read in times from a CSV file. I have converted the character time values to a POSIXt object using the

Dear List, I have a data frame of data taken every few seconds. I would like to subset the data to retain only the data taken on the quarter hour, and as close to the quarter hour as possible. So far I have figured out how to subset the data to the quarter hour, but not how to keep only the minimum time for each quarter hour. For example:

Dear R users, How to change lattice panel label/text from the automatically generated label (based on the conditioning) to our own set of label? for example: someStuff <- data.frame(area = rep(c("SOUTH", "NORTH", "EAST", "WEST"), each = 25), group = rep(c("A","B","C","D"), each = 5),

Hi, I have problem with upssched on windows. Upssched is not executed. I have 2 scripts, 1 for notification in upsmon and second for scheduling in upssched. Monitoring is working fine, script write to log. I'm using binary windows installer 2.6.5-6 from NUT. Here are my configs: --- nut.conf MODE=netclient --- upsmon.conf MONITOR ups_1000 at 192.168.3.95 1 <user> <password> slave

I have defined a daily timeseries for the 365 days of 2003 issuing: myts = ts(dati[,2:10],frequency=365,) > myts Time Series: Start = c(1, 1) End = c(1, 365) Frequency = 365 and mytime = as.POSIXct(strptime(as.character(dati[,1]),format="%Y-%m-%d")) contains the dates from "2003-01-01" to "2003-12-31" How can I combine mytime and myts in order to list

List, Consider the following data. gender mygroup id 1 F A 1 2 F B 2 3 F B 2 4 F B 2 5 F C 2 6 F C 2 7 F C 2 8 F D 2 9 F D 2 10 F D 2 11 F D 2 12 F D 2 13 F D 2 14 M A 3 15 M A 3 16 M A 3 17

Hello all, I don' like to have ( ) around my datetime, when send the csv file to somebody. How to remove it? >library(chron) >mydate <- c("2009-05-03","2009-06-07") >mytime <- c("10:30:00","10:20:00") >mydatetime <- chron(dates=mydate,times=mytime,format=c("y-m-d","h:m:s")) >write.csv(mydatetime)  

On Tue, 10 Mar 2020, Radoslav Chovan wrote: > I have problem with upssched on windows. Upssched is not executed. I > have 2 scripts, 1 for notification in upsmon and second for scheduling > in upssched. Monitoring is working fine, script write to log. I'm > using binary windows installer 2.6.5-6 from NUT. > --- upsmon.conf > MONITOR ups_1000 at 192.168.3.95 1 <user>

Due to the following lines in parse.format: else if (nf == 3) { sep <- "" fmt <- substring(format, first = 1:3, last = 1:3) } If a format code has 3 characters, it will not use a separator: > library(chron) > mytime = times('7:15:00') > format(mytime,'h:m') [1] "0715" - Phil Spector

ut 10. 3. 2020 o 16:34 Roger Price <roger at rogerprice.org> napísal(a): > > On Tue, 10 Mar 2020, Radoslav Chovan wrote: > > > I have problem with upssched on windows. Upssched is not executed. I > > have 2 scripts, 1 for notification in upsmon and second for scheduling > > in upssched. Monitoring is working fine, script write to log. I'm > > using binary

Prof. Harrell, My name is Armida Carbajal, I'm a graduate student intern at Sandia National Laboratories (SNL) and am conducting some research for my thesis project at the University of New Mexico in Statistics for SNL. My project entails a logistic regression and I wanted to create a classification table like the one found in SAS using the function CTABLE. I was running out of ideas on

Hello, I have two gamm question (I am using gamm in mgcv). 1. In have, say 5 time series. Monthly data, 20 year. The 5 time series are from 5 stations. The data are in vectors, so I have fitted something along the lines of: tmp<-gamm(Y ~ s(Year,by=station1)+s(Year,by=station2)+ s(Year,by=station3)+s(Year,by=station4)+

Dear R-forum, I am looking for a good resource/help on working with POSIXct values and controlling the pretty x-axis labels and tick marks for a data VS time plots. Specifically, I wish to do programming logic which creates different vertical ablines calculations based on the range of times which i am working on. The default plot results are adequate, but I would love to make explicit calls on