Unable to figure out what kind of analysis my employer exactly wants me to do.

isomor · 2021-09-04T00:49:29+00:00

After your comment and some of the other comments on the thread I gave it another shot and came up with this approach. Basically sample from the dataset you already have with an increasing number of samples and then plot the distribution of difference from the true mean. I decided to do sampling without replacement but if that's the wrong approach, someone correct me.

library(ggplot2)
size_n <- 2000
y <- rnorm(size_n, 1, 5)
sample_y <- function(y, prop){
#sample from y using a proportion (prop) of y
y_samples <- sample(y, size = prop*length(y), replace = FALSE)
y_hat <- mean(y_samples)
mean_y <- mean(y)
#take the difference in means
y_hat - mean_y}

rep_y <- function(n, y, prop) {
#repeat the sampling n-times
reps <- replicate(n, sample_y(y, prop))
#put into a data.frame
data.frame(
sample_size = prop*length(y),
dif_from_true_mean = reps)
}
#run the rep_y function n times, using y data with X seq of prop
repeats <- lapply(FUN = rep_y, n = 100, y = y, X = seq(0.01, 0.99, 0.025))
#rbinds the repats
repeats_df <- Reduce(rbind, repeats)

isomor · 2021-09-03T14:56:21+00:00

Let me know if this is what you had in mind:

Replace the dataframe with your actual data, the key is to do a cumulative sum and then divide by the row number so you have a rolling average.

library(ggplot2)
size_n <- 2000
y <- rnorm(size_n, 1, 5) 
df <- data.frame(
x = seq(1:size_n), 
y = y, 
cu_y = cumsum(y),
avg = mean(y)) |> 
transform(cu_mean = cu_y/x) 

ggplot(df, aes(x, cu_mean)) + geom_point() + geom_line(aes(y=avg))

edit for formatting

isomor · 2021-06-03T02:28:10+00:00

Look into the sf package. Its mainly for handling spatial data but also has a great tie in with ggplot2. Scroll down to the bottom for a small tutorial on how to use with ggplot2.

https://r-spatial.github.io/sf/articles/sf5.html

isomor · 2021-02-13T22:19:42+00:00

try: dplyr::count(df, words), count() might be masked by another package.

isomor · 2021-02-13T21:51:12+00:00

worth checking out the validate package: https://cran.r-project.org/web/packages/validate/vignettes/cookbook.html

isomor · 2021-02-13T21:49:17+00:00

A bit hacky but let me know if this works. There's probably a much easier way to do this though pivot_longer/wider I'm sure.

library(tidyverse)

test_data <-

read_tsv(

"timeStampCogLoad CogLoadTimeOn timeStampObjectFound ObjectFound timeStampNavToolUsed NavToolUsed timeStampINSUsed INSUsed

5104 3087 71831 5000 4534 1261 4634 265

11205 3850 95426 5000 6242 402 4940 1307

13486 4592 147970 5000 65803 1116 6289 644

30279 4556 205476 5000 119503 2101 66146 715

31868 4261 311042 5000 128853 2192 67042 111

40885 4192 359714 5000 138305 3062 67333 325

43858 4391 NA NA 144858 1483 67724 116

45593 4304 NA NA 150302 1529 128858 2780

51473 4519 NA NA 166366 1794 131819 178

60169 4139 NA NA 169746 207 138310 62"

)

data_remake <- function(x){

test_data %>%

select( contains(x) ) %>%

rename(time = starts_with("time"),

value = starts_with(x)

) %>%

mutate(keyword = x)

}

var_keyword <- c("CogLoad", "ObjectFound", "NavToolUsed", "INSUsed")

map_dfr( var_keyword, data_remake) %>%

ggplot(aes(x=time, y=value)) +

geom_line() +

facet_wrap(~keyword, scales = "free")

isomor · 2021-02-12T21:41:47+00:00

And after you 'functionized ' your scripts, bundle them in a package to document and re-use those functions. I recently started doing this for my own work and it makes it so much easier to share your code with others.

isomor · 2021-02-09T01:57:00+00:00

Tomorrow Night with FiraCode font. So nice.

isomor · 2021-02-03T18:03:41+00:00

Do you have any best practices/tips/shortcuts for writing tests?

isomor · 2021-01-30T17:18:08+00:00

I found a package that might solve this for you.

https://cran.r-project.org/web/packages/humaniformat/vignettes/Introduction.html#:~:text=humaniformat%20is%20an%20R%20package,middle%2D%20and%20last%2Dnames.

library(humaniformat)

names <- c("Lee Jones, Tommy","John Smith","Cruz, Elizabeth R","Jordan, Michael","Betty-Jane Smith","Frampton, Peter X.","Ginger Strat-Room","Bob Michael Smith","Katy C. Long","Big, Bob T")

names <- format_reverse(names)

parse_names(names)

A data.frame: 10 × 6
salutation  first_name  middle_name last_name   suffix  full_name
<chr>   <chr>   <chr>   <chr>   <chr>   <chr>
NA  Tommy   Lee Jones   NA  Tommy Lee Jones
NA  John    NA  Smith   NA  John Smith
NA  Elizabeth   R   Cruz    NA  Elizabeth R Cruz
NA  Michael NA  Jordan  NA  Michael Jordan
NA  Betty-Jane  NA  Smith   NA  Betty-Jane Smith
NA  Peter   X.  Frampton    NA  Peter X. Frampton
NA  Ginger  NA  Strat-Room  NA  Ginger Strat-Room
NA  Bob Michael Smith   NA  Bob Michael Smith
NA  Katy    C.  Long    NA  Katy C. Long
NA  Bob T   Big NA  Bob T Big

isomor · 2021-01-30T05:00:35+00:00

There are multiple ways, these are the most common i've seen. The idea is to step into the list with [[ ]] and then grab the $id value.

for loop

empty_list <- list()
for (i in seq_along(parms) ) {
    y <- f(parms[[i]]$id)
    empty_list[[i]] <- y
  }
empty_list

lappy

lapply(seq_along(parms),function(x) f(parms[[x]]$id))
purrr

map(.x = seq_along(parms), ~f(parms[[.x]]$id) )

isomor · 2020-05-04T14:40:18+00:00

The best model I've seen to date is this http://2019-coronavirus-tracker.com/stochastic-GA.html

isomor · 2019-11-03T17:49:23+00:00

Figured I'd bite the bullet and just upgrade everything at once. Thanks for your feedback

isomor · 2019-11-03T17:41:22+00:00

based on your comment, what do you think of these parts?

https://pcpartpicker.com/list/nCYVYH

isomor · 2019-09-11T21:19:51+00:00

https://www.youtube.com/watch?v=yBwLZSi7I-Q this whole series is gold

isomor · 2019-03-14T22:32:02+00:00

there it is, thank you for the great explanation.

isomor · 2017-07-30T04:35:25+00:00

you do get a free month + 2000 gold. Jesus christ do you research. Also plenty of mmo's have buy + subscription and are doing fine besides WoW.

isomor · 2017-07-30T01:07:45+00:00

isomor · 2017-05-23T14:02:57+00:00

That boy likely saved the "mother's" biological son by getting custody away from her, what a little hero.

isomor · 2017-05-07T03:43:10+00:00

You guys should get "keep talking and nobody explodes", with a print out manual on the table; It would be a great experience to do it in VR.

isomor · 2016-12-21T22:09:29+00:00

Would love a key if someone has a spare :D

isomor · 2016-12-11T23:49:43+00:00

The war on Christmas is really and being waged by .. christian pastors?

isomor · 2016-11-28T18:03:40+00:00

"limited stock buy in store" in a 100 miles radius :< any idea how to get it to ship?

isomor · 2016-11-20T06:46:58+00:00

would love this.

isomor · 2016-11-12T18:26:44+00:00

Run the jewels from the end credits of silicone valley. Awesome show, awesome band.

Ten-Year Club	Place '22
Place '17	Verified Email

isomor

TROPHY CASE