The BCEWithLogitsLoss combines the Sigmoid Activation with BinaryCrossEntropyLoss. If you want to test the outputs of your model later, you have to apply the sigmoid activation to get a logit in between you target label (0-1). So a large negative value would become a logit close to 0 after applying the sigmoid function.
But if all values are negative this implies that all outputs would be close to 0, which could be wrong! And after 1000 Epochs, with this small dataset, your model will most likely overfit!
How big is your loss? Have you tried calculating the accuracy, because this could be a more easily interpretablemetric?